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cetfws, 


This is a report on work performed during Phase I of the 
MAPS program from 21 September 1971 to 20 March 1972. 
The contract was administered by NASA-Langley Research 
Center under the AAFE program and monitored by Dr. J. M. 
Russell, III. 

The principal investigator was C. B. Ludwig, the co-investi- 
gators were W. Malkmus, M. Griggs, and E. R. Bartle. 

The following personnel have contributed to various subjects: 
W. Marggraf, J. Wang, C. N. Abeyta and R. M. Smith. 
Fruitful discussions were held with L. L. Acton and G. D. 
Hall of Convair Aerospace . 




j a r 

2'X\ • 


m,£ ATW 



ABSTRACT 


The present study has provided results upon which the 
design of a satellite "Remote Gas Filter Correlation (RGFC)" 
instrument can be based. Although a final decision about 
the feasibility of measuring some of the pollutants with the 
required accuracy is still outstanding and subject to further 
theoretical and experimental verifications, viable concepts 
are presented which permit the initiation of the design phase. 
The pollutants which are of concern in the troposphere and 
stratosphere were selected. The infrared bands of these 
pollutants were identified, together with the bands of inter- 
fering gases and the line parameters of the pollutants as well 
as interfering gases were generated through a computer 
program. Radiative transfer calculations (line-by-line) 
were made to establish the radiation levels at the top of the 
atmosphere and the signal levels at the detector of the RGFC 
instrument. Based upon these results the channels for the 
RGFC were selected. Finally, the problem areas, which 
need further investigations, were delineated and the supporting 
data requirements were established. 


iv 



CONTENTS 


FOREWORD iii 

ABSTRACT iv 

ILLUSTRATIONS viii 

TABLES xiv 

LIST OF SYMBOLS xv 

SUMMARY xviii 

1 INTRODUCTION 1 

2 MODEL INPUT (Task 1) 4 

2.1 THEORY OF SYMMETRIC AND NEARLY SYMMETRIC 4 

TOP MOLECULES 

2.2 DESCRIPTION OF SYMTOP PROGRAM 7 

2.3 RESULTS AND COMPARISON WITH LABORATORY DATA 8 

2.3.1 Carbon Monoxide 8 

2.3.2 Carbon Dioxide 1 1 

2.3.3 Sulfur Dioxide 14 

2.3.4 Nitrogen Dioxide 15 

2.3.5 Nitric Oxide 17 

2.3.6 Ammonia 18 

2.3.7 Nitrous Oxide 25 

2.3.8 Formaldehyde 2 8 

2.3.9 Methane 28 

2.3.10 Water 30 

REFERENCES 30 

3 MODEL STUDIES FOR GASEOUS POLLUTANTS (Task 2) 32 

3.1 THEORY OF RADIATIVE TRANSFER 32 

3.2 DESCRIPTION OF GAS FILTER CORRELATION INSTRUMENT 36 

3.3 DESCRIPTION OF THE POLA YE R PROGRAM 38 

3.3.1 The Basic POLA YER Program 38 

3.3.2 Test Calculation with POLA YER for the Solar Spectrum 39 

at 4.65p, 

3.3.3 Versions of POLAYER 41 

3.3.4 Future Development of POLAYER 43 

3.4 RESULTS 44 


v 



3.4.1 Carbon Monoxide 44 

3. 4. 1. 1 Comparison with Band Model Calculations 46 

3.4. 1.2 Influence of Stratospheric CO 46 

3. 4. 1. 3 Influence of Water Vapor 46 

3. 4. 1.4 Influence of Stratospheric H 2 0 on CO 49 

Measurements 

3. 4. 1.5 Model Calculations for CO with H 2 0 and N 2 0 50 

3. 4. 1.6 Model Calculations for CO with H 2 0, N 2 0, 53 

C0 2 and O3 

3. 4. 1.7 Comparison Between DV and Radiance 55 

Measurements 

3. 4. 1.8 The Radiance as a Function of Sun Angle 55 

3. 4. 1.9 Av as a Function of Sun Angle, Ground E mi ssivity 60 

and Temperature 

3.4.1.10 Influence of Variable Ground Emissivity on AV 60 

3.4.2 Carbon Dioxide @7 

3. 4. 2.1 Model Calculations for C0 2 at 1.6M- 67 

3. 4. 2. 2 Model Calculations for C0 2 at 2.1 y 67 

3. 4. 3 Sulfur Dioxide 69 

3. 4. 3.1 Model Calculations for S0 2 at 8.6M- 69 

3. 4. 3. 2 Model Calculations for S0 2 at 4M- 69 

3.4.4 Nitrogen Dioxide 7 q 

3. 4.4.1 Model Calculations for N0 2 at 7.6 M- 70 

3. 4. 4. 2 Model Calculations for N0 2 at 3. 3M- 70 

3. 4. 5 Nitric Oxide 76 

3.4.6 Ammonia 71 

3. 4.6.1 Model Calculations for NHg at 10. 5M- 71 

3. 4. 6. 2 Model Calculations for NH3 at 3M- 71 

3. 4. 7 Formaldehyde 72 

REFERENCES 72 

4 MODEL STUDIES FOR AEROSOL DETECTION (Task 3) 73 

4.1 EARTH-ORIENTED OBSERVATIONS 73 



4.1.1 Radiance Observations 73 

4.1.2 Polarization Observations 76 

4. 1. 3 Comparison of Radiance and Polarization Observations 80 

4.1.4 Observations Over Land or Clouds 81 

4.1.5 Problem Areas 87 

4. 1.5.1 Sim Glitter 87 

4. 1.5. 2 Aerosol Properties 88 

4. 1. 5. 3 Non-Lambertian Surfaces 88 

4. 1.5. 4 Surface Reflectance Gradient 88 

4.2 SUN-ORIENTED OBSERVATIONS OF AEROSOLS 89 

4.2.1 Cloud Interference 89 

4.2.2 A tmo sphe ric A ttenuation 91 

REFERENCES 95 

5 INVESTIGATION OF RGFC INSTRUMENT PARAMETERS (Task 4) 97 

5.1 OPTIMIZATION PROCEDURE 97 

5.2 LABORATORY CALIBRATION PROCEDURE 103 

5.3 IN-FLIGHT CALIBRATION 115 

5.4 INVESTIGATION OF RATIO TECHNIQUES 116 

5.4.1 Ratio of Signal Difference of One-Cell Pair to Signal 116 

5.4.2 Ratio of Signal Difference of Two- Cell Pairs 126 

6 STUDY OF MISSION REQUIREMENTS (Task 5) 132 

6.1 ATMOSPHERIC TEMPERATURE AND PRESSURE PROFILE 132 

6.1.1 Temporal and Spatial Variations of Temperature 132 

Profiles in the Real Atmosphere 

6.1.2 The Accuracy of Temperature Data 140 

6.1.3 Summary of Requirements for Atmospheric 143 

Temperature and Pressure Profile Data 

6.2 WATER VAPOR PROFILE 150 

6.2.1 Temporal and Spatial Variations of Water Vapor Profiles 150 

in the Real Atmosphere 

6.2.2 Water Vapor from SIRS (NIMBUS IV) 154 

6.3 CLOUD COVER 154 

6.4 THEORETICAL SENSITIVITY OF GFC INSTRUMENT TO 157 

UNCERTAINTIES IN ATMOSPHERIC PARAMETERS 

vii 



ILLUSTRATIONS 


2-1 Monochromatic transmissivity of the RO and R1 lines of 9 

the first overtone of CO for u=9. 0 cm-atm, T=300°K and 
p=l atm. 

2-2 Degraded transmissivity (slit width of 0. 11 cm 1 of the 10 

RO and R1 lines of the first overtone of CO for u=9. 0 cm-atm, 

T=30CPK and p=l atm (calculated spectrum solid, laboratory 
spectrum^ 6 ). 

Monochromatic transmissivity of the band system at 1. 6(j, of C0 2 13 

between 6330 and 6336 cm ^ for u=800 cm-atm, p=0.5 atm and 
T=296 J K. 

Degraded transmissivity (slit width of 0. 08 cm -1 ) of spectrum 13 

from Figure 2-3 (calculated spectrum solid, laboratory 
spectrum^ dashed). 

2-5 Monochromatic transmissivity of the S0 2 -vi band between 1140 

and 1180 cm -1 (showing band center and R- branch maximum). 

The conditions are the same as in Figure 2-6. 

2-6 Degraded transmissivity (slit width is 1 cm -1 ) from Fig. 2-5 

for u= .577 cm atm p e =l. 59 atm and T-298 U K (calculated 
spectrum solid, laboratory spectrum dashed). 

2-7 Monochromatic transmissivity of R(5 ) and R(5 ) lines 19 

of fundamental band of NO between 1896 and 1898 cm -1 . 

The conditions are the same as in Figure 2-8. 

2-8 Degraded transmissivity (slit width is 0. 13 cm -1 from Figure 20 

2-7 for u=0.263 cm-atm, p=200 torr and T=300 U K (calculated 
spectrum solid, laboratory spectrum^ 21 ) dashed). 

2-9 Monochromatic transmissivity of portions of the 10. 5p, band 21 

system of NH 3 . A segment of the Q-branch between 900 and 950 
cm -1 is shown on the left hand side and a segment of the R-branch 
between 1150 and 1200 cm -1 on the right hand side. The con- 
ditions are the same as in Figure 2-10. 

2-10 Degraded transmissivity (slit width of 8 cm 1 ) of a portion of 23 

the 10. 5p, band system of NH 3 between 900 and 1260 cm -1 for 
u=l. 617 cm atmgrpp p e =4185 torr and T=299 U K (calculated 
spectrum solid, laboratory spectrum^ 23 ) dashed). 

viii 


16 


16 


-3 


2-4 



2-11 Monochromatic transmissivity of NH 3 between 3100 and 3700 24 

cm -1 for u=6.23 cm atmg-pp. P e = 5090 torr » 211(1 T=299PK. 

2-12 Degraded transmissivity (slit width of 50 cm 1 ) from Figure 24 

2-11. Calculated spectrum solid; laboratory spectrum^ 23 ) 
dashed. 

2-13 Monochromatic transmissivity of NgO between 2204 and 2206 27 

cm" 1 for u=0.31 cm atm, p=2 torr, and T=220 U K. 

2-14 Degraded transmissivity (slit width of 0.05 cm -1 ) of spectrum 27 

from Figure 2-13 (calculated spectrum solid, laboratory 
spectrum^ 27 ) Dashed). 

2- 15 Monochromatic transmissivity (calculated) of HCHO between 29 

2880 and 2900 cm -1 (solid curve). Degraded transmissivity 
(measured by Hanst( 30 ) for u=0. 08 cm atm, p=0. 0079 atm and 
T=300 p K (dashed curve). 

3- 1 Upwelling radiation at the top of atmosphere. Contributions 35 

from reflected sun radiation E^(u)), atmospheric scattering 
E (uj), and thermal radiation Eq(cd). Ground temperature = 

3 (HP K, sun zenith angle =45° , particulate matter according 
to Elterman 1968. 


3-2 Transmissivity of atmosphere between 2135 and 2165 cm" 1 . 40 

Upper trace calculated mono chromatically, lower trace calcu- 
lated with slit function (dashed), experimental curve solid line. 

3-3 Flow diagram of SYMTOP/POLAYER logic. 43 

3-4 Comparison of linear and quadratic interpolation for the line-by- 45 

line calculations. 

3-5 Temperature, pressure and water vapor profiles used in the 47 

model calculations. 

3-6 Signal difference vs. CO concentration in troposphere for 48 

different amounts of water vapor. 

3-7 Signal difference vs. CO concentration with and without normal 51 

water vapor, indicating "negative correlation". 

3-8 Water vapor distribution according to Gutnick. 52 

3-9 Signal difference vs. CO in the presence of N 2 0 and H 2 0. 54 

3-10 Normalized signal difference and radiance vs. CO concentration 56 

with and without water vapor . 

3-11 Radiance at the top of atmosphere vs. ground emissivity for 57 

different sun zenith angles at a ground temperature of 270K 
and CO concentration = 0.1 ppm. ^ 


ix 



3-12 

3-13 

3-14 

3-15 

3-16 

3-17 

3-18 

3- 19 

4- 1 

4-2 

4-3 

4-4 

4-5 

4-6 

4-7 

4-8 

4-9 

4- 10 

5- 1 


Radiance at the top of atmosphere vs. ground emissivity 
for different sun zenith angles at a ground temperature of 
300K and CO concentration = 0. 1 ppm. 

Radiance at the top of atmosphere vs. ground emissivity for 
different sun zenith angles at a ground temperature of 330K 
and CO concentration = 0.1 ppm. 

Radiance vs. ground temperature for e = 1. 0 and e = 0. 8 for 
various sun zenith angles. 

Signal difference vs. CO concentration for different sun zenith 
angles and e = 0. 8 and e = 1. 0 at ground temperature = 270K. 

Signal difference vs. CO concentration for different sun zenith 
angles and e = 0. 8 and e = 1. 0 at ground temperature = 300K. 

Signal difference vs. CO concentration for different sun zenith 
angles and e = 0. 8 and e = 1. 0 at ground temperature = 330K. 

Signal difference vs. ground temperature at e =1. 0 and s = 0. 8 
for various sun angles and for CO concentration = 0.1 ppm. 

Signal difference vs. CO£ concentration for sun zenith angle = 
30° , ground reflectivity = 0. 02, ground temperature = 300K. 

Normalized upward radiance vs sun angle as function of 
wavelength. 

Normalized upward radiance vs aerosol mass loading as 
function of sun angle. 

Polarization vs aerosol content as function of nadir look angle. 

Reflected radiance vs albedo for different aerosol contents. 

Reflected radiance vs aerosol content as function of albedo. 

Polarization vs albedo for different aerosol contents. 

Polarization vs aerosol content as function of albedo. 

Schematic of occultation experiment. 

Transmission vs aerosol content as function of wavelength 
(grazing altitude 10 km). 

Transmission vs aerosol content as function of wavelength 
(grazing altitude 15 km). 

Normalized filter function used in the model calculations for CO, 
corresponding to OCLI filter 14-8806-760. 


58 


59 


61 

62 

63 

64 

65 

68 

75 

77 

79 

82 

84 

85 

86 
90 

93 

94 
99 


x 



5-3 

5-4 

5-5 

5-6 

5-7 

5-8 

5-9 

5-10 

5-11 

5-12 


Optimization of instrument response as a function of CO con- 100 

centration in instrument cell. Instrument parameters i = 1 cm, 

PT c 1 atm > t G = 288 * 16 ° K > c (w) = c '( cjd ) , given by Fig. 5-1. 
Atmospheric parameters: CO concentration in the troposphere = 

0. 1 ppm, ground temperature = 288. 16PK, lapse rate = -6.5°C/km. 
Arrow points to maximum response, i.e., CO concentrations = 

0.35, corresponding to a cell transmission of 0. 83. 

Sensitivity of instrument as a function of CO concentration in the 101 
troposphere with three CO concentrations in instrument cell as 
a parameter: 0. 05, 0. 35 and 1. 00. Other instrument and atmos- 
pheric parameters same as in Fig. 5-2. 

Sensitivity of instrument as a function of total pressure in instru- 104 
ment cell with nine CO concentrations in the troposphere. The 
figures on the curves refer to the multiplier, with which the con- 
centration of 0.1 ppm in the troposphere is multiplied. Arrows 
refer to maximum sensitivity. 

Sensitivity of instrument as a function of CO concentration in the 107 
troposphere for two different slit functions. The lower curve is 
for the slit function given in Fig. 5-1, the upper curve for a square 
slit function. If the upper curve is corrected to account for the 
larger area under the slit function, the points marked as (x) are 


obtained. 

Sensitivity of instrument as a function of optical thickness in the 109 
instrument "sample" (calibration) cell. 

Three normalized filter functions used in the computer Simula- 111 
tion studies. 

Sensitivity of the instrument to different ground temperatures 117 

and sun zenith angles as a function of ground emissivity at 
4.6p, for zero concentration of CO in the troposphere. 

Sensitivity of the instrument for the same conditions as in Fig. 118 

5-8, but referenced to DV for 0.1 ppm CO concentration. 

DV and DV/V as a function of CO concentration in the tropos- 119 

phere with and without water vapor. 

DV/V as a function of CO concentration in the troposphere for 121 

different ground temperatures and variations of atmospheric 
temperatures and water vapor concentration. 

DV/V as a function of ground temperature for various CO con- 122 

centrations in the atmosphere. 


xi 



5-13 

DV as a function of V for different conditions of CO 
concentration, surface temperature, atmospheric tem- 
perature and water profile perturbations, surface emissi- 
ties, and sun zenith angles. 

123 

5-14 

Perturbation of DV and V due to changes in the atmospheric 
temperature profile by ± 5°K and in the water vapor profile by 
± 50% in mixing ratio. 

125 

5-15 

Sensitivity of instrument as a function of total pressure in the 
instrument cell for different concentrations of CO in the 
troposphere, having a profile according to B37 (see Table 5-6). 

127 

5-16 

Sensitivity of instrument as a function of total pressure in the 
instrument cell for different concentrations of CO in the trop- 
osphere, having a profile according to B36 (see Table 5-6). 

128 

5-17 

Sensitivity of the instrument as a function of the total optical 
thickness of CO for the three profiles B20, B37 and B36, using 
cell pressures of 1 and 0.4 atm. 

130 

5-18 

The ratio of DV as a function of total optical thickness for 
three CO concentration profiles (see Fig. 5-17, using two cell 
pairs with pressures of 1 and 0.4 atm. 

131 

6-1 

Temperature and temperature gradient at 500 mb level for 
February 1971 in San Diego, California at 0000 and 1200 Z. 

133 

6-2 

Temperature and temperature gradient at 500 mb level for 
February 1971 in Washington, D. C. at 0000 and 1200 Z. 

134 

6-3 

Temperature and temperature gradient at 500 mb level for 
July 1970 in San Diego, California at 0000 and 1200 Z. 

135 

6-4 

Constant pressure chart at 500 mb level for February 1, 1971 
at 1200 Z. 

138 

6-5 

Temperature spatial distribution for February 1, 1971. 

139 

6-6 

Variability of relative humidity, temperature at 700 mb and 
thickness of 700-500 mb layer for 48 hours. 

151 

6-7 

Variability of temperature and water vapor mixing ratio for 
selected days in February 1971 in San Diego, California. 

152 

6-8 

Signal difference vs. fractional cloud cover for different 
amounts of CO concentration. Cloud altitude = 1 km. 

155 

6-9 

Signal difference vs. cloud altitude for 100% cloud cover 
for different amounts of CO concentration. 

156 


xii 



6-10 


158 


Ratios of DV/V vs. CO concentration for A (ground 
temp. = 300K), B (0.6 cloud at 1 km, T =291K), C (ground 
temp. = 285, cloud at 1 km, T = 284. 9K), D (cloud at 2 km, 

T = 278. 4K), E (cloud at 3 km, T = 271. 8K), F (ground temp. 

= 270K, cloud at 4 km, T = 265. 4K). 

6~ 11 Details of pressure and temperature profile in the troposphere. 159 

6-12 Signal difference vs. CO concentration for ground temperatures 160 

of 300 and 285K and atmospheric temperature variations. 

6-13 Radiance vs. CO concentration for same conditions as Figure 6-12. 161 

6- 14 Signal difference vs. CO concentration for different water vapor 162 

concentrations and different ground temperatures. 

6-15 Signal difference change per degree change per layer vs. altitude 164 

for no water and with water (Gutniek) for CO concentration of 
0.2 ppm and ground temp. = 300K. 


xiii 



TABLES 


1 Gases which are of Concern in the Troposphere and 
Stratosphere 

2 Infrared Bands of Pollutants and Interfering Gases 

3 Line Parameters Generated by SYMTOP for Pollutants 
and Interfering Gases 

4 Model Calculations Performed with PO LAYER 

5 Channels Selected for the GFC Instrument 

6 Problem Areas which Need Further Investigation 

7 Supporting Data Requirements 

2- 1 Types of Molecules Considered in Atmosphere 

3- 1 Difference of Signal Changes Due to CO with and without Water 

Vapor in the Stratosphere 

3-2 Signal Change for CO in the Presence of the Interfering 

Molecules H O and N O 
u Z 

3-3 Signal Change for CO for Varying and Constant Ground 
Emissivity 

3-4 Signal Change for NO as a Function of Concentration 
5-1 Effect of Instrument Cell Temperature Variations 
5-2 Effects of Lab Calibration Applied to Real Atmosphere 

5-3 Effects of Lab Calibration Applied to Real Atmosphere 

5-4 Effects of Lab Calibration Applied to Real Atmosphere 

5-5 Signal Change of CO for the Calibration Temperatures 

5- 6 CO Concentration for 3 Distributions 

6- 1 Verification of NMC Operational Temperature Analysis 


xviii 

xx 

xxii 

vii 

xxvi 

xxvii 
xxix 

5 

50 

53 

66 

71 

102 

110 

113 

114 

115 
126 
145 



LIST OF SYMBOLS 


A 
B 

B , 

C 

Ch(0) 


N 


C(<D) 

C'((1)) 


c 

E 

E 

E 

E 


G 

R 


E 


R0 


F 

V J »K) 

f(0) 

H 

s 

h 

J 

K 


Area of entrance pupil 
Rotational constant 
Rotational constant 

Ratio of self -broadening to nitrogen broadening efficiency 
Rotational constant 

Chapman function describing the path of solar radiation 
through the atmosphere for sun zenith angles 0 greater 
than 60° . 

Generalized instrument response function of the GFC 
channel containing gas cell. 

Generalized instrument response function of the 
GFC channel containing reference cell. 

Concentration 

Upwelling radiance 

Thermal radiation of ground and atmosphere 

Solar radiation reflected from ground 

Solar radiation scattered by the atmosphere without 
having been reflected by the ground. 

Solar radiation reflected from ground outside 
the field of view and scattered by the atmosphere 
into the field of view. 

Lower energy level for a given transition 

Rotational energy levels 

Equals sec 0 for 0 £ 0 ^ 60° and Ch(0) for 0 > 60° 

Sim irradiance at the top of the atmosphere 
Altitude of satellite 
Rotational quantum number 
Rotational quantum number 


xv 



LIST OF SYMBOLS 


k(oj) Monochromatic absorption coefficient in cm"" 1 atm"" 1 
i Instrument cell length 

n Number of scatterers 

N° Blackbody radiance 

p Partial pressure of pollutant gas 

P N Partial pressure of nitrogen gas 

P T Total pressure 

P e Equivalent pressure, defined by Eq. (5-5) 

Ho Instrument responsivity, wave number independent 

S Line strength 

T Temperature 

T^, Temperature of calibration source 

T g Temperature of ground 

u Optical depth in cm atm 

V Signal voltage 

AV, DV Signal voltage difference of gas and reference channels 
z Altitude 

a Line half width 


Y Function of sun zenith angle = 1 + f (0) 

e Emissivity 

9 Sim zenith angle 

h Asymmetry parameter 

Extinction coefficient for species i 
A. Wavelength 

-4 

U Micrometer = 10 cm 

v Symbol for transition 

o - Scattering cross section 

xvi 



LIST OF SYMBOLS 


Transmissivity through atmosphere = t t 

o p 

Transmissivity through atmosphere due to normal 
constituents only. 

Transmissivity through atmosphere due to pollutant 
species only. 

Transmissivity through gas cell of GFC instrument 
Transmissivity through reference cell of GFC instrument 
Solid angle in steradians 
Wavenumber in cm * 


xv ii 



SUMMARY 


The results of the first phase of the MAPS program entitled 1 'Experiment Definition 
Phase" are summarized in seven Tables. 

Table 1: Gases which are of concern in the troposphere and stratosphere . 

The gases given in this table have been selected on the basis of our work con- 
ducted for NASA and reported in CR-1380 (July 1969) and in its sequel, "Air 
Pollution Detection by Satellites" (to be published soon) and also on the basis of 
the study conducted by the Norfolk Working Group whose results were reported in the 


Table 1. Gases which are of Concern in the 
Troposphere and Stratosphere 



Troposphere 

Stratosphere 

C -compounds (Inorganic) 

co 2 ,co 

co 2 ,co 

S -compounds 

SO .H S 

SO 


2 2 

2 

N-compounds 

NO, NO , NH 

NO, NO . N 0,HN0 


2 3 

NHg 2’ 2 ’ 3 

HC -compounds 

HCHO,PAN,CH 4 

CH 4 , < HC> 

Oxidants 

0„ 

o 


3 

3 

Halogens 

I , HF,Br 


Water Vapor 


H 2° 

Particulates 

Total Aerosol 

Layered Aerosol 


Content 

Content 
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NASA publication SP-285 (1971). The regions of the atmosphere in which the 
pollutants are of concern are the troposphere and stratosphere. Two distinct 
observational modes apply to these two regions. For the troposphere, only a nadir 
experiment appears to be practical, while for the stratosphere, only an occultation 
experiment appears to be applicable. Although the task description of Phase I of 
NAS1-11111 included the pollutant modeling of both the troposphere and stratosphere, 
time and funding limitations have made it impossible for us to include the stratospheric 
modeling. In addition, the sensor and mission requirements for a sun occultation 
experiment are so different that a separate development is warranted. (This subject 
will be pursued further in Phase II) . 

The following discussions are limited to the nadir observations. Of the 14 

gases listed for the troposphere, the halogens, O , PAN and H S have not been 

3 2 

included in our model studies, because the concentrations and/or the band strengths 
of these pollutants (except O ) are so low as to make a remote measurement from 

O 

satellites extremely difficult. To observe O in the troposphere in a nadir experi- 

O 

ment is also nearly impossible because of the stratospheric ozone layer. 


Table 2: List of infrared bands of pollutants and interfering gases. 

Seven gaseous pollutants of importance in the trosposphere were investigated. 

The active infrared rotation vibration bands of suitable intensity were selected and 

the transitions of the interfering species were identified and listed in matrix 

form. It is readily apparent that the normal atmospheric constituents CO , 

2 

N 2 0, CH^ and H^O are present in almost every spectral region where pollutants are 
active. In addition, some isotopes and ozone are also interfering. The solar lines 
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lable 2. Infrared Bands of Pollutants and Interfering Gases 









































































are of importance only at the shorter wavelengths. Their line position and strength 
are rather well known from the solar tables and they wall be considered in the future 
data interpretation procedures. In the present modeling studies, they were not 
included. 

Table 3 : List of line parameters generated for pollutants and interfering gases . 

All the transitions listed for the pollutants and interfering species in Table 2 are 
considered for the model input. The status of generating the line parameters is 
given in this table and the qualitative agreement with available experimental data is 
indicated in the last column. The line parameters of nearly 30, 000 lines of nine 
molecules have been generated by the SYMTOP program. 


Table 4: List of model calculations performed. 

The line parameters of the molecules listed in Table 3 were used as input to the 
POLAYER program. The POLAYER program calculates the upwelling radiation at 
the top of the atmosphere, including the sun reflected radiance, the thermal emission 
of the earth's surface and the thermal emission of the atmosphere. In addition, the 
POLAYER program calculates the signal difference at the detector after the radiance has 
passed through both the vacuum and specifying cells. The radiative transfer calculations 
through a non-homogeneous medium were made with the following assumptions: 

1. The atmosphere i s plane -parallel, with 8 layers in the troposphere 
and 12 layers in the stratosphere; 

2. Atmospheric refraction can be neglected; 

3. The only consideration of particulates is in the transmission losses; 

4. The earth's surface is a "diffuse" reflector; 

5. The absorbing gases are in local thermodynamic equilibrium. 
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Table 3- Line Parameters Generated by SYMTOP for 
Pollutants and Interfering Gases 


Transition # Lines 


. 3 V y 3 

2 V V 3 

3V 1 _\j l +v 
2 2 3 

2v +v 
2 3 


2v -v 
1 2 


Benedict et 
SYMTOP 
SYMTOP 
SYMTOP 


SYMTOP 


SYMTOP 


Comparison 


Excellent 

Excellent 

Excellent 

Excellent 


Excellent 


Excellent 


v - 2v 
3 2 


SYMTOP 


Excellent 


SYMTOP 

SYMTOP 


SYMTOP 

SYMTOP 

SYMTOP 


Not Available 
Not Available 
Excellent 


2v +v 
2 3 


SYMTOP 


Excellent 


v +2v 
1 _ 2 


SYMTOP 


SYMTOP 

SYMTOP 

SYMTOP 


Excellent 


Excellent 

Excellent 

Excellent 


SYMTOP 


Excellent 


SYMTOP 


SYMTOP 


Very Good 


SYMTOP 


SYMTOP Very Good 

Gates et al. Excellent 
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Table b. Model Calculations Performed with POLA.YER 


1 

X 

(M) 

T 

Ground 

(°K) 

* Ground 

0 

Atm. 

Temp. 

Prof. 

Interf. 

Species 

Incl. 

Interf. 
Species 
Not Incl. 

PoU. 

Cone. 

ppm 

Poll. 

Cone. 

Prof. 

H 2° 
Cone. & 

Prof. 

Time 

per 

Run(min) 

CO 

4.6 

270 - 330 

0. 0-1.0 

0° -90° 

Sev- 

eral 

H 2° 

n 2 o 

C °2 

°3 

0-0.4 

3 

4 

16 

C °2 

■ 

300 

.98 

30° 

1 


H 2° 

CH 4 

320 

1 


50 

m 

300 

.98 

30° 

1 

H 2° 


300- 

400 

n 


16 

S °2 

— 

8.6 

300 

1.00 

0° 

■ 

H 2° 

n 2 ° 

°3 

.002 

i 

1 

75 

4.0 

300 

.8 

30° 

1 


n 2 ° 

CH 4 

.002 

i 


45 

no 2 

■ 

300 

1.0 

0° 

■ 

H 2° 

n 2 ° 

cir 

4 

.002 

■ 


61 

3.3 

300 

.8 

30° 

1 


CH 4 

H 2° 

.002 

■ 


40 

NO 

5.3 

300 

.8 

30° 

■ 

H 2° 

C °2 



■ 


23 

**3 

10.5 

300 

1.00 

0° 

■ 


H 2° 

o„,co„ 

3 2 

.02 

i 


15 

3 

300 

.8 

30° 

1 


c° 2> n 2 ° 

CH 4> H 2° 

.02 

i 


15 

HCHO 

3.5 

300 

.8 

30° 

1 


H 2° 

.002 

i 
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The table shows a list of conditions for which the model calculations were made. 
More extensive model calculations were made for CO than for the other gases. 
Nevertheless, the model calculations we did make for the other molecules were 
sufficient to establish the feasibility of their measurement and to determine the 
instrument cell parameters. Only one of the pollutants (NO) was found not feasible 
to be remotely measured by RGFC. The others were found to be either ’firmly" 
or "conditionally" feasible (definition given in next paragraph). It is essential that 
more model calculations be performed in order to ascertain fi rml y the ability of 
the RGFC instrument to measure all pollutant concentrations within given error 
limits. However, it is necessary to select the variation of the parameters judiciously 
in order to limit the model calculations to a reasonable number. Although we have 
incorporated many short cuts into the calculational routines (simultaneous radiance 
calculations for all concentrations of a given pollutant, establishment of linear 
relationships, etc.), the number of model runs for each pollutant is still very large. 
Assuming it is desired to calculate the signal levels for three ground temperatures, 
three ground emissivities, three sun angles, three atmosphere temperature profiles, 
three pollutant concentration profiles and four water vapor distributions, tiie total 
number of model runs amounts to 3 x 4 = 972 for each of the nine spectral intervals. 
Since some of the calculations involve several thousand spectral lines, and require 
computing times of up to an hour, it has become mandatory to consider band modeling. 
The band modeling will be based upon the line-by-line parameters and can be checked 
for accuracy. Thus, we will have a higher confidence level in these calculations 
than we would have without these checks with the line-by-line calculations. 



Table 5; List of channels selected for the GFC instrument. 

Based on the results of the model calculations, decisions about the gases to be included, 
number of channels, spectral regions and cell parameters have been made. Because 
the model calculations were not sufficient for some of the pollutants, the decisions 
were not 'firm" in all cases. This is indicated in the third column through the label 
"Cond", which means that the decision about the inclusion of the channel pair(s) is 

j|( 

"conditional", subject to the outcome of further model calculations (see next paragraph). 

In order to proceed with the design phase without the benefit of further model 
calculations in the immediate future, we assume that the feasibility and utility of 
channel pairs in question can be "firmly" established through later model calculations. 
Therefore, the design concept will be based on 13 channel pairs with the cell lengths 
established at 1 cm, 2 channel pairs with the cell lengths established at 20 cm for the 
CC> 2 measurement and two radiometer channels (indicated as one channel pair in 
Table 5), for the determination of the aerosol content in the atmosphere. 

In case the planned model calculations should indicate that it is not practical 
to include one or more of these channel pairs in question, the design of the RGFC 
may be changed to eliminate the "superfluous" channel pairs or to retain them for 
use as additional channels for the pollutants established as "firm". 


♦Although no model calculations for methane have been made yet, the measurement 
of this molecule is believed to be feasible. Recent studies by B. Weinstock and 
H. Niki (Science 176, 290, 1972) indicate that this molecule plays an important 
role in the balance of carbon monoxide. 

/. 
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Table 5. Channels Selected for the GFC Instrument 


Poll 

X 

pm 

Level 

of 

Conf 

No. of 
Channel 
Pairs 

J l 
cm 

CELL PARAMETERS 
c P>j> u 

atm cm atm 

C °2 

ED 

Cond 

2 


1 

2.5 

50 

CO 

4.6 

Firm 

2 

1 

.35 

1 

.35 

so 2 

8.6 

Firm 

2 

1 

.1 

1 

.1 


4 

Cond 

1 

1 

.2 

1 

.2 

no 2 

7.6 

Cond 

2 

1 

.1 

1 

.1 


3.3 

Cond 

1 

1 

.2 

1 

.2 

nh 3 

10.5 

Firm 

2 

1 

.2 

1 

.2 


3 

Firm 

1 

1 

.25 

1 

.25 

HCHO 

3.5 

Cond 

2 

1 

.2 

1 

.2 

CH 4 

3.3 


2 

1 

.4 

1 

.4 


7. 7 

p. XXV. 

1 

1 

.3 

1 

.3 

Aerosol 

Visi- 

ble 

Firm 

1 

— 

— 

— 

— 


Table 6; Problem areas which need further Investigations. 

This table lists the problem areas which became apparent for certain pollutants 
through the model calculations and which placed these pollutants in the category 
Cond' of Table 5. We believe that all of these problem areas can be satisfactorily 
solved through additional theoretical and experimental studies, conducted both in 
Phase II and in an independent program proposed for the ERTS-A Data Analysis. 


xxvi 





























TTAXX 


Table 6. Problem Areas which Need Further Investigation 




Statement of the Problem 

Indicated Solution 

o 

o 

to 

2.1 

1. The unknown dependency of the CC> 2 
transmissivity on the Mie scatter- 
ing source term 

2. The influence of the variation in 
ground reflectivity and atmospheric 
scattering on high accuracy reqmts 

1. Perform Monte Carlo calculations of 
multiple Mie scattering with the inclusion 
of C0 2 absorption 

2. Conduct model calculations with two 
cell pairs to minimize ground and 
atmospheric effects 

S °2 

4 

Present calculations indicate negative 
signals, even with no interfering gases, 
apparently due to difference in X depend- 
ency between calibration blaclcbody and 
sun reflected radiance 

Establish firmly the cause for the negative 
signal through additional computer model- 
ing: investigate simulation of the X depend- 
ency of sun reflected radiance during calib- 
ration. 

N °2 

7.6 

Present calculations indicate strong 
influence of water vapor absorption 

Determine whether measurement uncertainty 
is within accuracy requirements for concen- 
tration determination 


3.3 

See under SO - 4^m 
^ 

See under SO - 4M-m 
2 

HCHO 

3.5 

See under SO - 4(^m 

See under SO - 4nm 

Aerosol 

Visible 

1. Sun glitter from ocean surface 

2. Aerosol Properties 

3. Non- Lambertian land surface 

4. Surface reflectance gradient 

1 . Establish utility of two wavelengths 
measurements 

2. Establish theoretical and/or empirical 
relationships through a field test pro- 
gram (ERTS-A) 

3. Assumption of small influence to be 
checked in proposed field test program 

4. See under 3 























This latter program would provide invaluable data, because theoretical results will 
be compared with field measurements taken at the ground, from aircraft and 
satellite. 

The problem areas identified for the measurement of CO arise mainly from 

<u 

the fact that the requirement of at least 1% accuracy in the determination of the con- 
centration is much higher than for the other pollutants. More details are given in 
Section 3.4.2 of this report. 

Table 7: Supporting data requirements. 

This table gives a preliminary list of supporting data requirements. Most of these 
requirements are within the present state-of-the-art, such as ground brightness 
temperature, atmospheric temperature and pressure profiles, and cloud cover 
amount. Other supporting data such as water vapor profile and cloud top height are 
more uncertain at present. 

The most significant result of the study conducted under Task 5 of Phase I 
is the conclusion that the accuracy of the temperature profiles available from existing 
objective analyses by the National Meteorological Center (NMC) is sufficient for 
satisfactory pollutant data analysis, thus eliminating the need for a separate sounder 
on the same spacecraft. Unfortunately, the same cannot be said about the measure- 
ment of the water vapor distribution. Only low altitude distributions (up to the 500 mb 
level) are available from the objective analysis of the observations (NMC) at synoptic 
times (0000Z and 1200Z time). Current weather forecasts do not include relative 
humidity, nor are observations made at altitudes higher than the 500 mb level. The 
planned synchronous meteorological satellites (SMS) are scheduled to provide water 
vapor data at the higher altitude. 

The other constituents of the atmosphere, whose concentrations must be known 
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Table 7. Supporting Data Requirements 


Parameter 

Accuracy 

Requirement 

Provided on Same 
Satellite Platform 

Provided from 
Different Sensor 
Platform 

Potential Technique 

Ground brightness 
temp, at \ 

± 2°K 

One GFC channel for 
each spectral inter- 
val. Accuracy better 
than ± 2°K 



Atmospheric Temp. 
Profile 

± 2°K 


Objective analyses: 12 
hour forecasts from NMC 
within an accuracy of ± 2° K 
(RMS) 


Water Vapor 
Profile 

Not fully 
established 


Objective analyses: Mean RH 
in 1000-500 mb layer within 
an accuracy of ± 50% 

Profile measurements 
extending to higher 
layers thru SMS-(OES) 

Cloud Cover 
Amount 

±10% 

Vidicon, having the 
same fov 



Cloud Top 
Height 

Not yet 
established 

Estimate through 
temp, measurement 


o 

7615 A oxygen line 
determination 








for the pollutant data analysis, do not need to be measured. The assumption of the 
"standard atmosphere" is sufficient. 

The knowledge of cloud top height is desirable when pollutant data are to be 
obtained above the cloud cover. The absolute temperature of the cloud top, and thus 
its height, may be inferred from a measurement of the brightness temperature and 
an assumed emissivity. This subject needs further study, together with the utility 
of obtaining pollutant data in the presence of partial cloud cover. 
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INTRODUCTION 

Air pollution has long ceased to be only a local problem. The millions of tons of 

pollutants emitted into the atmosphere every year by industrial nations affect the 

entire atmosphere beyond the boundaries of local districts, states, and countries. 

The best-known examples of global pollution are the increase of carbon dioxide and 

of particulate matter, both of which have an influence on the global heat balance. 

But there are other lesser known phenomena, such as the formation of particulates 

from gaseous pollutants and the nature of removal mechanisms for some of the 

pollutants, which must be studied on a global basis. 

Satellites can provide the global coverage. Instruments on-board unmanned 

or manned spacecraft can measure the vertical total burden of pollutants above an 

area on the surface. When these data are continuously obtained and augmented with 

data from monitoring stations located on or near the surface of the earth, the dynamics 

of contaminated air masses, the distribution and dissipation of air pollutants and the 

long-rangp effect on the global climatology can be studied. 

Satellites can also provide the platform for sensing the upper atmosphere on 

a global scale which cannot be done as economically from other platforms, such as 

balloons. A knowledge of the chemical composition and reaction schemes of that 

atmospheric region is of great importance. Based on previous studies and the 

study by a Working Group of about 50 leading U. S. authorities in the field of remote 
(4) 

sensing of pollution , the pollutants which are of concern on a global basis are: 

CO The global mean concentration of CO is increasing and it has an effect 

2 » & 

on the global climate; 

CO The sources and sinks are not all known and the long range trends in 

global loading are not established; 
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S0 2 Accumulation has toxic effects; in addition it contributes to particle 
formation, thus influencing the global heat balance; 

H S Contributes to the formation of SO n through oxidation; 

NO x Accumulation has toxic effects and is a vital ingredient in the formation 
of photochemical smog. It also contributes to the formation of particles; 
nitrous oxide may be important in the chemistry of the stratosphere; 

NHg Contributes to the formation of particles; 

<HC > Are ingredients of photochemical smog and contribute to particle 
formation; CH 4 may contribute to the formation and decomposition 
of CO; 


Probably influences the chemistry and particle formation in the stratosphere; 

Changes in the stratospheric ozone layer have great effects in the amount 
of uv radiation reaching the earth's surface; ozone in the troposphere is an 
ingredient in photochemical smog; 


H 2° 

Others 


Is important in the chemistry of the stratosphere and particle size 
distribution; 

Fluorocarbons and halogens are toxic and are probably accumulating in 
the troposphere. 


The measurement of some of the above pollutants by a satellite -borne infrared 
RGFC instrument can be ruled out on the basis of low concentrations and/or weak 
or non-existing infrared rotation-vibration bands. The pollutants which do not fall 
into these categories and which are the subject of the present study are CO , CO, CH., 
S °2’ N0 2’ m s> HCH0 and particulates. These pollutants are of concern in the 
troposphere and the measurements are made with nadir experiments. The pollutants 
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which are of concern in the stratosphere are O^, HNO^, NO, NO^, N^O, CH^, SO^, NHg, 
H O and particulates. Their observation is made through limb or sun occultation 

A 

experiments. Although line parameters for most of these molecules have been 
generated in the present program because they were needed also for the nadir 
experiment, model calculations for the sun occultation experiment have not been 
performed as yet. 


REFERENCES 

1. C. B. Ludwig, R. Bartle and M. Griggs, NASA-CR-1380 (1969). 

2. Study of Critical Environmental Problems (SCEP) MIT Press (1970). 

3. Study of Man's Impact on Climate (SMIC) MIT Press (1971). 

4. Remote Measurement of Pollution, NASA SP-285 (1971). 
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2 

MODEL INPUT 
(Task I) 


2. 1 THEORY OF SYMMETRIC AND NEARLY SYMMETRIC TOP MOLECULES 


The molecules of interest to the present program are either symmetric top or nearly 
symmetric top molecules (with the exception of H O). A symmetric top molecule 
(i.e. , a molecule for which at least two of the principal moments of inertia are 
equal) may be treated analytically much more readily than an asymmetric top 
molecule. An asymmetry parameter h may be defined^: 


2B - (A+C) 
A-C 


(2-1) 


where the rotational constants A, B, and C are inversely proportional to the 
principal moments of inertia. The parameter h equals -1 for a prolate symmetric 
top, 0 for a most asymmetric top, and +1 for an oblate symmetric top. A list of 
the types of molecules considered is given in Table 2-1. 

The energy levels for a rigid symmetric top molecule are given by 

F y (J, K) = u> v + BJ(J+1) + (A-B) K 2 (2-2) 

where J and K are the rotational quantum numbers and v represents the set of 
vibrational quantum numbers. For a non-rigid top, (1) can be written more 

p 

precisely as a multiple power series in J(J+1) and K , with coefficients different for 
each vibrational state. For a transition between any two states, the wavenumber of 
the transition can be expressed simply as the difference between the energy levels; 


UJ (v / , J / ,K / - v, J,K) = F v ,(j',K') - F y (J,K) 


(2-3) 


4 



Table 2-1 Types of Molecules Considered in Atmosphere 


Linear Molecules 

CO 

NO 

C °2 

N O 
2 

Spherical Top Molecules 
CH 4 

Other Symmetric Top Molecules 
NH 

3 

Nearly Symmetric Top Molecules 
NO H = 9939 

Ct 

SO H = -. 9416 

Ct 

H 2 CO h = -.9623 

Asymmetric Top Molecules 
H 2 0 h = -. 4377 


Selection rules for J and K are obeyed: for a perpendicular band (change of 
electric dipole moment perpendicular to the top axis, e. g. , band of NH^ at 6.1 Mm), 

AJ = -1,0, + 1 ; AK = -1, + 1 . (2-4) 

For a parallel band (change of electric dipole moment parallel to the top axis, 
e. g. , v band of NH at 10. 5 Mm), 

& o 

AJ= -1,0, + 1 ; AK = 0 . (2-5) 
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By starting with J = 0, K = 0, the wavenumbers of all rotational transitions may be 
calculated. Relative strengths of the lines are known from quantum mechanics: 
e. g. , in a perpendicular band, when AJ = +1 and AK = +1, the line strength is pro- 
portional to^ 


. = (J+2+K) (J+l+Kl 

KJ (J+l) (2 J+l) ( 2 ~ 6 ) 

and the Boltzmann factor and statistical weight of the lower state. 

Similar expressions exist for all combinations of AJ = -l, o, +1 and AK = -1, +1. 
For a parallel band, when AJ = +1 ( and AK = 0), the analogous expression is 

A = i - J+1 > 2 - K 2 
KJ (J+l) (2 J+l) 

Similar expressions exist for AJ = -1,0 and AK = 0. 

Thus it is possible to generate a systematic list of transitions and relative line 
strengths for a symmetric top molecule. Certain molecules such as NH^ and CH 

having particular symmetry properties, require special modifications to calculate the 
statistical weights. 

Line parameters are outputted in a systematic pattern relative to the rotational 
quantum numbers, but are not sequenced by wavenumber. A separate sorting sub- 
routine can be used to select the lines lying in specified wavenumber intervals and sort 
these by wavenumber for input to the POLAYER program. 

Certain molecules (e.g. , SO^ and NO^) are not symmetric tops, but are very 
nearly so. The first-order effect of asymmetry on the spectrum is to change the 
wavenumbers of the transitions slightly. In such a case, SYMTOP calculates the 
lme parameters for the equivalent symmetric top, and makes first-order corrections 
for asymmetry to the wavenumbers of the transitions. 
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2.2 DESCRIPTION OF SYMTOP PROGRAM 

Since many of the species of interest are symmetric top molecules or slightly 
asymmetric tops, we have developed a program which calculates the energy level 
structure of a symmetric top molecule, and from the selection rules and matrix 
elements for transitions between these energy levels, calculates the wavenumbers and 
(relative) strengths of the allowed vibration-rotation transitions. A value of the total 
band strength is required to produce dimensionalized line strengths. As the transi- 
tions are calculated, the spectral data are written on magnetic tape. This tape may 
be used directly as input for program PO LAYER (to be described later) or it may be 
desirable to perform a presorting of the tape by wavenumber, possibly also merging 
one or more other tapes containing data for other molecular species. 

In several instances, existing tabulations of line parameters were available 
from the literature. This was the case for three H O band systems, as well as for 
the CO fundamental band. The spectral data were transcribed to punched cards, 
which were acceptable as input data to the earlier forms of POLAYER or as input 
to the presorting routines used for the later versions of POLAYER. 

The general form of SYMTOP which was developed required modifications for 
almost every different molecule to which it was applied. In general, intensity 

alternation results from statistical weight factors which depend on the symmetry 

16 

of the molecule and the spins of identical nuclei. CO^ , for example, has missing 

12 13 

alternate levels (i.e. , the statistical weight factor = 0) for both the C and C 

16 18 

species. However, CO O lacks this symmetry, and does not have such intensity 
alteration. In the case of NH , levels with K=3N +1 and 3N + 2 have statistical weight 

O 

of 1, while it is 2 for K=3N. 
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2. 3 RESULTS AND COMPARISON WITH LABORATORY DATA 
2.3.1 Carbon Monoxide 

CO-4. 6M - Line positions (cu) , strength (S), half-width (a) and lower energy level 

(F) of 145 lines for the fundamental band of C 12 0 and C 13 O have been taken from 
( 2 ) 

Benedict et al. Since the parameters S and ot were derived from experimental 
data, the accuracy is given by the experimental uncertainty which is estimated as 5%. 
Good agreement was also obtained between our calculated and a measured sun 
spectrum (see Section 3. 3). 

CO - 2.34 Unfortunately no listing similar to that for the fundamental is available 
for the first overtone band (0 -2). A SYMTOP program was used to generate line 
locations, line strengths, and lower state energy levels for this band. The molecular 
constants were taken from the experimental data of Benedict et al. Data for 160 
lines have been generated and stored on tape for both C 12 0 and C 13 0 species. A 
deck of punched cards with the same spectral data as that recorded on tape in addi- 
tion to line half-width taken from data of James and Plyler (3) was also produced. 

The generated line positions were compared with the observational data of Connes^. 

It was found that the value of u) g of C 13 0 presented by Benedict et al. was misprinted, 
which led to a discrepancy of 40 cm 1 . The discrepancy has been corrected in the 
program. Of 26 observed lines being compared, the worst discrepancy is . 0038 cm" 1 , 
which is within the measurement accuracy. 

The generated data were further verified by comparing with a laboratory 
spectrum through the program POLAYER. The band strength of 2. 06 (cm^atm^STP) 
was taken from experimental data of Kerb et al. (5) The spectrum chosen for com- 
parison was taken from the experimental data of Plyler and Thibault (6) , whose 
experimental conditions are pathlength =9.0 atm cm, temperature = 300° K and 
pressure - l atm. Figure 2-1 shows the monochromatic transmissivity a n ri 
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Wavenumber 


Figure 2- 


. Degraded transmissivity (slit width of 0. 11 cm~l 
of the R0 and R1 lines of the first overtone of 
CO for u= 9.0 cm-atm, T= 300° K and p= 1 atm 
(calculated spectrum solid, laboratory spectrum(®) 





Figure 2-2 shows a comparison of the generated spectrum (solid curve) assuming 
a slit width of 0. 11 cm 1 and the laboratory spectrum (dashed curve). The agreement 
of these two curves is so close that they are nearly indistinguishable. 

2. 3. 2 Carbon Dioxide 

C0 2 ”1.64 The linear molecule version of the SYMTOP program can be applied to 

the molecule CC> 2 with only minor changes. Because of its symmetric structure, 

alternate lines are missing from the spectrum. The 1.64 CO band system, comprised 

of four fundamental transitions, has been calculated. Data for 1634 lines have been 

generated and stored on tape for the fundamental and first and second "hot" transitions 
12 13 

°t C0 2 (1.1% of C0 2 in the normal atmosphere) whose intensities are 

greater than one thousandth of the most intense line in this region. The parameters 
chosen to calculate the spectral data were taken from the literature; slight adjust- 
ments were made by a curve-fitting technique to improve the agreement with 
experimental line position measurements. We have compared the calculated locations 
with the wavenumbers measured by Courtoy Of 732 measured lines, 465 are 
found to be within 0. 01 cm" 1 , 722 within 0.05 cm" 1 and all within 0. 1 cm" 1 . The 
most intense lines are found to be within 0.01 cm" 1 . The line position fit is thus 
considered quite satisfactory. 

/g\ 

A paper by Boese et al; provides the quantitative measurements of line 

intensities for the fundamental transitions. The total band intensity in the program 

was set equal to the observed value 2369 x 10~ 5 (cm“ 2 atm" 1 STP) and the proper 

ratios 0.0536, 0.4432, 0.4517 and 0.0515 were assigned to the (30^1), (22^1), 

(14 1) and 06 1) bands respectively. Then a comparison was made with the 

intensities of individual lines. Of 149 measured lines being compared, 99 are found 

to be within the stated experimental error. A typical ejqperimental error is about 

5% of the observed value. The worst discrepancy from the mean measured value 
—5 —2 —1 

is 4. 2 x io (cm atm ) 2 qq ° k ( 15 % °f the line strength) compared with stated 
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experimental error of ± 2. 0 x 10~ 5 (cm^atm" 1 )^ R (8% of the line strength). 

Boese et al. also present calculated line strengths to compare with their measured 
values, and for this line they show a similar discrepancy of 13% of the mean 
measured value . 

The spectral data stored on the tape are further verified by comparing a 
spectrum generated by POLAYER with the laboratory spectrum of Boese et al. The 
laboratory conditions are pathlength = 800 atm cm, pressure = 380 mm Hg and 
temperature = 296° K. A triangular shaped slit function is assumed in the program 
with slit width chosen to be 0.08 cm . Figure 2-3 shows the undegraded spectrum 
from 6330 to 6336 cm \ In Figure 2-4 the solid curve shows the calculated degraded 
spectrum while the dashed curve represents the measured spectrum. The hot lines 
are clearly shown in the spectrum and the magnitude of absorption is well reproduced. 
We believe that the slight discrepancies between the two curves are ascribable to 
the assumption of a triangular slit function in computation. 

C °2 ~ 2<1M ‘ The Parameters u>, S and F of 1173 lines of the 2 v j + v g band systems 
have been generated through use of the linear molecule version of program SYMTOP. 
A half-width of a = o. 064 cm 1 was assumed. The generation of data for this band 
system was necessary because of the negative results of the modeling studies for 
the 1. 6M* band. The molecular constants for the 2. in band were taken from the work 
of Courtoy . A curve-fitting technique was applied to improve the computed line 

positions. The resultant accuracy of the generated spectral lines is of the order of 
0.02 cm \ 

C Q 2 ~ 5Mi The Parameters cu, s, and F have been generated by SYMTOP for 585 
lines of this perpendicular band system which overlaps both the carbon 

monoxide 4. 6n band and the nitric oxide 5. 3n band. This perpendicular band 
(V 1 +V 2 ) iS different ™ structure from the 1. 6n (parallel) band, most notably with the 
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TRANSMISSIVITY 



Figure 2-3. Monochromatic transmissivity of 
the band system at 1 . 6n of C0 2 
between 6330 and 6336 cm" 1 for 
u=800 cm-atm, p=0.5 atm and 
T= 296°K. 


Figure 2-4. Degraded transmissivity (slit width 
of 0.08 cm -1 ) of spectrum from 
Figure 2-3 (calculated spectrum solid, 
laboratory spectrum( 8 ) dashed). 




appearance of a strong Q branch. The transitions included in the program (00% - 
30 0, 00 0 - 11 0, 01 0 - 12 0 and 01 0 - 20%) are those which have been detected 
and identified in the solar spectrum. The molecular constants were taken from the 

data of Benedict published in "The Solar Spectrum from 2.8 to 23. 7 Microns " 

(9) ’ 

by Migeotte et al. We have reproduced the calculated values (typically to about 

0.03 cm ) listed in the solar spectrum except for the one transition (01% - 12%) 

with the odd lower state quantum number J, apparently due to misprints of the 

molecular constants. For this transition, a different set of molecular constants was 

taken from the work of Courtoy (?) ; the resultant accuracy of the generated spectral 

lines is comparable to the accuracy of the lines generated for the other transitions 

(typically about 0. 03 cm ). The values of the band strengths were taken from 
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C0 2 ~ 2 - 7M, » 10tx Lhie Parameters of two more parallel bands of carbon dioxide 
have been generated using the same SYMTOP program as that used to generate the 
1. 6n C0 2 band system. The 2. 7 m- and the 104 band systems are required because 
they overlap the ammonia bands (3. 04 and 10. 54). For these band systems, the 
molecular constants were taken from Migeotte et al. (9) Only the portion of the 2. 74 
band system detected in the solar spectrum (i.e. , wavenumber less than 3561 cm' 1 ) 
was generated and 150 lines were stored on tape. In the 104 region 427 lines were 
stored. For these two cases we are able to reproduce the calculated values listed 
in the solar spectrum which has a typical accuracy of about 0. 01 cm' 1 . The 
strengths of these three bands were taken from Goody^ 19 \ 

2.3.3 Sulfur Dioxide 

S °2 8,6M< The P arame ters ©, S and F of 6311 lines of the v band of SO have 

been generated based on the theory of slightly asymmetric top molecule. Values of 
the molecular constants were taken from Morino et al. (11) 
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A half-width or = 0.08 cm 1 atm 1 was assumed. Considerable effort was made to 

improve the representation of the spectra of slightly asymmetric molecules, such as 

SO Perturbations of the energy levels for the lowest values of K have been 
& 

included, with a resultant improvement in the representation of the spectrum. 

A portion of the 8. 6M- band including the band center and the R-branch maximum 
is shown in Figure 2-5, and the same spectrum viewed with a triangular slit of 
width 1 cm * is compared with the results of Burch^^ in Figure 2-6. 

(13) 

Tables of SO energy levels have been obtained from G. Steenbeckeliers ' ' 

of Belgium. Comparison of these tables shows that some of the perturbed energy 
levels are very accurately calculated by our program, while certain others are not. 
This discrepancy should be investigated, and should lead to a further improvement 
in the calculated spectrum. 

We have been informed verbally that the SO line parameters in preparation 
by R. F. Calfee may be soon available. Consequently, we have suspended further 
work on the refinement of our SOg line parameters, pending clarification of the 
availability of Calfee f s data. 

2.3.4 Nitrogen Dioxide 

NO —7. 6M- The SO version of SYMTOP has been used to generate a set of 

u & 

parameters for 2583 lines of the v band of NO . The molecular constants were 

(14) 2 

taken from Abe et al. Because of the large numbers of lines (and consequent 

heavy overlapping) we have neglected the level- splitting effects of the unpaired 

electron, which were considered in the case of NO. We have found no good spectra 

for comparison. We are, however, planning on measuring our own laboratory 

spectra of NO,^. Until better laboratory data do become available, we will not 

attempt further refinements in the calculation of the NO line parameters. 

2 
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Figure 2-5. Monochromatic transmissivity of the 
SOo-vi band between 1140 and 1180 
cm 1 (showing band center and R- 
branch maximum). The conditions 
are the same as in Figure 2-6. 


Figure 2-6. Degraded transmissivity (slit width is 

1 cm -1 ) from Fig. 2-5 for u= .577 cm atm 
P e = 1. 59 atm and T = 298° K (calculated STP 
spectrum solid, laboratory spectrum dashed). 




The NO molecule is very similar in structure to the SO molecule, and the SO 
2 2 2 

version of SYMTOP was used directly for NO . A slight amount of line splitting 
resulting from the odd electron in the NO molecule is of very minor importance 
because of the heavy line overlapping which occurs in the spectra of these molecules. 
A literature search indicated that insufficient data for NO are available for our 
needs* 18 *. 


2.3.5 Nitric Oxide 

NO - 5. 3P Parameters for 484 spectral lines in the NO fundamental band system 

have been generated by a modified version of the SYMTOP program for linear 

molecules. Because of the impaired electron, the ground electronic state of NO is 
2 2 2 

a n state. The two levels n, and n 3 describe the state of the molecule, and 

t 2 

divide the band into two subbands. Each line of the band is further split in two by 

A-doubling. Accordingly, the SYMTOP program had to be modified for this molecule. 

(19) 

The molecular constants were taken from Herzberg to generate the spectral data. 

The line locations were compared with the experimental values of James^^. Of 98 

measured lines, 35 lines are found to be within 0.05 cm 1 , 58 lines within .1 

cm , and the worst discrepancy is . 192 cm \ Although the calculated values of 

James are in better agreement with the observed values, it is not clearly mentioned 

in the literature how the set of values was computed. However, the accuracy of the 

generated data is considered satisfactory to serve our purpose. The computed 

line strengths were also compared with the measured values of Abels and Shaw^ 21 \ 

-1 -2 

Of 73 measured lines, 47 lines are found to be within .05 (atm cm ), 71 lines 

-1 -2 

within .1 (atm cm ) OAA0 . The comparison shows excellent agreement. 

oUl) K. 
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The generated spectral data were used as an input to the program POLAYER to 
simulate a standard spectrum. Figure 2-7 shows part of the unsmeared transmissi- 
vity. The real spectrum to be compared was taken from the paper by Abels and Shaw 
with laboratory conditions: temperature =300°K, pressure = 200 mm Hg and path- 
length = . 263 atm cm. By comparing the line locations obtained by James and the 
spectrum recorded by Abels and Shaw, it was found that the calibration of the Abels 
and Shaw spectrum was off by approximately 0.1 cm" 1 . In Figure 2-8 the solid 
curve represents the generated spectrum from POLAYER and the dashed curve shows 
the laboratory spectrum with corrected calibration. Other portions of the spectrum 
were also compared and the evidence shows that they are in good agreement. 

2.3.6 Ammonia 

NH 3 -I0.5ix An effective representation of the NH^ 10. SM> band system has been 
produced. SYMTOP has generated a total of 1050 lines in this spectral region. We 
have compared the calculated locations with the wavenumbers measured by Garing et 
al. Of 463 measured lines, 188 are found to be within 0. 1 cm" 1 , 351 within 0. 5 
cm , and 415 within 1.0 cm . The smallest discrepancies are quite strongly 
correlated with the most intense (hence, the most important) lines. It may be 

possible to refine these results even further; however, at the present, this is con- 
sidered quite satisfactory. 

Further verification of the spectral output from SYMTOP is provided by 
comparing a spectrum generated by POLAYER with a laboratory spectrum. 

Quantitative measurements usually have limited spectral resolution, but a good 
check is provided on the gross features of the spectrum. 

Graphs of segments of the calculated spectrum of NH^ (portions of the Q and R 
branches) are shown in Figure 2-9. These graphs are output from POLAYER and 
represent the NH g spectrum as viewed with no instrumental distortion (i.e. , mono- 
chromatic ally). The conditions of pathlength and pressure were chosen to simulate 
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Figure 2-8. Degraded transmissivity (slit width is 0. 13 cm" 1 

from Figure 2-7 for u= 0.263 cm-atm, p =200 torr 
and T= 300° K (calculated spectrum solid, laboratory 
spectrum dashed). 
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Figure 2.9 Monochromatic transmissivity of portions of the 10. 5 M' band system 

of NH 3 . A segment of the Q-branch between 900 and 950 cm -1 is shown 
on the left hand side and a segment of the R-branch between 115 0 and 

1200 cm ^ on the right hand side. The conditions are the same as in 
Figure 2-10. 



(23) 

one of the cell measurements of France and Williams . IMortunately, the authors 
did not state the slit width that they used. Graphs of transmissivity for an assumed 
slit width of 8 cm 1 were generated by PO LAYER and are shown in Figure 2-10. 

The comparison is quite good as seen in Figure 2-10. The spectrum presented 
in France and Williams' paper is extremely small and their wavenumber scale 
cannot be read to an accuracy of better than 2 or 3 cm” 1 . The comparison might 
be improved by assuming a slightly narrower slit width or a more realistic slit 
function than the assumed triangular shape. 

NHg - 3. 0M- The same version of the SYMTOP program used to generate the 

spectral data for the 10M- (v^ ammonia band has been applied to generate data for 

the 3M- 0 ^) region. The molecular constants were taken from the experimental data 
(24) 

of Benedict et al. , and the generated line positions were compared with their 
observations. 1015 lines were recorded on the magnetic tape. Of 257 observed 
lines being compared, 195 are found to be within 0.5 cm" 1 , 225 within 1 cm" 1 . The 
comparison shows the agreement being as good as that in the 10M- band. 

A simulated spectrum was generated to match the experimental conditions of 
(23) -1 

France and Williams v assuming a slit width of 50 cm . It is unfortunate that no 

higher resolution spectrum such as that shown in the 10M- region is available for 
comparison. Figure 2-11 shows the unsmeared case while the comparison of the 
simulated spectrum (solid curve) is reasonably good as seen in Figure 2-12. 

NH - 2. 3M> The ammonia 2. 3M- region consists of a parallel band v +v 
■■ 1 ■ 1 2 

(w = 4293. 716 and 4320. 06 cm ) and a perpendicular band v +v (© = 4416.908 and 

2 3 0 

4434. 610 cm ). The band strengths of the parallel and perpendicular bands are 
•2 -1 

2. 9 and 19. 7 (cm atm STP), respectively, taken from a detailed study of these 

(25) 

bands performed by Benedict et al. The set of molecular constants presented 
by them has been used to reproduce the line parameters and 3007 lines were 
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Transmissivity 
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DO 

CO 



Wavenumber (cm 1 ) 


Figure 2-10. Degraded transmissivity (slit width of 8 cm -1 )of a portion of the 10. 5 P 

band system of NH 3 between 900 and 1260 cm -1 for u= 1.617 cm atm gT p 
P e =4185 torr and T=299°K (calculated spectrum solid, laboratory 
spectrum^ 22 ) dashed). 
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Figure 2-11. Monochromatic transmissivity of NH g between 3100 and 3700 
cm -1 for u= 6.23 cm atmg Tp , p g = 5090 torr, and T= 299° K. 
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Figure 2-12. Degraded transmissivity (slit width of 50 cm -1 ) from Figure 2-11. 

Calculated spectrum solid; laboratory spectrum ( 23 ) dashed. 




recorded on tape . Unfortunately, the generated data do not reproduce their computed 
values to the same degree of accuracy achieved in the other NH bands. The reason 

O 

for these discrepancies is not clear at this moment. Comparison of the generated 
data with their observations has been made on the v +v parallel band. Of 336 

X z 

measured lines, 311 are found to be within 0.5 cm , all within 1 cm . Although 
the generated data are not as good as the computed values presented by Benedict 
et al. , the agreement is considered acceptable for this program. 

No complete comparison has been made with the Benedict et al. observations 
of the v +v perpendicular band. Of about 100 lines being compared, the majority 
of the lines are found to be within 1 cm . The line parameters could be further 
improved and an accuracy of less than 0.5 cm 1 could be expected, although the 
computed values by Benedict et al. show somewhat better agreement. However, since 
ammonia is only a possible interfering species in this region, the generated data are 
considered to be of acceptable accuracy. 

2.3.7 Nitrous Oxide 

N 04. 5M- The wing of this band overlaps the CO fundamental band. The 
z 

parameters for 976 lines in the 4.5M- band of N O (line strength, position, and energy 

z 

level) were recorded on tape. For those lines whose wavenumbers are less than 
2240 cm , a deck of punched cards (868 cards) was also produced. The parameters 
chosen to generate the tape are taken from a paper by Tidwell et all 26 ^ The computed 
wavenumbers were compared with the solar spectrum of Migeotte et al^ Of 201 
measured lines being compared, 179 lines are found to be within 0.05 cm"" 1 , 199 
within 0.1 cm 1 , and the worst discrepancy is 0.18 cm 1 . 

A spectrum generated by POLAYER is compared with a laboratory spectrum 

(27) 2 1 

produced by Plyler et al; The band strength of 1920 (cm" atm STP) was adopted 
/oft) 

from Yale et al. The laboratory conditions are: temperature = 220° K, 
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pressure = 0.00263 atm and pathlength = 0.31 atm cm. Figure 2-13 shows the calculated 
spectrum in high resolution from 2204 to 2206 cm \ Comparison of the calculated 
smeared (solid curve) and the measured (dashed curve) spectrum is given in Figure 
2-14. Absolute values of absorptivity were not given for the original laboratory 
spectrum, but the spectrum appears generally in good agreement. It appears that the 
experimental slit width (not specified) is slightly greater than the value of 0.05 cm" 1 
assumed in the computation. The apparent inconsistencies at 2204.2 cm 1 and 
2205. 8 cm 1 are believed to be due to measurement noise. 

N^O - Other Bands The nitrous oxide vibration-rotation bands which overlap 

bands of other species SO (4.04, 8.64), NO (7.64), NH o (3.04,10.54) include the 

following transitions: v +2v , 2v , 2V +v , v +v and 2v All of the bands 

1^1 2 O 1 3 2 

mentioned above are parallel bands. The SYMTOP program used to generate the line 
parameters of 4.54 band was applied directly to generate the same line para- 
meters of these bands. Of the v and 2v g transitions, molecular constants are 
available from the solar spectrum of Migeotte et al! 9 * For the v transition, the 
computed values listed in the table were reproduced from the program to within the 
given accuracy (~0.01 cm 1 ). However, for the 2v transition, the set of molecular 
constants shown in the table was unable to reproduce the same set of line parameters, 
apparently due to misprints of the molecular constants. A curve-fitting technique 
was applied to improve the computed values. The resultant accuracy in the positions 
of the generated spectral lines is comparable to the values listed in the table (typically 
about 0. 01 1 ). 

The molecular constants for the rest of transitions mentioned above are 

available from the data of Tidwell et alP^ They also computed the ling positions 

and compared with their observations. We are able to reproduce their computed 

values with a typical accuracy of about 0.01 cm \ 

A total of 1308 lines were stored on a single tape for all transitions. The 

band strengths were taken from Goody ^ except one transition 2v +v whose 

2 3 

value is not available and is estimated for this program. 
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Figure 2-13. 


Monochromatic transmissivity of N O 
between 2204 and 2206 cm -1 for u= 2 
0.31 cm atm, p =2 torr, and T=220°K. 


Figure 2-14. Degraded transmissivity (slit width 

of 0.05 cm - *) of spectrum from Figure 
2-13 (calculated spectrum solid, 
laboratory spectrum^ 27 ) dashed). 


2.3.8 Formaldehyde 


HCHO - 3 . 5M- The formaldehyde molecule (HCHO), which is an asymmetric top, has 
been modeled by the SYMTOP program. The asymmetry parameter * is nearly -1 
(-.9623), so that HCHO does not deviate greatly from being a symmetric top. 

There is a great lack of quantitative spectra for making comparisons. 

(9) (29) 

Herzberg' shows a spectrum taken by H. H. Nielsen' , but the experimental 

conditions specified are clearly incorrect. We have a spectrum taken by P. L. 

Hanst^ 30 ^ with a resolution of apparently better than 1 cm” 1 . 

The Nielsen and Hanst spectra do not seem to be completely consistent, even 
allowing for difference in resolution. 

A sample portion of the calculated high-resolution HCHO spectrum is shown 
in Figure 2-15. The large absorption peaks are Q-branches, which are well 
reproduced in wavenumber location. The finer R-branch structure in between the 
Q-branch peaks is also approximately reproduced. No attempt was made to match 
the resolution of the Hanst spectrum, which is shown in Figure 2-15 for comparison. 


2.3.9 Methane 


ch 4 Methane is a spherical top molecule, and it is this fact which results 

in considerable complication in calculating its spectrum. The symmetric top 
approximation (Eq. 2-2) predicts that all rotational levels of the same J value will 
have the same energy level. This increased degeneracy results in the appearance 


of higher-order perturbations which separate the energy levels in a manner di ffi cult 
to treat analytically. 


First-order corrections have been incorporated into SYMTOP and the program has 
generated the correct total strength and spectral location of the unsplit P-and R - 
branch lines from 2870 to 3170 cm . The problem remaining is that of representing 
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WAVENUMBER CCM- 1 3 


Figure 2-15. Monochromatic transmissivity (calculated) of HCHQ between 2880 
and 2900 cm 1 (solid curve). Degraded transmissivity (measured 
by Hanst^ 30 )) for u= 0.08 cm atm, p= 0.0079 atm and T = 300°K 
(dashed curve). 





the fine structure of these "lines" which are split into several components because 
of the molecular symmetry. Approximate expressions for the splitting will be 
developed. 


2.3.10 Water Vapor 

H^O-l. 6 , 2. 7 & 6 . 3M- Extensive tabulations of water vapor lines have been 
(31 32) 

published by NBS ’ . At present we have data on punched cards for nearly 

1000 lines in the 2 . 7M- and 6 . 3^ bands. 
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3 


MODEL STUDIES FOR GASEOUS POLLUTANTS 

(Task 2) 

3.1 THEORY OF RADIATIVE TRANSFER 

The radiation emergent from the atmosphere, designated here as E(uu), results from 
the operation of several different processes. The underlying surface emits thermal 
radiation; the atmosphere absorbs and emits radiation at a rate dependent on its 
local composition and temperature. The surface also reflects downwelling radiation, 
consisting mainly of solar radiation. Scattering of radiation by the atmosphere also 
takes place. 

The actual calculation of E(cu) is an extremely formidable task. The major 
complexities arise from the presence of the scattering processes, which alter the 
direction of propagation of individual photons in a manner which can only be predicted 
statistically. It is thus desirable to separate E(cu) into its components as far as 
possible in order to make the problem more manageable. The quantity E(cu) may be 
written as a sum of component terms in the following form: 

EW = V*) + E r W + E 0 <«) + E^c) (3-D 

in which E q (io) represents the thermal radiation emitted by the underlying surface 
and the atmosphere, E R (to) represents the incident solar radiation reflected by the 
surface, E^(oj) represents the radiation scattered by single and/or multiple- 
scattering processes in the atmosphere without having been reflected from the 
surface, and E r0 (oj), another atmospheric scattering term similar to E (u>), but for 
radiation which has undergone a reflection from the surface. 

The functional dependence of these components of E(u>) is indicated below 
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E G (<U)= E G < e ’ T ’ U) ’ T S ’ T A ) 


A 


E R (“ ) ) = E R ( p , 0 , T , u> ) 

E 0 (“>) = E 0 (9,na,r,u)) 

E R0 (a ’ )=E R0 (P ' 9 ' na ' T ' a,) 

where T , T are the surface and atmospheric temperatures, e, p are the surface 
emissivity and reflectivity, 9 is the sun zenith angle, (no) is the product of number 
of scattering centers (n) and scattering cross section (ct), r is the transmissivity of 
atmosphere, referring generically to all sources of attenuation: molecular and 
particulate absorption and molecular and particulate scattering. 

The equation of radiative transfer, without scattering and without solar radia- 
tion, for thermal radiation emerging from a plane -parallel atmosphere with an 
underlying surface of emissivity e(w) can be written 

E(<d) = E (<u) 

= e(u))N° (w,T s )t(u),0)+ C N°| id, T(z)j dz (3-3) 

*'0 

0 

in which N (u), T) is the Planck blackbody function, T is the surface temperature, 

s 

h is altitude at top of atmosphere, t(uu 9 z) is the monochromatic transmissivity of 
the atmosphere at wavenumber between altitude z and top of atmosphere, T(z) is 
the temperature of the atmosphere at altitude z. 

The general expression for the transmissivity is given by 

t (id* z) = II T (cd,Z) = exp 
i i 

where (<jo» z ') is the extinction coefficient for species i which is given by the sum of 


d 

-1 


E k (u) > z / )n. (z ') dz * 
i 1 1 


(3-4) 


(3-2) 
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the scattering and absorption coefficients, and n is the number density of contributing 
molecules of species i in the optical path. 

The contribution from sunlight reflected from the surface becomes significant 

at shorter wavelengths. Tins contribution is given by the component E (a>) : 

R 

e r W - cos 6 H g ( 0 )> ( T (W) j Y ( 3 _ 5 ) 

where [l-e(io)] is the ground reflectance of a diffuse surface, 0 i S the sun zenith 
angle, H g (u>) is the sun irradiance at the top of the atmosphere, t(ib) = <r((ju, 0) is the 
transmission vertically through the atmosphere, and y = 1 + f(9) with f(9) = sec 9 
for 0 £ 0^ 60° and Ch 9for 0 > 60°(Ch 9 is the Chapman function). 

In order to gain some insight into the radiance levels which can be expected 
at the top of the atmosphere, we estimate the contributions of the different terms 
for some arbitrary but realistic conditions. For the first term in Eq. (3-1), E , we 
assume a surface temperature of 300° K, with an emissivily of 1 for oceans and 0.6 
for land (see Figure 3-1). For simplicity in the figure, we have not included the 
effect of atmospheric absorption and emission. In general, this inclusion would 
lower the radiance level in the region of atmospheric absorption bands. 

For the second term, E R , we assume a sun zenith angle of 45° , a reflectivity 
of 0.02 for oceans and select 0.4 as an upper value of land surfaces, and we also 
assume that the transmission is given by the extinction due to Rayleigh and Mie 
scattering, and ozone absorption. The extinction is taken from the Elterman^ model. 
For the third term, E^, we refer to the calculations of Plass and Kattawar* 2) , who 
have calculated the reflected radiance for a zero albedo as a function of 
wavenumber, sun angle and different amounts of Rayleigh and Mie scatterers. How- 
ever they did not include in their multiple scattering calculations the effect of 
molecular absorption (except for Cy. The fourth term, E^, has not been calculated 
explicitly for an atmosphere including particulates. However, it is contained in the 
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Wavelength (H) 


Figure 3-1. Upwelling radiation at the top of atmosphere. Contributions from 
reflected sun radiation E R (uu), atmospheric scattering E 0 (u>), and 
thermal radiation Eq(uj). Ground temperature = 300° K, sun zenith 
angle = 45° , particulate matter according to Elterman 1968. 
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calculations of Plass and Kattawar for surfaces with a given constant albedo, extend- 
ing to infinity. The contribution of this term, especially from surfaces with an albedo 
different from the surface which is in the field-of-view, will be the subject of 
further investigation. 

It is seen that the amount of sun radiation scattered from the atmosphere is 

comparable to the amount reflected from the ground. The relative magnitude depends 

upon the ground reflectivity and the number of scatterers. Thus, in a "clean" 

atmosphere with the standard Elterman distribution of aerosols and a high ground 

reflectivity of 0.4, the contribution of the atmospheric scattering(E )is about 4 per 

0 

cent of the amount reflected from the ground (E ) for \> 24. (For X< 2H the con- 
tribution of E 0 increases rapidly and exceeds E r for \ < 0. 4P-). On the other hand, 
when the atmospheric content of aerosols is larger and the surface reflectivity is low, 
such as for water, the contribution of E^ becomes comparable to the amount of E r . 
Thus, it becomes essential to know the term E^ before the reflected sun radiation 
can be utilized for a measurement of a particular pollutant. 

In the thermal region, the situation is analogous. The emission from the 
ground is dependent upon the ground brightness temperature. The species in the 
atmosphere absorb and emit depending upon the local temperatures and optical thick- 
ness. Thus, it becomes essential to know these ground and atmospheric parameters 
before the thermal region can be utilized for a measurement of pollutants. 

3.2 DESCRIPTION OF GAS FILTER CORRELATION INSTRUMENT 

The basic principle of operation of a gas filter correlation (GFC) instrument is to 

pass an undispersed beam alternately through one of two paths terminating on the 
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same detector: the first, through a cell containing a sample gas, and the second 
through a vacuum cell and an adjustable aperture which is used in zeroing the 
instrument. The output of the instrument is the difference of the signals from the 
two channels. 

A detailed analysis of a characteristic GFC instrument was presented in 
Reference 3. The precise form of the equations depends on the particular geometry 
of the instrument, since each mirror reflectivity, window transmissivity, etc., 
must be considered. However, all these items, along with instrument responsivity 
may be lumped together into certain generalized instrument response functions. 

The signal change AV is expressed by the equation 


AV = A Or 


/e ( ») [c(») 


T g (W) - C (CD) 



in which A Oand R Q represent the instrument throughput and responsivity. The 
quantities C(cd) and C 7 (ou) represent generalized filter functions and contain trans- 
missivities and reflectivities of various components, as well as the transmissivities 
of filters which have been selected to restrict the operation to spectral regions where 
absorption bands of the species of interest exist. 

The quantities T G ( tja ) and represent the monochromatic transmissivity 
through the gas cell and the transmissivity of the adjustable aperture. This aperture 
is set by using a calibration blackbody source and adjusting the aperture until a zero 
signal is outputted by the instrument. The value of t is given by 

A 


T 

A 


/ 


0 

N (CD, T ) C(tD) T (U)) d(D 
c G 


f 

•'Acd 


N° (cd T ) C‘ 
c 


(CD) dcD 


where N°(cd, T ) is the Planck blackbody function at the calibration blackbody 
c 

temperature T . 
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3.3 


DESCRIPTION OF THE POLAYER PROGRAM 
3.3.1 The Basic POLAYER Program 

The designation POLAYER refers to a family of programs for calculating atmospheric 
radiance monochromatic ally. The calculational procedure is basically the same 
in all versions; the differences are primarily in the input/output format. The opera- 
tion of the basic POLAYER program will be described first. 

The atmosphere is divided into a number of layers (perhaps 10-20), each 
of which is assumed to be homogenous in temperature, pressure, species concentra- 
tion, etc. The earth is described by a Lambert surface whose emissivity is inputted 
as a tabulated function of wavenumber. The reflectivity, assumed to be one min us 
the emissivity, is used to determine the reflected sunlight, but not (because of the 
computational difficulties) the reflected atmospheric radiance. 

The program provides calculations over a specified wavenumber 
range. At the initial wavenumber, the equation of radiation transfer is evaluated 
at several points j at the top of the atmosphere (i.e., in front of the instrument), 
at the exit of a gas-filled instrument cell, and at the exit of a vacuum cell con- 
taining an adjustable aperture. 

An increment of wavenumber is determined by the program; this increment 
is determined by the closeness of spectral lines. When lines are widely spaced, a 
coarser mesh is permitted; as a spectral line is approached, the mesh becomes 
fmer in order to reproduce the spectral structure about the line center. This internal 
selection procedure can be over-ridden by specifying a minimum wavenumber 
increment; this procedure is useful, for example, in simulating cell measurements 
when some minimum pressure (e. g. , of the order of one atmosphere) in the line of 
sight is known and thus rapid spectral changes will not occur within spectral intervals 
smaller than some specifiable minimum (e.g. , 0.02 cm"" 1 ). 
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As the program proceeds through the selected wavenumber interval, various 
integrals are calculated by the trapezoidal rule. These integrals include those of the 
radiances previously described as well as certain other useful integrals of cell 
transmissivities, atmospheric transmissivities, etc. Aside from these integrals, all 
other calculated quantities are erased from memory when a new wavenumber is 
selected. 

At each wavenumber and in each atmospheric layer, the absorption coefficient 

of each spectral line is calculated and stored in memory. The lines are assumed to 

be of Voigt shape (i.e. , a convoluted Doppler-Lorentz shape). A convenient 

(4) 

approximate form developed by Whiting' 'has been used. This form has been shown 
to reproduce the Voigt shape very closely under most conditions, and with an 
absolute maximum error of 5% occurring at near-zero pressures. 

A major part of the calculational effort goes into the evaluation of these 
absorption coefficients. Since these are erased after they are used, we have made it 
possible to calculate a set of results for different pollutant concentrations, either 
in the atmosphere or in the instrument cell. Thus a set of several calculations may 
be made at only a very moderate increase in computing cost. This feature is very 
useful for conducting sensitivity studies. 

The program also simulates the calibration of the instrument cell by an on- 
board blackbody of specified temperature and the setting of the vacuum cell aperture 
to the balance condition as described in Section 3.2, 

The integrals are weighted by an instrument response factor which is read in 
as a tabulated function and incorporates the detector response and the transmission 
functions of the filters being used. 

3.3.2 Test Calculation with POLAYER for the Solar Spectrum at 4. 65M- 

The basic POLAYER program was tested through a calculation of the solar spectrum 

between 2136 and 2166 cm . The line parameters of the fundamental band of 

C O , C O and the v fundamental band of H O were inputted^** The other 

« 2 
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atmospheric constituents CO and N O were not included in this calculation. The 
results of the calculations are compared with the measurements of the solar spectrum 

(7) 

by Migeotte et al. The wavelength interval chosen involves the three plates 28, 29 

and 30, for which the sun elevation was 25-32° , 19-25° and 24-30° , respectively. 

We have used a mean value of 25° for the calculation with POLAYER. In addition 

to the atmospheric absorption lines, the experimental spectrum contains the 

Fraunhofer lines of CO. The concentration of CO used in the calculation was 0. 1 ppm 

(7a) 

in the troposphere and 0.001 ppm in the stratosphere. The water vapor concentration 
at the time of the measurements was listed as 58% relative humidily. In the first 
test calculations, the standard Gutnick concentration profile was used. However, 
the comparison between the calculated and experimental results showed that the 
Gutnick water vapor concentration was too high by a factor of 5. The slit width was 
listed as 0. 2 cm 1 , which was also used in the calculations in conjunction with a 
triangular slit. The results of the calculations with and without a slit function are 
shown in Figure 3-2. The comparison of the calculated spectrum (degraded by the 
slit function) with the experimental spectrum shows good agreement, except for a 
number of weak lines, which are primarily N O absorption lines. 

Li 

3.3.3 Versions of POLAYER 

One version of POLAYER has been adapted specifically to the drawing of graphs 
(e.g. , of transmissivity vs. wavenumber). Either (or both) the true monochromatic 
transmissivity or a transmissivity degraded by a slit function of triangular shape 
can be generated. This version is most useful for verifying the spectral data written 
on tape by comparing an experimental spectrum with the synthetic spectrum generated 
for the same set of experimental conditions. 

Another version is designed to calculate the signal outputted by a GFC 
instrument without the plotting capability previous ly described (which requires 
significant core storage). 
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Different versions also exist to accept spectral line data in varying degrees 
of refinement. One version (the earliest stage of development) will accept randomly 
listed (tape or card) spectral line data, sort these by wavenumber, and store them in 
memory. On very large programs, since all the spectral data are stored in memory 
all the time, this results in excessive demands on core storage. 

A modification was introduced to reduce memory usage. A separate program 
is used which merges the spectral data from a number of tapes, sorts them, and 
writes them on a new tape. This tape is then used with a modified version of 
POLAYER, which reads in and stores data for a limited number of lines at a time, 
then reads in new data as needed, while erasing the old from memory. In this way, 
the required field length is kept to a minimum. At the present time three operative 
versions of the POLAYER program are: 

POLAYER I: Includes an automatic plotting routine used for the comparison 

of calculated and experimental data; and includes a sorting 
routine for input tapes. 

POLAYER II: Automatic plotting routine removed; but includes the same 
sorting routine as POLAYER I. 

POLAYER III: No plotting and sorting routine; accepts only one sorted input 
tape, which is the output of a separate sorting routine. 

Sorting programs have been developed which accept several tapes, each con- 
taining spectral line parameters for only one molecule (unsorted as to wavenumber) 
from program SYMTOP and create a new tape used as input to POLAYER in. On 

this tape, all lines are sorted by wavenumber, and the line strengths are entered in 
absolute units. 

Since it is no longer necessary to perform the line sorting in POLAYER, only 
a limited number of lines need to be stored in memory at one time, greatly decreas- 
ing the field length requirements, and, consequently, the cost of POLAYER runs. 

The present status of the logic schematic diagram is shown in Figure 3-3. 
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Figure 3-3. Flow diagram of SYMTOP/POLAYER logic 


3.3.4 Future Development of PO LAYER 

As stated in 3. 3. 3, the line-by-line calculations are very time-consuming when 
several thousand lines are involved. The central processor times become 
prohibitively long when many model calculations are required. Thus, it becomes 
imperative to investigate different approaches. Before we resort to the band-model 
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approach, different integration schemes are being investigated. In the present one, a 
linear interpolation between adjacent points on the spectral line contours is made, 
which requires very small steps when the radius of curvature is small, as it is near 
the line centers. An alternate method is to use a quadratic fit between three selected 
adjacent wavenumbers to achieve a better fit with coarser spacings. Different selection 
criteria for the points to determine the parabolae should be used for line centers, 
line wings in the presence and absence of neighboring lines, etc. As an example, 
computations for two neighboring lines, one from CO and the other one from H O, 
have been performed with the integration schemes. The results are shown in Figure 
3-4. The results in terms of line contours are the same but the number of points 
considered is reduced from 49 in the linear interpolation scheme to 25 values in the 
quadratic fit. 

In the present program, the number of contributing lines is fixed, e. g. , at 
20 or 40, depending on the spacings and strengths. It is realized that the selection 
of the number of contributing lines should be more flexible and should be based on 
the actual wing contribution. This procedure may not decrease the integration time 
in certain cases, but the calculations would be more precise. 

3.4 RESULTS 

In the following model calculations, using Eq. (3-1), the contribution of the scattering 
terms E^ and E^ r have not been considered, because the required Monte Carlo 
calculations, which include the absorption by the pollutant molecules, have not been 
performed yet. 

3.4.1 Carbon Monoxide 

\_ ■ - — - 
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Figure 3-4. Comparison of 


■■■ 



3. 4. 1.1 Comparison with Band Model Calculations. We have compared the line-by- 
line calculations with band model calculations. The signal change in radiance inte- 
grated over the band plass between clean and polluted atmospheres has been calcu- 
lated in the line-by-line case to be about 5. 75 per cent. The band model calculations 
are made in spectral intervals of about 0.1 |i. If the band model results are averaged 

over the band pass, the mean signal change is about 20 per cent less than the mean 
* 

signal change of the line-by-line calculation. In this case, then, the band model 
calculations are conservative. More test cases, also with other molecules, have to 
be run to establish the utility of band model calculations in the majority of sensitivity 
calculations. Of course, to calculate the possible influence of other interfering gases 
and the degree of "randomness of correlation", only line-by-line calculations can 
be used. 


3. 4. 1.2 Influence of Stratospheric CO. The total atmosphere is divided into 20 
layers, 8 for the troposphere and 12 for the stratosphere. The distribution of CO 
was assumed to be 0.025 to 0.4 ppm in the troposphere and 0.001 ppm in the strato- 
sphere.^ ^y restricting the atmosphere to 8 layers (troposphere only), the signal 
difference is affected by less than 0. 3 per cent for 0. 1 ppm and less than 0.05 per cent 
for 0.4 ppm CO in the troposphere. Thus, it was decided to do the present calcula- 
tions with only 8 layers since the computing time is cut down by more than a factor 

of 2. Of course, for a detailed investigation to be conducted later of the influence 
of the stratosphere (especially of stratospheric H O) and of different CO mixing ratios 
in the stratosphere, the program will be rim with 20 layers. 

3.4. 1.3 Influence of Water Vapor. The atmospheric profile used in these r un s 
together with the water vapor distribution is shown in Figure 3-5. The temperature 
and pressure profile is the Standard Atmosphere and the water distribution (fractional 
concentration) is taken from the Gutnick^ model called "standard wet" from here 


♦"Signal Difference" is denoted interchangeably as DV or AV. 
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Fi^ire 3-5. Temperature, pressure and water vapor profiles used in the 
model calculations. 


47 




!»■■■■■ 


iHsssasss 




riifliimni 

hifntmmmamm 


It am 


I Nil 


Figure 3-6. Signal difference vs. CO concentration in troposphere for different 
amounts of water vapor 
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on. The calibration temperature and ground temperature are 300° K and the ground 
emissivity is unity. A total of 90 CO lines and 108 H O lines were included. A 
comparison of the results between three different water concentrations and the dry 
atmosphere is shown in Figure 3-6. As expected, the AV is reduced by the presence 
of water vapor in the atmosphere. This reduction between dry and standard wet 
atmosphere amounts to about 23% for a CO concentration of . 1 ppm. The reduction 
is caused firstly by the transmission loss due to H O and secondly by the "negative 
correlation, " which can be studied when (CO) =0. In that case, Av is negative as 
seen on Figure 3-6. It means that the water lines are not completely randomly 

distributed with respect to the CO lines, but correlate with the spacings between the 
lines of CO. The relative effect of both mechanisms on the reduction of AV is 
shown in Figure 3-7. If the AV originating from CO and the standard wet atmosphere 
is divided by q (as calculated by the program), the dashed line is obtained. One 
can see that the transmission loss is the predominant cause at the higher CO con- 
centration and the "negative correlation" is the primary cause at the lower 
concentration. 

3.4. 1.4 Influence of Stratospheric H^O on CO Measurements. Test runs were made 

to study the influence of stratospheric water vapor especially in view of neglecting 

the 12 stratospheric layers. The water distribution was taken from the Gutnick^ 

model up to 31 km and a constant mixing ratio above 31 km as shown in Figure 3-8. 

The results show that the influence of the stratospheric water vapor is negligible. In 

-7 2 

Table 3-1, the DV in 10 w/cm sr are given for the 8 and 20 layered wet atmospheres 
and for five different CO concentrations in the troposphere and constant 0.01 ppm 
evenly distributed in the stratosphere. 
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Table 3-1. Difference of Signal Changes Due to CO with and without Water 
Vapor in the Stratosphere 


CO 

ppm 

Troposphere 

Only 

With Strato- 
sphere 

Difference (%) 

.025 

2. 1400 

2.1825 

2 

.05 

3.9028 

3.9281 

.6 

.1 

6.1361 

6.1482 

.2 

.2 

8.6514 

8.6563 

.05 

.4 

11.2950 

11.2960 

.01 


It should be noted that these results are conservative since this model assumes 

( 9 ) 

relatively high mixing ratios for H O. The mixing ratio given by Mastenbrook 
is a factor of 10 smaller at 15 km altitude. 

3.4. 1.5 Model Calculations for CO with H O and N O. The standard atmospheric 

profile with a ground temperature of 300° K and a calibration temperature of 300° K was 

used to calculate the signal difference for different amounts of CO in the presence 

of H^O and N 2 0 having a constant mixing ratio of 0. 25 ppm. The value for zero 

concentration corresponds to a (negative) correlation of either H O or N O or both 

2 2 

with the CO spectrum in the instrument cell. These values are shown for 10 ~^ DV 
in w/cm 2 sr in Table 3-2. 
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Figure 3-7. Signal difference vs. GO concentration with and without normal 
water vapor, indicating "negative correlation. " 
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Figure 3-8. Water vapor distribution according to Gutnick 



Table 3-2. Signal Change for CO in the Presence of the Interfering 

Molecules HO and N O 
2 2 



CO Only 

H 2 0 Only 

N 2 0 Only 

CO & H O & N O 
2 2 

For (CO) = 0 

0 

- .42 

-.06 

-.45 

For (CO) = .1 

8.06 

6.19 

7.15 

5.88 

% 

0 

-6.8 

0.8 

-7.6 


The sensitivity curves for detecting CO in the presence of H O and N O are shown 

2 2 

in Figure 3-9. 

3. 4. 1.6 Model Calculations for CO with H O, N O, CO and O . For the complete 

2 2 2 3 

evaluation of the influence of interfering gases on the signal differences for CO, CO 

2 

and O must be included. Although the line parameters for the v +v band of CO 
a 12 2 

have been generated, the model calculations have not been performed yet. The 
line parameters for the 2v^ band of O^ must be generated yet. Although this band is 
quite weak, it will influence the signal difference for CO to some extent because 
of the stratospheric ozone layer. It should be noted that this band is of little 
significance when measurements are made from aircraft platforms. 
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3.4.1. 7 Comparison Between DVand Radiance Measurements. For a concentration 
range from 0 to 0.4 ppm of CO in the dry troposphere, theDV normalized at the value 
for (CO) = 0.1 ppm is shown as the upper curve in Figure 3-10. The lower curve is 
the result for the standard wet atmosphere with N O, also normalized. As mentioned 
previously, the reduction in signal difference is partly due to transmission loss and 
partly due to "negative correlation. " One can see that the signal is much more 
dependent upon the CO than the interfering gases H O and N O. 

The solid line labeled "V" is the radiance integrated over the wave n umb er 
interval (2070 to 2220 cm 1 ) as calculated for the dry atmosphere, normalized at 0.1 ppm. 
The dashed line labeled " V " is the result for the standard wet atmosphere with N 0 0. One 
can see that in this case the radiance is a slowly varying function of CO concentration, 
but that the radiance is strongly influenced by the interfering gases. 

3.4. 1.8 The Radiance as a Function of Sim Angle, Ground Emissivity and Temperature. 

One of the outputs of the POLAYER program is the integral of the radiance at the 

entrance of the instrument, i.e. , / E(o>) dw, where E(u>) is given by Eq. (3-3) and 

Acu 

Aiu by the limit values of the spectral bandpass of interest. A square filter function is 
assumed, When e =1.0 for all u), there is no reflected sun energy. For any 
e / 1.0, the reflected sun energy is dependent upon the sun zenith angle 9, i.e. , 9=0 
when the sun is overhead, and 9= 90 ° represents night time observation. Since E(cu) 
is a linear function of e, only one calculation per sun zenith angle was performed. 

If the calculation is performed for e= 0, the dependency of the ground temperature 
is removed. Thus, a set of straight lines can be generated from calculations for 
e = l with the ground temperature as a parameter and for e = 0 with the sun zenith 
angle as a parameter. 

Several runs with e = l, 0.9, 0.8, 0.6 and 0.4 were performed and the 

linearity verified. Representative plots for T_ = 270, 300 and 330° K with a 

GR 

constant CO concentration of 0.1 ppm are shown in Figures 3-11 through 3-13. 


55 



Normalized DV and V 
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Figure 3-10. Normalized signal difference and radiance vs. CO concentration 
with and without water vapor. 
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Figure 3-11. Radiance at the top of atmosphere vs. ground emissivity for 
different sun zenith angles at a ground temperature of 270K 
and CO concentration = 0.1 ppm 
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Figure 3-12. Radiance at the top of atmosphere vs. ground emissivity for 

different sun zenith angles at a ground temperature of 300K 
and CO concentration =0.1 ppm 
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Figure 3-13. Radiance at the top of atmosphere vs. ground emissivity for 
different sun zenith angles at a ground temperature of 330K 
and CO concentration =0.1 ppm 
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g 

The ordinates are labeled RR1I x 10 , which is the computer symbol for 
E(co) da: w/cm sr. The atmosphere profile was kept the same. As the 

temperature difference between ground and atmospheric temperature increases, 
the radiance becomes larger. Note that different scales in the figures are used. 

In Figure 3-14 a cross plot of the radiance versus the ground temperature for e = l 
and e = 0. 8 and for four different sun zenith angles is presented. This plot illustrates 
the rapid increase for increasing temperature differences between ground and 
atmosphere. 

3. 4. 1. 9 AV as a Function of Sun Angle, Ground Emissivity and Temperature. Similar 

to the Figures 3-11 through 3-13 are the graphs for AV as a function of e, T and 0. 

G 

Straight lines are obtained which are cross plotted for AV as a function of CO con- 
centration with e and 6 as parameters. In Figures 3-15 through 3-17 are shown the 
results for 270, 300 and 330° K and a ground emissivity of 1 and 0. 8, representing 
typical values for water and land surfaces, respectively. One can see that for a low 
ground temperature the AV is increasing with decreasing sun zenith angle, but not so 
for high ground temperatures. This is illustrated in Figure 3-18, where AV is plotted 
as a function of T qr with c and 9 as parameters. The line for e = 1. 0 crosses over 
at the higher temperatures. The night time observations are represented by the 
9= 90° curves, which are always lower. 

For the cases where the ground temperatures are different from the 
temperature of the calibration source, deviations in AV are observed. Since these 
deviations are related to the instrument performance, they are discussed in Section 5. 

3.4. 1. 10 Influence of Variable Ground Emissivity on AV. The effect of a variable 
ground emissivity versus an averaged (w independent) emissivity on AV was 
investigated. The emissivities were taken from the measurements by Hovis (10) 
for two typical surfaces, silica sand and Pawnee Grassland soil. The spectral 
emissivity for the silica sand varies from 0. 98 at 5M- to 0. 78 at 4. and the average 
emissivity is found to be 0. 8426. In Table 3-3 we have tabulated the results (AV) 
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Figure 3-14. Radiance vs. ground temperature for e= 1.0 and e= 0.8 for 
various sun zenith angles 
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Figure 3-15. Signal difference vs. CO concentration for different sun zenith 
angles and e= 0.8 and s = 1.0 at ground temperature = 270K 
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Figure 3-16. Signal difference vs. CO concentration for different sun zenith 
angles and e = 0. 8 and e = 1.0 at ground temperature = 300K 
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Figure 3 17. Signal difference vs. CO concentration for different sun zenith 
angles and e = 0. 8 and e = 1. 0 at ground temperature = 330K 
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Figure 3-18. Signal difference vs. ground temperature at e = 1.0 and e = 0.8 
for various sun angles and for CO concentration =0.1 ppm 
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for the silica sand surface, using a dry atmosphere, and a ground temperature and 
calibration temperature of 300° K: 


Table 3-3. Signal Change for CO for Varying and Constant Ground 
Emissivity 


Sun Angle 8 

= 0° 



CO 

e(w) 

e 

Difference (%) 

.05 

-7 

5.58 X 10 

-7 

5.53 X 10 

.9 

.1 

-7 

8.16 X 10 

-7 

8.11 X 10 

.6 

.2 

-6 

1.10 X 10 

l.io x io " 6 

0 

.3 

— fi 

1.27 X 10" 

1.27 X 10~ 6 

0 

.4 

1.39 X 10“ 6 

1.39 X 10“ 6 

0 

Sun Angle 8 

= 60 ° 



.05 

5.08 X 10' 7 

-7 

5.03 X 10 

1 

.1 

7.41 X 10~ 7 

-7 

7.39 X 10 

.3 

.2 

-7 

9.99 X 10 

-6 

1.00 X 10 

.1 

.3 

1.16 X lo" 

1.16 X 10~ 6 

0 

.4 

r* 

1.27 X 10" 

1.27 X 10 

0 


Similar results are obtained with the Pawnee Grassland soil, i.e. , the 
difference between spectral emissivities and averaged emissivity input is insignificant. 
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3.4.2 C arbon Dioxide 


3.4. 2.1 Model Calculations for CO at 1.6M-. Calculations for the condition of a 

standard atmosphere with 320 ppm of CO ocean surface with a reflectivity of 0. 02 

and a sun zenith angle of 30° show that the sensitivity of the GFC instrument for 

measuring a 1 % concentration change in the 1.6M- system (22^1) is too low. The band 

is weak so that the average transmission through the instrument cell for 1 atm CO 

and an optical length of 10 cm and 100 cm is 0.9995 and 0.9897, respectively* which 

is too far removed from an optimized condition. To optimize the cell conditions, an 

3 

optical thickness of 10 cm atm would be required, which is impractical for space 

application. The signal change for 320 ppm CO in the atmosphere is about 
-9 2 

6 x lo w/cm sr using a cell transmission of 0.9897. This signal change is further 
reduced by about two orders of magnitude when changes of 1% in CO concentration are 
to be measured. 

3.4.2. 2 Model Calculations for CO^ at 2. 1M>. The model calculations with the 
20^1 band system at 2M- show the potential of measuring the CO concentration to 

Ct 

within one per cent. For this band the optical thickness needs to be only about 

50 cm atm for an optimized condition, resulting in an average cell transmission of 

0. 76. However, a cell length of 50 cm would be impractically long. Thus, the 

present calculations are performed with a cell length of 20 cm and a total pressure 

of 2.5 atm. More model calculations have to be performed to investigate the influence 

13X6 

of the higher than atmospheric pressure. The use of the isotope C O , which is 

12 16 ^ 

established to have an abundance of 1. 1 per cent of that of C O , was considered. 

2 

A stronger band could be used and, thus, the cell length could be reduced. However, 

the two strong bands are located at 2. 7^ and 4. 3H, which are heavily overlapped by 

water vapor and influenced by thermal emission, respectively. More work needs to 

be done in this area. The dependence of the signal change on the concentration is 

shown in Figure 3-19 normalized to the signal resulting from a concentration of 300 ppm. 

-7 2 

The signal level for 325 ppm of CO is about 2 x 10 w/cm sr. 

u 
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[Ili] 



Figure 3-19. Signal difference vs. CO 2 concentration for sun zenith angle 
30° , ground reflectivity = 0.02, ground temperature = 300K 


0 - 
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However, these results are obtained without introducing uncertainties in the 
ground and atmospheric variables. In addition, the atmospheric scattering term was 
not included in the calculations because the functional relationship of the absorption 
by CO is not known at present and must be established through further Monte Carlo 
calculations. 

3.4.3 Sulfur Dioxide 

3. 4. 3. 1 Model Calculations for SO at 8. 6M- . A normal U. S. wet atmosphere with a 

Li 

ground temperature of 300° K, and a SO concentration of 2 ppb was assumed to 

determine the signal difference. A square filter from 8.4034^ (1190 cm ) to 

9.1743M' (1090 cm 1 ) was used in the instrument. In this region, about 4000 lines of 

SO and 37 water lines are included. Additional interfering lines from N O and O 
& 2.3 

have not been considered yet. The calculations show that the signal difference is in 

2 

the order of 1. 1 x io w/cm sr for a cell condition of 5 cm length, 100 per cent SO 
concentration and 1 atm total pressure. However, these cell conditions are not yet 
optimized. We estimate that the product of cell length x concentration should be 
approximately .1 cm-atm for optimum cell conditions. More model calculations need 
to be performed to establish the final cell parameters and to determine the change of 
signal difference as a function of different SO concentration profiles and the utility 
for using two cell pairs. 

3.4. 3.2 Model Calculations for SO at 4^. Using a ground temperature of 300° K, a 

ground emissivity of 0. 8, a sun zenith angle of 30° and no interfering gases, a 

—8 2 

signal change of -1. 3 x io w/cm sr due to a SO concentration of 2 ppb was 
calculated. The negative sign appears to be the result of the wavelength dependency 
of the reflected sun light (-6000° K), which is very much different from the one of 
the 300°K calibration source used for balancing the two channels. More analytical 
work needs to be done to understand cle arly what that means in terms of data 


69 



interpretation. Also, the interfering gases N O, CH and perhaps N must be 

2 4 2 

included in future calculations. 


3.4.4 Nitrogen Dioxide 

3. 4.4. 1 Model Calculations for NO at 7.6^. A normal U. S. wet atmosphere with a 

ground temperature of 300° K and a NO concentration of 2 ppb was used to determine 

the signal difference. A square filter from 7. 1429 H (1400 cm" 1 ) to 7.6923M- 

(1300 cm 1 ) was used in the instrument. In this region, about 1500 NO lines are 

2 

included, in addition to water and NO lines. Additional CH lines have not been 

2 4 

included yet. The calculations show that the signal difference is in the order of 
-7 2 

4. 5 x io w/cm sr for a cell condition of 3 cm length, 1.66 atm total pressure 
and 50 % NO^ As in the case for S0 2 , the cell conditions of 1 cm length with 

10 percent of N0 2 at 1 atm total pressure will be near optimum. The difference 
between measuring N0 2 and S0 2 is in the interference of H 2 0, which is much 
stronger in the case of NO than for SO . However, we believe that the effect of 

<£ £i 

H 2 0 can be reduced significantly by the ratio technique. (See Section 5.4. ) More 
model runs are required to firmly establish this. 

3.4.4. 2 Model Calculations for N0 2 at 3. 3M-. The conditions used for the model 

calculations of SO at 4M- have been used here. The calculated signal difference is 
-8 . 2 

-2 x io w/cm sr, which is slightly more negative than the result obtained by for 

S0 2 . The apparent reason for the negative sign is the same as given for SO - 4|i. 

2 

Further model calculations are needed, which must include the interfering molecules 

CH and HO. 

4 2 

3.4.5 Nitric Oxide 

This molecule is very heavily overlapped by H 2 0 and previous band model calculations 
had indicated that it will not be possible to observe NO in the atmosphere. The 
present line-by-line calculations tentatively confirm this conclusion. Using a U. S. 
wet standard atmosphere, a ground emissivity of .8 and ground temperature of 300°K 
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the signal differences for different amounts of NO are the following: 

Table 3-4. Signal Change for NO as a Function of Concentration 

NO(ppb) _0 OJ 1 10 

4.4807 -4.4699 -4.3733 -3.4677 

1 .9976 .9760 .7739 


DV X 10 
Normalized 


The negative sign arises from the presence of water vapor due to the negative 

correlation of its lines with the spacings between the NO lines. This effect may be 

reduced considerably by the ratio technique. However, the signal change due to a 

-9 2 

10% concentration change at 1 ppb is less than 10 w/cm sr, which appears to be 
insufficient. 

3.4.6 Ammonia 

3. 4. 6.1 Model Calculations for NH at 10. 5 P. The calculations show that the measure 

O 

ment of NHg at 10. 5 M- in concentrations of less than 1 ppb is feasible. Calculations 
were made assuming a ground temperature of 300° K and ground emissivity of 1. 0 

and the standard dry atmosphere. The signal difference for 20 ppb of NH is 

-5,2 3 

2. 33 x 10 w/cm sr. More detailed calculations have to be made in order to 

establish the sensitivity as a function of different concentration profiles, and ground 

conditions and the effect of the interfering gases N O, CO , H O and O . 

Cl CI Cl o 

3.4. 6. 2 Model Calculations for NH„ at 3M-. In the calculations, a ground temperature 

O 

of 300° K and a ground emissivity of 0. 8 was assumed. The sun zenith angle was 

assumed to be 30° . The signal difference for 20 ppb of NH without interfering gases 

-7,2 3 

was calculated to be - 2.1 x io w/cm sr. More model calculations, which must 
include the interfering species CO , N O, CH and H O, are needed to establish 

Ci Ct 4 Cl 
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the influence of the atmospheric and ground variables and the X dependency of the 
reflected sun radiation on the balance condition. 

3.4.7 Formaldehyde 

The present calculations tentatively show that the measurement of 2 ppb of formaldehyde 

-7 2 

at 3.5^ appears to be feasible. The signal difference is about - 2 x 10 w/cm sr 

with optimized cell conditions but without the inclusion of interfering gases. More 

model calculations are needed to determine the influence of H O and HDO and the X 

2 

dependency of the reflected sun radiation. 
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4 


MODEL STUDIES FOR AEROSOL DETECTION 
(TASK 3) 

The objective is to determine the feasibility of measuring the atmospheric aerosol 
content by one or more separate radiometer channels to be included in the GFC 
instrument. 

4.1 EARTH-ORIENTED OBSERVATIONS 

Previous work at Convaii^ad indicated that satellite measurements of (a) the 

intensity of backscattered radiance, and (b) the polarization of backscattered radiance, 

would yield information on the aerosol content of the atmosphere. Calculations of 

(2) 

these parameters exist in the literature and further radiance intensity calculations 
were performed for us under another NASA contract, No. NAS1-10466, by Plass and 
Kattawar; these calculations were used to assist in determining the feasibility of 
measuring the atmospheric aerosol content with the GFC instrument. 

4.1.1 Radiance Observations 

The calculations by Plass and Kattawar made for Convair were designed to determine 
the dependency of the backscattered (reflected) radiance, as a function of wavelength, 
on the aerosol content of the atmosphere with an underlying ocean surface, assumed 
to be a specular reflector. The aerosol distribution and the total number of 
aerosols were varied. It was found that the upward radiance depends strongly on the 
total number of aerosols, but not on their vertical distribution. This is an important 
result. The total number of aerosols is a preferred parameter to measure since it 
is directly related to the mass loading of the atmosphere. 
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The previously published calculations of Plass showed that the upward radiance 
is much more sensitive to the aerosol content when the surface albedo is zero than 
when it is unity. Thus, the calculations performed for us by Plass and Kattawar 
were confined to an underlying ocean surface, which has a low albedo (. 02) for high 
sun elevations. In addition to providing more sensitivity for an aerosol measurement 
the ocean surface has known reflectance properties (see discussion of this assumption 
in Section 4.1.5). 

A detailed discussion of these calculations is given in our report "Air Pollution 
Measurements from Satellites" being prepared under NASA contract No. NAS1-1-466, 

Figure 4-1 shows the variation of normalized upward radiance with m- q (cosine 

of solar zenith angle) as a function of wavelength for the normal aerosol content 

(1968 Elterman)and three times normal. The calculated points show some scatter 

about the smooth curves due to the statistical nature of the Monte Carlo method. The 

curves are not extended into the M- Q = .95 to 1.0 region where the specular reflection , 

from the ocean surface would be seen by the observer. 

* 

The three wavelengths considered show comparable sensitivity to aerosol 
changes. However, the normalized radiance is less at the longer wavelengths, and 
since the incoming solar flux decreases at longer wavelengths, the absolute radiance 
level decreases rapidly with increasing wavelength (the absolute radiance at 1.67 m- 
is about 6% of that at 0. 7M-). Thus 0. 7M is the preferred wavelength of observation, 
and although there is some ozone absorption, the calculations show a negligible 
effect of ozone variations. 


( 2 ) 

♦These were selected because previous scattering calculations ' were made at 
these wavelengths and because the region \>0. 7 p, is little influenced by ozone 
absorption and Rayleigh scattering. 
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Nor malized Radiance 



Figure 4-1. Normalized upward radiance vs sun angle as function 
of wavelength. 


The angular distribution of the backscattered radiance is not a sensitive 
indicator of aerosol amounts, especially in the region M-= 0.5 to 1.0 where observa- 
tions from a typical polar orbit of 600 nmi altitude would be constrained due to the 
fact that a look angle of M>= 0. 5 is approximately tangent to the earth's surface. 

Also, non-vertical look angles would have a greater probability of having clouds in 
the field of view. 
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The smoothed data in Fig. 4-1 are replotted in Fig. 4-2 to show the relation- 
ship between the upward radiance, normalized to unit incident solar flux, and the 
aerosol content of the atmosphere for various sun angles. A simple linear relation- 
ship is shown to exist between radiance and the mass loading. These straight lines 
are based on only two values of mass loading, but a linear relationship may be 
established by considering the published data of Plass and Kattawar for zero albedo, 
which are shown in Fig. 4-2. 

It is concluded then, for the ideal model conditions that the aerosol content of the 
atmosphere can be determined from radiance measurements from a satellite looking verti- 
cally down at the ocean at a wavelength of 0. 7 p,. For the real atmosphere and ocean certain 
problems arise, as discussed in Section 4. 1.5 so that more than one wavelength of obser- 
vation will be needed. 

4.1.2 Polarization Observation^ 

The polarization of radiation scattered by molecules is greater than that scattered by 
particles. Thus the polarization of radiation scattered by the atmosphere is reduced 
as particles are added to the atmosphere. It is this difference in polarization from 
that of the theoretical Rayleigh atmosphere which must be analyzed to determine 
the aerosol content. The aerosol content cannot be deduced directly due to the 
complexities of the coupling of Rayleigh and Mie multiple scattering and the presence 
of non— Lambertian surfaces, and must be determined by comparing the measurements 

with calculations for model aerosol atmospheres. 

( 2 ) 

Sekera ’ has discussed the theory of this type of observation, and Hariharan^ 
has performed aircraft measurements, making some deductions about the aerosol 

size distribution. However, the models used did not account for multiple scattering. 

( 2 ) 

Plass and Kattawar have made Monte Carlo calculations, including multiple 
scattering,of the polarization of radiation scattered by realistic model atmospheres 
with various aerosol contents at wavelengths from 0.27 m- to 1. 67 m-. The wavelength 
most sensitive to aerosol changes is 0.7M; at shorter wavelengths the Rayleigh 
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Normalized Radiance 



4,2 

Aerosol Mass Loading 10 (J-g/m 


Figure 4-2. Normalized upward radiance vs aerosol mass loading 
as function of sun angle 
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optical thickness is much greater than the aerosol thickness and the molecular 
polarization dominates, and at longer wavelengths the polarization is small with 
corresponding small effects by aerosols. 

As in the case of radiant Intensity observations the effect of aerosols decreases 
as the albedo increases, so that measurements of polarization over low-albedo sur- 
faces are desirable. Unfortunately, as found by Coulson et al. <5) low reflectance 

surfaces are highly polarizing; this would make the modeling and data interrelation 
difficult. 


For an overhead sun, the most sensitive look angle, from the calculations of 
Plass and Kattawar, is given by U=0. 35. However, as pointed out earlier, for 
NIMBUS-type orbital altitudes, observations are limited to it. 0.5 to 1. The 
sensitivity decreases as n increases, i.e. , looking towards the horizon is preferred. 
However this look angle will almost certainly include clouds in the field of view 
which would complicate the data interpretation. 

The results of Plass and Kattawar for a wavelength of 0. 7P are presented 

in Fig. 4-3. The curves illustrate the decrease of sensitivity as the look angle 

approaches the vertical, although there is little difference between u= . 55 and 

.15. It should be noted teat these curves are for zero albedo; there are currently 

no calculations available on tee polarization of the upward radiance over an ocean 

surface (since polarization calculations are considerably more complex than cal- 

ations of radiance only), so the discussions here are limited to calculations for 

ero albedo. Thus the values of polarization and their sensitivity to aerosol content 

are higher than are to be expected in the real atmosphere with an underlying non-zero 
albedo. 

The calculations of Plass and Kattawar have been made mainly for an overhead 
sun, whereas for satellite observations it is of more interest to consider the 
observer looking vertically down with a variable sun angle. Unfortunately a simple 


78 



Polarization 


I ■■■■■! '■■■■■■■■■ Hill ■■■■■■■■■■ ■■■■■■■■■■ 


■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■ 

isssssi* ■■■■■■■■■■■■ ■■■■■■■■■■■■■■■■■■■■ ■■■■■■■■■■■■ 

■■■■■■■■ \m ■■■■■■■■■■ dbbbbbbbbb ■■■■■■■■■■ ■■■■■■■■■■■■ 
■■■■■■■■i ■■■■■■■■■■■■■■■■ ■■■■■■.■■■■■■■■a. ■■■■■■■■■■ 

mi ■■■■■ in ■■■■■iiiii ■■■■■■■■■■■■ 

■■■■■■■■■ t ■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■ ■■■■■■■■■■■■ 

aaaaaaaaav ■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■...■■..■■■. 

■ ■■■■■■■■K «■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■ 
■■■■■■■■■a <■■■■■■■■■■■■■■■■■■■■■■■■■■■■■ ■■■■■■■■■■■& 
■■■■■■■■■a i ■■aaBiaiiiMaaaHMBBaaaaaaBiBmiiUMaa 

& laiaia an BBaaBaaaaa aiaaa uaaBBaaBBBBBBB am 



■■■■BBBBBBBl ■■■BBBBBBBBBBBBBBB BBBBB BBBBB IBBIBBBBBB IB 
BBBBBBBBBBBa 1 BBBBBBB BBBBBflBBBfl BBBBBBBBBB BBBBB BBBBB ■■ 
aaBBB BBBBBBBl BBBBBBB IflBBaaBBBIflBBBBBBBBBIBBBBaaaBBBB 

■ ■BBBBBBaflBflllflBBaaBIBBaBBaBflBIBBBIBRBBa BBBBBBMHM 

■ ■■afl BBBBB BBa> IB BBBBB ■■■■■■■■■■ BBBBB BBBBB BBIBBBBB 
■■■■■■■■■■■■■I IBIBBB BBBBB BBBBB BBBBB BBBBB BBBBBBBB 
BBBBB BBBBB BBBB 1 BBBBB BBBBB BBBBB BBBBB BBBBB BBBBBBBbIM 
IBBBBBBBBBBBim BBBBB BBBBB BBBBBBBBBBBBBBB BBB»B 8 iaaBB Baa 


■BBBB 

BBBB 

BBBB 


S HflBBBBBBBB BBBBB BBBBB BBBBBBBBBB BBBBBBBB 
■bS5bB|BBB BBBBB BBBBB BBBBB BBBBB BBBBBBBM 
■BBBBfliBBB BBBBB BBBBB BBBBB I BBBB BBBBBBBBBlI 
[flfl BBBBBBBB BBBBBBBBBB BBBBBBBBBB BBBBBBBBBI I 
BBBBB BBBBB BBBBBBBBBB BBBBB BBBBB BBBBBBBBBI I 
BBBBBBBBBB BBBBBBBBBB BBBBB BBBBB BBBBBBBBBI I 
BBBBB BBBBB BBBBBBBB BBBBBBB BBBBB BBBBBBBBBI I 
BBBBB BBBBB BBBB BBBB BBBBBBB BBBBB BBBBBBBBBI ( 
BBBBBBBBBB BBBBBBBBBB BBBBBBBBBB BBBBBBBBBI I 
BBBBBBBBBB BBBB BBBB BBBBBBB BBBBB BBBBBBBBBI I 
BBBBBBBBBB BBBB BBBBBBB BBBB BBBBB BBBBBBBBBI I 
BBBBB BBBBB BBBBBBBBBB BBBBBBBBBB BBBBBBBBBI I 
BBBBB BBBBB B BBBB BBBBB BBBBB BBBBB BBBBBBBBBI I 
BBBBBBBBBB BBBBBBBBBB BBBBBBBBBB BBBBB BBBBI I 
BBBBB BBBBB BBBBB BBBBB BBBBBBBBBB BBBBBBBBBI I 
BBBBB BBBBB BBBBB BBBBB BBBBB BBBBB BBBBB BBBBI I 
....... BBBBBBBBBB BBBBBBBBBB BBBBBBBBBI 

BBBBB BBBBBBBBBB BBBBBBBBBB BBBBB BBBBB BBBBBBBBBB BBBBBBBBBI 


BBBBBBBBBBBBBBB BBBBBBBBBB BBBBBBBBBB BBBBB BBBBB BBBBBBBBBI I 


BBBBBBBBBB 


I BBBBB BBBBB BBflflflBk'BBB BBBBBBBBBB BflflBBBBBBBBBBBBBBfliMV 
BBBBB BBBBB BBBBBBB BBBBBBBB BBBBB BBBBBBBBBB BBBBB BBBBBBB 

52555 25552 25555 55i.'!5 55555 5555! 55555 5555! 55515 55555 5! 


■ BBBB BBBBBBBBBB BBBBBBBBBB B.l 


I BBBBBBBBBB BBBBBBBBBB Bk BBBBB 
I BBBBBBBBBB BBBBBBBBBB BBBBBBB 


^HBBaBBBBB BBBBBBBBBB 
BIB BBBBBBBBBB BBBBBBBBBB B 
BBBBBBBB BBBBB BBBBBBBBBB B 




BBBBBBBBBB BBBBB BBBBB BBBBB BBB.'B BBBBB BBBBB B 
BBBBBBBBBB BBBBBBBBBB BBBBB BBBB.-' BBBBBBBBBB! 
BBBBBBBBBB BBBBBBBBBB BBBBBBBBBB ^BBBBBBBBBB 
■BBBB BBBBB BBBBBBBBBB BBBBBBBBBB B^BBBBBBBBS 


BBBBBBBBBB 

BBBBBBBBBB 

BBBBBBBBBB 


■■■BBBBBBBB 


BBBBB BBBBBBBBBB BBBBBBBBBB B« 
BBBBB BBBBBBBBBB BBBBBBBBBB ■■ 
B BBBBB BBBBBBBBBB BBBBBBBBBB bB 


■BBBB ■■■■■■■■■■ ■■■■■■■■■■■■ 
BBBBI BBBBBBBBBI BBBBBBBBBB ■■ 
■BBBBB BBBBB BBBBB BBBBBBBBBB ■ ■ 

■■■■■■■■■■■■■■■■■■■■■■■■■■■I 

■■■■■■■■■■■■■■■■■■■■■■■■■■■I 

■BBBBB ■■■■■■■■■■■■■■■■■■■■■ ■ 


■ ■■■■■■■I I 
^^■■■BB BBBBI 

■ ■■■■■■■■■■□I 

■ ■■■■■■■■■■■I 1 


BBBBB ■■■■■■■■■■ 
!■■■■■■■■■■■■■■ 
■■■■■■■■■■■■■■i j 


■BBBB BBBBBBB BBl I 


■ ■■■ 

■ ■ ■ b^H 

1BBBB ■■■■■■■■■■■■■■■■■■■■■■■■■ ■■■■BB.-’BBB ■■■■■■■■■■■■ 
L'BB.B BBBBB ■■■■■■■■■■ ■■■■■■■■■■■■■BBBBBBB ■■■■■■■■■■■■ 
B. 1 BBB BBBBB BBBBB BBBBB BBBBBBBBBI BBBBB BBBk'B BBBBB BBBBBBB 

■ k'BBBBBBBB BBBBB ■■■■■■■■■■ ■■■■■■■■■■■■■ BkJ BBBBB ■■■■■■■ 

■■BBIBBBBB ■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■"■■■■■■! 

BBBBB. BBBBBBBBBBBBBBB ■■■■■■■■■■ ■■■■■■■■■■ ■■■■■■■*”■■■ 
BBBBB BBBBB ■■■■■■■■■■■■■■■■■■■■ ■■■■■■■■■■ ■■■BBBBBBB BBl 

■ ■■■■■BBBB B. 1 BBBBBBBB BBBBBBBBBB BBBBBBBBBB ■■■■■■■■BBBBB 

:sss:::ss::ssss^::s::::ssss:s::::s::::s::::::ss:::s| 



■■■■■■■■■■ ■■■■■■■••■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■ ■■ 
BBB ■■■■■■■ ■ ■■■■■■■!>■ ■■■■■■■■■■ ■■■■■■■■■■■■■■■■■■■■■■ 

■ B. ■■■■■■■... ■ ■■■B.k'BB ■■■■■■■■■■■■■ ■■.■■■■HMaM 

■ ■.....■■■ BBBBBBB.. ■ .'BBBBBBBBB BBBBB BBBBBBB 

J 0 ggg e a m n a_g a rs 5 b a BBaa^Hama a biibb 5 555 b asa 


BBBB BBBBB ■■■■■■■■■■■■■■■■■■■■ ■■■■■■■■■■ I 
■■■■■■■■■■■■■■■■■KBB BBBBB BBBBB BBBBBBBBBB BBBBB BBBBBBBBBI I 

■ ■■■■■■■■■BBBB ■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■ ■■■■■■■■■■ I 
..■■■...■■■■■■....■.....■■■■■■■■■■■...■■■■■.■■■■■■■■■I I 
BBBB ■■■■■■■■■■ ■■■■■■■■■■ ■■■■■■■■■■ ■■■■■■...■ B.B.B BBBBI I 

■ ■■■■■■■■■BBBB ■■■■■■■■■■ ■■■■■■■■■■ ■■■■■■■■■■ ■■■■■■■■■! I 

■ ■■■■■■■■■■■■■■■■■■■BBBB ■■■■■■■■■■ ■■■■■■■■■■ ■■■■■■■■■! I 
BBBB ■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■ ■■■■■■■■■! I 
BBBB ■■■■■■■■■■■■■■■■■■■■ ■■■■■■...■ BBBBBBBBBI I 
...B ■■■■■■.... B. ........ ....... BB. BBBBB BB..B BBBBBB..BI I 

............. BB. ...... BB B....B...I I 

...B ■■..■■■■■■ a....... .B........E. ■...■...■■ ■..■■■■■■■ I 

■ ..■.....■....■.....■■.■■.■■■■■■■■■■■..■...■■..■■■■■■I I 

■ ■■■■■.■■■BBBB ■■■■■■■.■■■■■■■■■■.■ ■■■■■■■■■■ ■■■■■■■■■! I 
|BHHH||BI ■■■■■ ■■■■■■■■■■ BBBBB BBBBB ■■■■■■■■■■ ■■■■■■■■■! I 

B... .......... .......... .......... ......... ■ ....... BBl I 

■ ■■■■■■■■■■■■■ ■■■BMMMMMMUMMMMMMMMMMMMmJ 
.'■■■■■■■■■■•BBBBB 

■ • ■ ■ ■ ■ a ■ ■ ■ ■ ■ ■ ■ a ■ 

■ ■■.' .............................. B........B ......... I I 

.... k'BBBBBBBBB .................... .......... .........I I 

... ■■■.''■■■■■BB ■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■ ■■■■■■■■■! I 

BBBB BBBw •■■■■■■■■■■■■■■■■■■■■■■■■■ BBBBB BBBBB ■■■■■■■■■! I 

..■.■.....'*■........■....■...................■....■■■■1 I 

.... ■ ■■■■ ■k.'IBB B...B....B ■■■■■■■■■■ BB... B.B.B BBBBB. BBBI I 
...■■...■B.B.?. I 

.■■■■...■■...■ •■■■■■■■■■ ■........■■■..■■■■■■ ■■■■■■■■■! I 

.BBB. ..■■■■■■■ ■^■BBBBBBB BBBBB BBBBB ■■■■■■■■■■ ■■■■■■■■■■ I 

■ ■■■BB. ■■■.... BBBBBBBBBB BBBBB BBBBB BBBBBBBBBB BBBBB BBBBI I 

■ ■■■■■ ..■■■■■■ ■■■.•**..■■ ■■■■■■■■■■ ■...■■...■ ■■■■■■■■■! I 

BBBI ■ ...■■...■■■■■■■>-'•.■ BBBBB BBBBB ■■■■■■■■■■ ■■■■■■■■■! I 

■ ■■■■■■■■■■■■■ ■■■■■■■»..'■■■■■■■■■■■■■■■■■■■■■ BBBBB BBBBI I 

■ ■■■ ■■■■■■■■■■ ■■■■■ BBBBk. •BBBB BBBBB BBBBBBBBBB BBBBBBBBBI I 
.BBB ■■■■■BBBBBBBBBBBBBBB ■».•■■■■■■■ B.B.B BBBBB BBBBBBBBBI I 

■■■■ ■■■■■■■■■■ BaBBBB.BBBBBfe>IBBBBBIB|B|BBB|BBBBBBBBB||J 

■BBBBB ■■■■■■■■■■ ■■■■■■■■■■ BBBBh'BBBBn 


■■■■■■BBBBBBBB ■■■■■■■■■■ ■■■■■■■■■■ 9 

B. . BB. BB BBBBB BBBB. 4 MHHRMIIIIM 

BBBB ■■■■■■■■■■ ■■■■■■■■■■ BBBBBBBBBB BBBBBBBBBB ■■■■■■■■■ 




■ ■B.B...B B.l 


. ......B....I 


BB... ■..■■! 
■■...■■■■El 


B.B.. ■■••■■■■■■ ■ 
B.B.B B...w ■■■■■■ 
BBBBB ■■■■■•:'■■■■ 
■■■■■■■■■■■■*•■■ 
BBBB. B.B.B B.BBlwB 
. ■■■■■... ■■■..■■ 


...B ■■■■■■■■■■ ■...■■■■■■ B.......B 

■BBS ■■■■■■■■■■■■■■■■■■■■ ■■■■■■... 
..■.■■...■■■■■■■■.■■■■.. ■■■■■■■■■ 
.... .......... .......... ......... I 

IB... ■...■..... ■...■■■.■■ ......... 1 

....■■.■■....■■.....■.■.■........I 

■^•■■■.■BBBBBBB ■■■■■■■.. ■■■■■■■■■■! 


■■...■..BBl I 


■ B.B.B B. BBl I 

■ ■■■.■■...1 I 

I ■■■■■BBBBI I 


>■ ■■■■■■■■■■ ■■■■■■■■■! I 



^...BB.B.B B........B ......... ■ ■....■.... ■..■■■■..■■ 

a ■■■■■■■■■ ■■■■■■■■■■ ■■■■■■■■■■ ■■■■■■■■■■ ■■■■■■■■■■ ■ 

■ k’B. ■■■■■■ ■■■■■■■■■■ .......B.B ■....■.... ■..■■■■■■■ ■ 

■B •.■■...■■..■■■■■.■ a....... MB ■■■■■■■■■■ ■■■■■■■■■■■ 

■ ■^ VBBBB|| 1 U||||| 11 | 1 |BBBIBII ■■■■■■■■■■■■■■■■■■■■■ 

■ ■■>'■ ■■■^^^^^^^^^^■■■BBBBB ■■■■■■■■■■ ■■■■■■■■■■■ 

■ ■■■•■■■^^^^^^H"li^HBBBBBBB ■■■■■■■■■■ ■■■■■■■■■■ ■ 

■ ■■■■•■■SV^^^^H*£^Hb ■■■■■■ ■■■■■■■■■■■■■■■ ■■■■■■ 

■ ■■■Bk«aH«'^£-^^^.j.^HBBBBBBB BBBBB ■■■■■■■■■■ ■■■■■■ 

■■■■■■»> 

BBBBB BB. ^^^^^^^^^^HBBBBBBB BBBBBBBBBB BBBBB BBBBB B 

■b. ■■■■■. ' rmnvvnnnBBBB. .■■■■■.■...■■■■.■■..■■■ 

■■■■■■■■■k BB. ■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■ B 

■■■■■■■■■■ k •■■■■■■■■ B...B....B BB... B.B.B ■■■..■.... ■ 


BBBB ■....B.Bk ' B...B..... B........B ■....■...■ B...B...BI I 

■ ■■■■■■BBBBBBB BBBBB ■....■■■■■■■■■■■■■■■ ■■■■■■■■■■ I 

■ B.B ■■■■■■■■■■ ■ ■••?■■■■■■ B...B..BBB ■■■■■■■■■■ ■■■■■■■■■> I 

.... ■■■■■....■ ■■■..^•... B... ...... B........B BMMfififiBMB 

■■■■ BBBBBBBBBBBBBBB ■Ik^lIBBBIBBBia ■■■■■■■■■■ 

■ ■■■■■■■■■■■■■ ■■■■■■■■■■ •:•■■■■■■■■ ■■■■■■■■■■ 

■ BBB ■■■■■■■■■■ ■■■■■■■■■■ BBk^B ■ ■■■■■ 

■ ■■■ ■■■■■■■■■■ ■■■■■■■■■■ ■■■■*-''■■■ ■ 

■ ■■■BBBBBBBBBB ■■■■■■■■■■ ■ ■■■■■■*.'' • Bl 
■■■■■■■■■■■■■■ ■■■■■■■■■■ ■■■■■■■■■• J 
■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■I 

BBBB ■■■■■■■■■■ ■■■■■■■■■■ ■■■■■■■■■■ ■< 

■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■BBBB' 

■■■■■■■■■■■■■a ■■■■■■■■■■■■■■■■■■■■ bH 

■■■■ ■■■■■■■■■■ ■■■■■■■■■■■■■■■■■■■■■■■ 

■■■■ ■■■■■■■■■■ ■■■■■■■■■■ ■■■■■■■■■■ IBB 
■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■ 

■■■■■■■■■■■■■■ ■■■■■■■■■■■■■■■■■■■■■■■ 


■ ■■■■■■■■■■■■■■■B.B.' «■■■■■■■■■■■■■■■■■■■■■ 
■...■■...■■...■■.■BB k^a.B ■■■■■■■■■■ ■■■■■■■ 
■B. ■■■■■■■ ■...■■■■■■ BBk. '*■■■■■■ ■...■■....■■ 
■■■■■■■■■■ ■■■■■■■■■■ ■■■■•'■■■■ ■■■■■ ■■■■■ BBl 
■BBBB B.B.B ■■■■■■■■■■ ■■■■■■■^•B ■■■■■■■■■■■■ 

■ !•■■■■ ■ ■«■•■ 

■■■■■■■■■■ ■■■■■■■■■■■■■■■■■■■■ BhC^BB. 


.. ..■■■■■ 

.......BB 

■■■■■■■■■ 

■■■■■■■■■ 


BB. ■■...■• I 


■■■■■■■ 

■■■■■■■ 

.■...■■HHHMB 

■ ■■■■■■■.B.B. BBl I 


■■■■■■■■■■■■■■■■■■■■■■■■■■■■■BBBBBBBB ■■■■■■ 

■ ■■■■■■■BB. ■■■■■■■■ ■■■■■■■■■■ B.B.. BB. ■ B.B.B M _ . 

B. ....... BBBBB. B.BBBB. ■■■■■■■ ■....■........■iiiiiigiiil 

... ■■■■■■ BB. BB... B.B BBB.. BBBB. BHnHHMHMBBd 

■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■ 

■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■I 


■ ■■....■...■■..■■■■■■I I 


B.B.B B.B.B ■■■■■*-Z' , "B 

■■■■■■■■■■■..■.■...■■..■■■.■.■■...■■.■.■■■.■■■■■■■a 

BBBBBBBBBB ■■■■■■■■■■ ■■■■■■■■■■ ■■■■■■■■■■■■■■■■■■■■■ 


■B.B....B B...BBB. 
■..■■■B.B IBB. BBB. 
■ ■■■■BBBB BBBBBBBB 


Aerosol Content 


Figure 4-3. Polarization vs aerosol content as function of 
' nadir look angle 
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reciprocity relationship does not exist for polarization^ 

However some inference may be made. If a relationship for polarization similar to 
that for intensity existed, viz. , 

W(0, W = I (0,H o ,«> o ; !>,<») 

where I (0,^,0;^^) is the radiance at the observation angle specified by M-, 0 due 
to an incoming solar beam specified by M- ,0 , then the values of polarization and the 
sensitivity to aerosol content shown in Fig. 4-3 would be reduced; if the polarization 
were unchanged by interchanging the observer and source, then the sensitivity would 
be unchanged. Thus it appears that calculations for vertical observations will not 
indicate more sensitivity than that shown in Fig. 4-3. 

4.1.3 Comparison of Radiance and Polarization Observations 
A comparison of the accuracies of the radiance and polarization observations of 
aerosols for the model atmospheres may be made by considering the sensitivity curves 
shown in Figure 4-2 and Figure 4-3. It is reasonable to assume that the radiance 
can be measured to ±1%, so that from Figure 4-2 it is found that the error in aerosol 
content is about ±1. 3%. The polarization of backscattered radiation can probably 
be measured to an accuracy of ±1% in the range 5-75%^ ^ so that the error in near- 
normal values of aerosol content is about ±7. 5% for a look angle of n = . 55. It 
should be noted that this large error is somewhat underestimated, since it is for zero 
albedo; in addition, it will probably be even larger if the observer and sun are inter- 
changed to consider more realistic conditions for satellite observations. 

It should also be remembered that calculations of polarization for various 
vertical distributions of aerosols have not been performed, so it is not known whether 
polarization, like radiance, is independent of the distribution. Of course, the error 
estimates for both techniques assume that the relationships based on the calculations 
of Plass and Kattawar in Fig. 4-2 and 4-3 are correct. 
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It is concluded from the preceding discussions, based on model calculations, 
that satellite measurements of radiance are preferred to measurements of polariza- 
tion to determine aerosol content. This is in contradiction to the often repeated 
statement that polarization measurements offer the best means of monitoring aerosol 
content. Such a conclusion may have been reached on the basis of ground-based 
observations of sky polarization. However, the results of Plass and Kattawar clearly 
show that the backscattered polarization is significantly less sensitive to aerosol 
changes, particularly at look angles close to the horizon. It is possible that 
polarization observations would provide information on the aerosol size distribution 
and vertical distribution; however this is not apparent from the published literature. 

4-1.4 Observations over Land or Clouds 

The previous discussions showed that the aerosol content may be determined with an 
underlying surface of known low albedo, such as the ocean. With higher albedos 
such as for land (-0.18) and clouds (-0.6) it is difficult to determine the aerosol 
content accurately without precise knowledge of the surface reflectivity. 

Plass and Kattawar ^ ^ calculated the backscattered radiance (at 0. 74 ) as a 
function of albedo for the Kondratiev et al. aerosol vertical distribution. (The 
aerosol optical thickness for this model is slightly less than that for the Elterman 
distribution). These results are shown in Fig. 4-4 together with curves for different 
aerosol contents based on later calculations by Plass and Kattawar (2) . 

These latter calculations, based on the Elterman distribution, were made for 
albedos of A=0 and 1. 0, so that the curve-shapes were made to follow the 1968 
calculations for intermediate values of albedo. The different aerosol contents are 
noted on the curves as N (the normal Elterman distribution), 3N, N/3 and 0 
(Rayleigh atmosphere). It is clear that as the albedo increases, the sensitivity of 
the radiance to aerosol content decreases. 
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Reflected Radiance 



Figure 4-4. Reflected radiance vs albedo for different aerosol 
contents 
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A cross-plot of the data in Fig. 4-4 is given in Fig. 4-5. , and shows that large 
errors are possible in measuring the aerosol content if the albedo is not known 
precisely. For example, over land, a normal aerosol content above a true albedo of 
0. 2 would be interpreted to be about 2. 7 times greater if the albedo were incorrectly 
estimated to be only slightly lower at 0.18. 

The polarization calculations of Plass and Kattawar (1970) are shown as a 
function of aerosol content and albedo in Figs. 4-6 and 4-7. In Fig. 4-6, only the 
data for the normal aerosol content were calculated for several values of albedo; the 
other curves are based on values at A =0 and 1. 0 and the shape of the curve N. It is 
seen that at albedos greater than 0. 4 (e. g. , clouds or snow) there is essentially no 
sensitivity to aerosol content. At the same time, the polarization does not change 
much with albedo so that no information is obtained on the albedo from a polarization 
measurement. Using the same example as above for the radiance measurement, it 
is seen that the estimated aerosol content would be only 1. 3 times greater if the 
albedo were incorrectly estimated to be 0. 18. Thus it appears that a polarization 
measurement would give a better estimate of the aerosol content. However, when 
the measurement errors in the polarization (±0.01) and radiance (±1%) are considered 
in addition to the albedo error in the above examples, then the techniques are 
comparable, both giving the aerosol content within about a factor of 3. 

From the above discussions, measurement of aerosol content above high 
albedo clouds will be difficult. However it is probably generally true that above clouds 
the aerosol content will be close to that of Elterman's clean air model. Over land, 
where the albedo is lower (~0. 18) , the measurement of absolute aerosol content is 
not accurate enough due to albedo errors; however, it will be possible to make 
accurate relative measurements over a period of time above any particular area. 

From such a series of measurements it'will be possible to determine the surface 
albedo by assuming that the lowest radiance value corresponds to the Elterman 
clean air condition. 
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Figure 4-5. Reflected radiance vs aerosol content as function of albedo 
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Polarization 



Albedo 


Figure 4-6. Polarization vs albedo for different aerosol contents 
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Figure 4-7. Polarization vs aerosol content as function of albedo 
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It may be possible to improve the aerosol determination over high albedo 
surfaces by making radiance measurements as a function of wavelength or by making 
simultaneous radiance and polarization measurements. The investigations ^°^using 
ERTS-A data and ground-truth data will be invaluable in resolving these problems. 

4.1.5 Problem Areas 

The proposed method of measuring aerosols from radiance observations is based on 
theoretical calculations which use a model atmosphere, a model aerosol distribution, 
and assume a smooth specular-reflecting water surface. In practice, of course, these 
model conditions are never realized so that deviations from the theoretical relation- 
ships are to be expected. We have proposed to investigate these deviations under an 
ERTS program* \ which will include overflights by aircraft to obtain additional 
ground-truth data. We believe that only the combination of a theoretical study and 
field measurements will make a solution of some of the problems possible. The 
anticipated problems and their possible solutions are as follows: 

4. 1. 5. 1 Sun Glitter: If the ocean were perfectly smooth as assumed in the calculations 
an image of the sun would be seen at the specular reflection angle, and the only 
upwelling surface radiation observable at other look angles from space would be the 
diffuse sky radiation reflected from the ocean surface and the radiation scattered 
up from below the ocean surface. As the smooth ocean surface is increasingly 
disturbed, a glitter pattern becomes increasingly larger about the specular point. 

At sun zenith angles greater than about 30° the glitter effect has been considered 
negligible (except for very rough seas) at the nadir point. However, recent 
measurements by Hovis (private communication from R. S. Fraser , NASA, 

Goddard Space Flight Center ), suggests that this assumption is not correct, so 
that the ocean surface radiance at the nadir is not known accurately. 

This problem might be overcome by making observations at two wavelengths, 
assuming that the spectral variation of the surface radiance is known. The choice 
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of wavelengths must be made carefully since the spectral distribution of the radiance 
does vary due to ocean properties such as chlorophyll content, suspended matter and 
depth. 

4. 1.5. 2 Aerosol Properties. In order to make model calculations, ass um ptions 
must be made about the aerosol properties, such as size distributions, vertical 
distribution, sphericity of the particles and the refractive index. All of these 
assumptions affect the phase function to be used in the multiple scattering calcula- 
tions. The variability of each of these parameters in the real atmosphere should be 
investigated. Some work along these lines has been reported by Rosenberg 

who suggests that the Junge aerosol model is not applicable to the atmosphere - in the 
case of humid atmospheres, due to condensation on the particles, and for dry 
atmospheres because the particles are non-spherical. 

4.1.5. 3 Non- Lambertian Surfaces. Most calculations of scattering properties of the 
atmosphere use an underlying Lambertian surface. However, the calculations for 
the ocean treat the surface correctly as a specular reflector for both the direct and 
sky radiation, and include the upwelling radiation scattered by the water from beneath 
the surface. Raschke uses a two-layer model of the atmosphere -ocean system, 
and a rough ocean surface (equivalent to a wind speed of about 1.4 m sec to calculate 
the radiance outside the atmosphere, but does not show a quantitative comparison 
with calculations for a smooth surface. 

4. 1. 5. 4 Surface Reflectance Gradient. The calculations of Plass and Kattawar 
assumed an underlying surface of constant reflectivity extending to infinity. However 
in the case of two adjacent surfaces with different reflectivities, radiation reflected 
from one surface will be scattered into the atmosphere above the other surface so 
that the apparent radiance above that surface is increased over the calculated theoreti- 
cal value. This is probably a second order effect if the radiance is measured at 
points well removed from the boundary of the two surfaces. This will be verified 
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by comparing the theoretical calculations with aircraft measurements 9"^. 

4. 2 SUN-ORIENTED OBSERVATIONS OF AEROSOLS 

Satellite measurements of aerosols in the stratosphere are well suited to the 

occulta tion method, (illustrated in Figure 4-8), in which observations are made of 

the sun's radiation as the satellite enters or leaves the earth's shadow. From the 

observed attenuation along successive light paths, and allowing for refraction, the 

vertical distribution may be determined, if the aerosols are assumed to be uniformly 

distributed in the horizontal. These observations are generally restricted to the 

stratosphere and higher due to the presence of clouds and high attenuation (scattering 

and absorption) in the troposphere; these problems are discussed in some detail below. 

The technique has not yet been applied to aerosol measurements, although 
(13) 

Pepin v ’ has satellite measurements scheduled. The only satellite occultation 

measurement to date appears to be that of ozone by Rawcliffe et al. and 
(15) 

Miller and Stewart v . Their measurements used the full disc of the sun (0. 5° 
angular diameter), so that the vertical resolution in the deduced ozone was only about 
16 km. Improvement of this height resolution would require the use of a sun-pointing 
device to observe a small area of the sun. To obtain a 2 km resolution, one -eighth 
of the disc's vertical diameter would be used, requiring extreme pointing accuracy, 
which will be very difficult to achieve on a satellite. 

4.2.1 Cloud Interference 

An obvious problem in using occultation methods to probe the lower atmosphere is 
the presence of clouds. Most clouds are confined to the troposphere, the height of 
which varies from about 17 km in the tropics to below 10 km in polar regions. If it 
is assumed that the clouds reach up to 10 km, it may be calculated that a cloud-free 
region nearly 700 km wide is required so that the line-of-sight from a 300 km orbit 
to the sun tangent to the earth's surface is unobstructed. Such a large clear area will 
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Figure 4-8. Schematic of oc 
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experiment 



be very rare; even if the atmosphere were cloud-free, the strong attenuation (scatter- 
ing and absorption) of the radiation through the lower atmosphere would preclude 
occultation probing at these low altitudes (see below). 

Even if the observations are confined to the stratosphere, the problem of clouds 
is not entirely eliminated. Clouds do occur in the stratosphere. In addition to the 
occasional penetration of the tropopause by thunderheads, there are, at higher 
altitudes, nacreous (or mother-of-pearl) clouds and noctilucent clouds. Noctilucent 
clouds occur at about 80 km and are generally observed at high latitudes in the 
summer. The nacreous clouds are formed in the 20-30 km region and observed in 
middle latitudes in the winter. This latter altitude region is exactly where the SST 
(supersonic transport) commercial aircraft will be flying, so with the addition of 
particulates and water vapor from the SST exhaust products we might expect more 
observations of nacreous clouds in future years. These stratospheric clouds present 
a problem to occultation measurements since they will attenuate the sun's radiation 
by an unknown amount resulting in misinterpretation of the observations. However, 
if the interest is in the global loading of aerosols in the stratosphere, the data may 
be averaged over all latitudes to minimize errors due to stratospheric clouds. 

4.2.2 Atmospheric Attenuation 

The attenuation of the sun's radiation varies as a function of wavelength, and as a 
function of the grazing altitude of the radiation. Both gaseous absorption and mole- 
cular scattering must be considered. In the ultraviolet and visible region, ozone 
absorption is important, and in the infrared beyond about 0. 8 M- absorption by other 
atmospheric gases, particularly water vapor, is important. Rayleigh scattering is 
negligible in the infrared, but is a significant factor in considering shorter wave- 
lengths for occultation. For aerosol observations an atmospheric "window" must be 
chosen. 
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Calculations of the sensitivity of various wavelengths to changes in aerosol 
content may be made based on the 1968 model atmosphere of Elterman. Fig. 4-9 
and Fig. 4-10 show results of simplified calculations for occulatation grazing 
heights of 10 km and 15 km respectively. The calculations are simplified to the 
extent that the sun is considered as a point source, so that the variation of transmission 
of the light rays from different parts of the sun is neglected. In addition, refraction 
of the light rays in the atmosphere is ignored. 

Thus the transmission t is given by 

t = exp(-2t(au) Ch(90°)) 


where 

t( (jo) is the wave number-dependent vertical optical thickness above the 
grazing altitude and is the sum of contributions from Rayleigh 
scattering and aerosol scattering. Gaseous absorption is negligible 
in atmospheric windows at high altitudes. 

Ch(90° ) =35, is the Chapman function for zenith angle 90° . 

Ch(9) replaces the usual secant 6 when 9 is large. 

The factor 2 accounts for the two equal optical paths BD and BE in Fig. 4-8. 

The results show, as expected, that the shorter wavelengths are more 
sensitive to aerosol variations than the infrared region. Observations of attenuation 
at the shorter wavelengths would have to take into account the contribution of 
Rayleigh scattering, which is well known. In the infrared, the Rayleigh scattering 
is negligible, and, at the higher altitudes above about 12 km, the water vapor 
continuum, absorption is negligible. Thus, no corrections to the observed data 
would be required at high altitudes in infrared window regions. However, the 
increased sensitivity, the well-established Rayleigh correction, the higher sun 
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Aerosol Content 


Figure 4-9 Transmission vs aerosol content as function of 
wavelength (grazing altitude 10 km) 
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Transmission 



0 N 2N 

Aerosol Content 

Figure 4-10 Transmission vs aerosol content as function of 
wavelength (grazing altitude 15 km) 
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radiance, and more sensitive detectors at shorter wavelengths would point to choos- 
ing the visible or ultraviolet region for aerosol observations by occultation. 

It is seen from Fig. 4-9 and 4-10 that at shorter wavelengths the transmission 
changes by 10-20% for a 10% variation in the normal aerosol content. Thus the 
sun's spectral radiance with and without atmospheric attenuation must be measured 
very accurately to monitor stratospheric aerosol changes due to SST operation. It 

/I \ 

has been estimated that in the 1985-1990 period, assuming operation of 500 SST's, 
the aerosol concentration at SST altitudes will increase by about 10% on the average. 
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5 


INVESTIGATION OF RGFC INSTRUMENT PARAMETERS 

(Task 4) 


5.1 OPTIMIZATION PROCEDURE 

The basic POLAYER program is used to optimize the instrument cell parameters of 
the RGFC instrument. For a given condition in the atmosphere, i.e. , E ((D) at the 
entrance of the instrument, the parameters c l (concentration of pollutant gas)x 
cell length in cm),p T (total pressure in atm, using nitrogen as foreign gas broadener) 
and T (temperature in degrees Kelvin of the gas in the cell) are varied until a 
maximum sensitivity is obtained. These parameters are related through the instru- 
ment cell transmission, t g (<d), and the optical thickness, u by 


V»> = 


exp 


-k ((d) u 


(5-1) 


where k(uj) is the monochromatic absorption coefficient in cm ^atm \ and u 
is the optical thickness, given by 


u = p T c l cm atm 


(5-2) 


The absorption coefficient for a spectral line with Lorentz shape 

Sq'/tt 


k((D ) = 


((D-(D^)^ + 


(5-3) 


-1 -2 

where S is the line strength in atm cm , a is the line half-width at one-half of 
the line strength. The line half-width a is related to pressure, temperature and 
foreign gas and self broadening through 
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(5-4) 


a = a p y~ZT 
o e v T 

where °' o is the line half-width at 1 atm and standard temperature T for nitrogen 
broadening, T is the temperature and P g is the equivalent pressure, given by 


P B N +P N 
P T [c (B-l) + lj 


(5-6) 


where p and ?N are the partial pressure of the pollutant gas and nitrogen gas, 
respectively, and ^is the ratio of self-broadening to nitrogen-broadening efficiency. 


For a given P T and T , an optimum sensitivity of the instrument exists for a 
particular value of the product c A. Calculations of AV were performed by varying 
the CO concentration in the instrument cell and maintaining the other instrument 
parameters, l = l cm, p T = 1 atm, and T = 288. 16 e K for 0. 1 ppm CO uniformly 
distributed through the troposphere, the earth's surface temperature equal to 288. 16° K 
and the standard US atmospheric temperature profile.* The instrument response 
function C(u)) and C V) are assumed to be equal and are given by the curve shown in 
Figure 5-1. The results, shown in Figure 5-2 indicate a maximum sensitivity is 
obtained for C qq = 0.35; the corresponding is = 0.832. 


Additional calculations were made using the optimum value for c (=0. 35), 
cf. Figure 5-2, and the non-optimum values for c cQ = 0.05 and 1.0. The calcula- 
tion was made for CO concentrations ranging from 0 to 0.4 ppm uniformly distributed 
through the atmosphere. The remaining instrument parameters and earth surface 
temperature and atmospheric temperature profile were the same as those indicated 

in Figure 5-2. The results shown in Figure 5-3 indicate that a value of c = 0 35 

CO 

IS an optimum va lue over the atmospheric CO concentration range from 0 to 0.4 ppm. 


This temperature was selected to be equal to the 
standard atmosphere. 


surface temperature of the 


U. 


S. 
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Normalized Filter Transmission 



Figure 5-1. Normalized filter function used in the model 
calculations for CO, corresponding to OCLI 
filter 14-8806-760. 
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Figure 5-2. Optimization of instrument response as a function of CO concentration in 

instrument cell. Instrument parameters: ■£= 1 cm, p_ = 1 atm, = 288. 16° K, 
c(tu) = c (uj), given by Fig. 5-1. Atmospheric parameters: CO concentration 
in the troposphere = 0.1 ppm, ground temperature = 288.16°K, lapse rate = 
-6.5°C/km. Arrow points to maximum response, i.e., CO concentration =* 
0.35, corresponding to a cell transmission of 0.83. 
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Figure 5-3. Sensitivity of instrument as a function of CO concentration 
in the troposphere with three CO concentrations in instru- 
ment cell as a parameter: 0. 05, 0. 35 and 1. 00. Other instru- 
ment and atmospheric parameters same as in Fig. 5-2. 
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We have assrmed that variations in the temperature of the correlating gas in the 
instrument cell, T , will have only small effects on the observed signal output. To 
verify this assumption Run B-09 was re-run but T was set at 270, 300 and 330° K. 

The results of these calculations are shown in Table 5-1 for three different atmospheric 
CO concentrations. 

Table 5-1. Effect of Instrument Cell Temperature Variations 


t g ,°k 

(C0, atm' ppm 

7 2 

10 DV ,W/cm -sr 
c 

% Error 

300 

0.025 

1.7510 




0.100 

4.4431 

— 


0.400 

8.0675 

— 

330 

0.025 

1.7587 

-0.44 


0.100 

4.4319 

0.25 


0.400 

7.9179 

0.19 

270 

0.025 

1.7350 

0.91 


0.100 

4.4330 

0.23 


0.400 

8.1840 

-1.44 


The per cent error 4s calculated according to 

DV (T = 300° K) - DV (T ) 


error = 


( 


DV(300° K) 


x 100^ 
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The calculated values of DV , where DV = DV(corrected) = DV-DV for(CO) =0, 

c c atm 

indeed indicate that variations in T have only a small effect on DV . Since we 

c 

anticipate controlling T = 300° K to within at least ± 1°C, the maximum effect due 
to variations in T will be about . 05 per cent. 

The influence of the total pressure in the instrument cell was studied in computer 
run B20, in which the optical thickness (c l) was kept constant. The results of DV 
as a function of total pressure in the cell with CO concentration in the atmosphere 
as a parameter are shown in Figure 5-4. The maximum sensitivity indicated by 
arrows is not very strongly dependent upon the total pressure in the cell. 


5. 2 LABORATORY CALIBRATION PROCEDURE 

In this section, we investigate whether a laboratory calibration of the instrument 
response using a calibration cell which represents a one-slab atmosphere can be 
applied to pollutant measurements in the real, multi -slab atmosphere. 

The signal difference is given by 

AV = A vJe (id) [R (u>) t q ((D) - R>) t a ] do (5-6) 

Aid 

where Av is the signal output of the instrument (volts), AO is the etendue of the 
2 / 

instrument (cm -sr), and R(uo) and R (u>) are the responsivities of the two cells in 
the instrument (volts/watt). When the instrument is calibrated, E(co) is replaced by 
N (T,<d)t, where N (T,to) is the radiance emitted by a calibration source and t is the 
transmission through a calibration sample of the gas of interest. 

Setting 


R^a 

R(oj) 



(UJ) « T 


A 


(5-7) 


103 







we have 


AV = 


R o f E (UJ) ^ " T A da> 

Auj 


( 5 - 8 ) 


where is the overall instrument responsivity over Acju. Thus, AV is proportional 

to R . 
o 

Similarly, when the instrument is calibrated by placing known concentrations 

of CO in the sample cell having the transmission t (uj) and using a calibration source 

c 

with a radiance, N (T,uj), 


AV cal = AQ f N ° (T,UJ) V W) t g ( u>) - R >) T A 1 du) 

Auj 

( 5 - 9 ) 

,R o / N ° (T ,w) t ( W ) [t («j) - t # ] duj 
Auj U A 

Thus 

f E (^) fR(UJ) T (U)) - R' (UJ) T ] dU) 

H o « f - 

f TO !> (uj) - T \ 1 duJ 
Auj u a 


f N ° (TO) T „ (W) rR(< 0 ) T (UJ) - R>) T ] duj 

Auj— - A 

f N ° TO T _ («>) Cr («) -t' ] duj 

Auj u a 


( 5 - 10 ) 


Computer calculations have been made in which the instrument is calibrated 
to obtain its overall responsivity and this responsivity applied to a signal from a 
real atmosphere. These results are then compared with the exact results that 
would be obtained if an overall responsivity were not used. 

Before these results are shown, a brief discussion of the computer simulation 
is pertinent. From the equations presented, it is clear that if a laboratory 
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calibration is made for varying optical thicknesses of CO, as u (CO) -*o 
d(DV)l - f o o 

=R 0 7 N < T S .»)*>• 

Ju=0 Aw 

and, knowing the blackbody calibration source temperature, the overall responsivity 
R o may be determined - However, for the computer simulation, one of the problems 
is in selecting the "correct" value of Atu over which to perform the calculations. Of 
course, in calibrating the real instrument this is automatically taken care of. 

For the computer simulation, one choice is to perform the calculations over a 
spectral band-pass defined by the cut-off frequencies. Another choice is to select 
a band-pass defined by the values when the transmission through the two instrument 
channels is equal to one-half. Both of these choices have been used in the computer 
simulations, as will be shown, and the results obtained by using these two choices 

should be representative of the values that would result from a real calibration of the 
instrument. 

The special case of 

R (ip) = R' ((d) = C (u>) * C' ((D) 

was investigated initially; i.e. , Run B-09 shown in Figure 5-3. The results of Run 

B-09 are re-presented in Figure 5-5. An identical computation except C(co) = 1.0 

over Aid is also plotted, Run B-18. These results are slightly different from those 

presented previously. The DV offset for (CO) =0 caused by the difference 

atm 

between the ground and instrument reference temperature has been corrected by DV 
(corrected) = DV - DV for <CO) >tffi -o. An effective R q (or C) le computed for 
( c O) at m =0, 1 P pm b y 

- = DV jRun B-09) 
o DV (Run B-19) 0,516 

The DV for Run B-18 is multiplied by R q (= 0. 516) and plotted as (X) on Figure 5-5. 

The results indicate that using the mean responsivity, R q , creates an error of less 

than one per cent in determining the concentration of CO in the atmosphere up to a 
level of 0.3 ppm. 
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Figure 5-5. Sensitivity of instrument as a function of CO concentration in the troposphere 
for two different slit functions. The lower curve is for the slit function given 
in Fig. 5-1, the upper curve for a square slit function. If the upper curve is 
corrected to account for the larger area under the slit function, the points 
marked as (x) are obtained. 
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A calibration of the instrument has also been simulated with the identical 

instrument parameters used in Run B-09, but a calibration source temperature of 

313. 2° K and a sample gas temperature of 288. 2°K were inputted. The results of 

these calculations are plotted in Figure 5r6. Run B-16 assumed the same instrument 

response function as Run B-09 and Run B-17 assumes C(w) =1. o over Aw. This mean 

responsivity, R , obtained by instrument calibration is R = 0.500: this differs 
° o 

from the atmospheric R q by less than 2 per cent. Changing the calibration source 

temperature to 290° K and the sample gas temperature to 260° K resulted in a value 

of R q = 0.512, which differs from the atmospheric R by less than one per cent; 

these temperatures more closely approximate the ground temperature and effective 

atmospheric temperature used in Run B-09. A series of calculations were made 

for the special case of C(w) = c'(<v) and the results are summarized in Table 5-2. 

Runs B-09, B-18 and B-43-2 give the values for DV computed through the real 

c 

atmosphere. The other results are for a one-slab calibration. The spectral response 
function for both channels, C(w) = C^w), is shown in Figure 5-7. 

A s indicated previously, the overall responsivity (indicated by C in Table 5-2) 
is obtained by a computer simulation of a one-slab calibration. For example, DV 

c 

for Rim B-16 was computed using a blackbody calibration source temperature, 

T° = 313. 2° K, and a temperature of the gas in the sample calibration cell, T = 
b sc 

288. 2° K, and the spectral response functions, C(to) = C '(id); dv for Rim B-17 was 

c 

computed using identical conditions except C (to) = C / (to) =1 from tu = 2070 cm” 1 
-1 1 
to a > 2 - 2220 cm . is computed as a function of CO concentration by 

DV (B-16) 

C = - 

1 DV c (B-17) * 

which results in ^ 0.5. Defining the spectral bandpass by the values when the 
transmission is one-half, i.e. , w = 2 111 cm" 1 and to = 2176 cm" 1 , gives values 

X £ 

for C g «0.98. Multiplying the responsivity, C, with the results of the atmospheric 
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Table 5-2. Effects of Lab Calibration Applied to Real Atmosphere 


_2 

DV w/cm -sr ; t = .8178 ; C(u)) = c'(tu) ; <# = 2143 cm' 

* A O 


< CO >atm 

PPM 

B-09 

10 7 dv 

C 

B-18 
10 7 dv 

C 

B-43-2 

7 

10 DV 

c 

B-16 

10 7 dv 

C 

B-17 

7 

10 DV 

C 

B-43-1 
10 ‘dv 

c 

C 1 

C 2 

B -32-1 

10 7 dv 

c 

B-32-2 

10 7 dv 

c 

C 3 

.500 x 

DV(B-18) 

7 

X 10 

.981 x 
DV(B-43-2) 

xlO 7 

.512 

DV(B-18) 

7 

xio 

0 

0 

o 

0 

0 

0 

0 

— 

— 

0 

0 

— 

0 

0 

0 

.025 

1.746 

3.32 

1.806 

3.561 

7.08 

3.617 

.5030 

.9845 

2.242 

4.32 

.5191 

1.660 

1.772 

1.700 

.050 

2.915 

5.58 

3.006 

5.802 

11.58 

5.903 

.5026 

.9829 

3.599 

6.97 

.5164 

2.790 

2.949 

2.857 

.075 

3.772 

7.25 

3.897 

7.370 

14.75 

7.505 

.4997 

.9820 

4.527 

8.81 

.5138 

3.625 

3.823 

3.712 

.100 

4.442 

8.57 

4.581 

8.556 

17.13 

8.718 

.4995 

.9814 

5.220 

10.20 

.5118 

4.285 

4.497 

4.388 

.125 

4.991 

9.65 

5.140 

9.505 

19.05 

9.6G9 

.4989 

.9810 

5.770 

11.30 

.5106 

4.825 

5.042 

4.941 

.150 

5.455 

10.57 

5.616 

10.294 

20.66 

10.502 

.4983 

.9802 

6.224 

12.22 

.5093 

5.285 

5.509 

5.412 

.200 

6.212 

12.09 

6.388 

11.555 

23.21 

11.796 

.4978 

.9796 

6.943 

13.67 

.5079 

6.045 

6.267 

6.190 

.300 

7.318 

14.31 

7.524 

13. 337 

26.81 

13.620 

.4975 

.9792 

7.941 

15.71 

.5055 

7.155 

7.381 

7.327 

.400 

8.118 

15.92 

8. 331 

14.571 

29.30 

14.884 

.4973 

.9790 

8.616 

17.08 

.5045 

7.960 

8.173 

8.151 

C<«) 

150 cm 


n 

75cm~* 

150 cm 


n 

75 cm - * 



150 cm -1 



3.5% 

- 1.2% 

1.2% 


T ° = 288. 

u 

2°K 


T ° = 313. 2°K 
o 







DV Error 

at 0. 1 ppm 



T ~ 2G0‘ 
atm 

K 


T = 288. 2°K 
EC 



i 








Atmosnhere 

Calibration 



Calibration 


































2] 


Wavenumber, 




computations, which were made using corresponding rectangular spectral response 
functions, simulates the application of the calibrated response to the case of a real 
atmosphere. The values thus obtained may be compared with the "true" values 
that would result if an overall responsivity had not been used. For example, at 

C CO = 0,1 ppm C i °' 5) X DV C for Run B “ 18 (= 8.57 x 10 -7 ) gives 4.285 x io" 7 
which compares with DV c for Run B-09 (= 4.442 x io" 7 ); or an error of 3.5% in DV 

c 

would result from applying the calibrated instrument responsivity. 

Also presented in Table 5-2 are the results obtained for the calibrated mean 
instrument responsivity when the calibration is simulated using temperatures more 
closely approximating the real atmospheric temperatures. These results indicate 
errors in DV c will be less by about a factor of three. The results presented in 
Table 5-2 indicate better agreement between "true" values for DV and "one-slab 
calibration" values if the half-width of the spectral response function is selected for 
the computer simulation instead of the cutoff frequencies. 

Computer simulations have also been performed for cases when C(cu) ^ C / ((ju). 
Two extreme cases were selected. In both cases C(u>) was maintained as indicated 
by Figure 5-6; but, c'(u)) was shifted with respect to C(u>) by plus 20 cm" 1 in one case 

and by minus 20 cm in a second case. The results are presented in Tables 5-3 
and 5-4. 

As in the case where C(<u) = C '(uj), good agreement is obtained between the 
"true" values for DV c and the valued obtained when a calibrated mean instrument 
responsivity is applied. The best agreement results when the integrations are 
performed over the spectral intervals defined by the one-half transmission limits, 
namely, for Table 5-3, when = 2091 cm" 1 and o> 2 = 2176 cm” 1 and for Table 5-4 
when a) = 2111 cm 1 and cu = 2196 cm" 1 . 

To summarize, the computer calculations indicate that a one-slab calibration 
of the instrument’s responsivity may be applied to signals obtained by the instrument 
when observing a real atmosphere. This means that a detailed knowledge of the 
spectral responses of the two instrument channels is not necessary. 
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Table 5-3. Effects of Lab Calibration Applied to Real Atmosphere 


2 i 

DV , W/em - sr ; t = .7642 ; C(u>)^c'(w); u) (c) = 2143 cm ; <o (C ') = 2123 cm' 

C A 0 0 


(CO) atm 
PPM 

B-10 
10 7 dv 

C 

B-25 

10 7 dv 

c 

B-45-2 

10 7 dv 

V 

B-21 

10 7 dv 

c 

B-22 

io 7 dv 

c 

B-45-1 

10 7 dv 

c 

6 i 

S 2 


.771 x 
DV(B-45-2) 
7 

XI 0 

0 

0 

0 

0 



0 

— 

— 

0 

0 

.025 

1.592 • 

3.13 

2.111 

3.266 

6.67 

4.28 

.4896 

.7631 

1.518 

1.628 

.050 

2.663 

5.25 

3.503 

5.322 

10.92 

6.94 

.4874 

.7669 

2. 546 

2.702 

.075 

3.446 

6.83 

4.523 

6.750 

13.90 

8.78 

.4856 

.7688 

3.312 

3.488 

.100 

4.056 

8.08 

5.314 

7.820 

16.14 

10.14 

.4845 

.7712 

3.919 

4.098 

.125 

4.553 

9.10 

5.S55 

8.669 

17.96 

11.22 

.4827 

.7726 

4.414 

4.592 

.150 

4.972 

9.96 

6.478 

9.368 

19.47 

12.11 

.4812 

.7736 

4.831 

4.996 

.200 

5.651 

11.39 

7.349 

10.474 

21.81 

13.50 

.4802 

.7759 

5.524 

5.668 

.300 

6.631 

13.45 

8.598 

12.008 

25.12 

15.43 

.4780 

.7782 

6.523 

6.631 

.400 

7.330 

14.94 

9.474 

13.043 

27.33 

16.73 

.4772 

.7796 

7.246 

7.306 






| 






CjWCCu) 


M 

L 


am 

L i 



3.4% 

-1.0% 










T° = 288 
G 

2“K 


T °= 313. 
s 

2° K 







T ~ 260° K 

atm 


T = 288. 
sc 

2° K 







Atmosphere 

Calibration 
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Table 5-4. Effects of Lab Calibration Applied to Real Atmosphere 


DV , w/c 
c 

2 

;m -sr ; 

t = .8753 
A 

; C(u>) A C>) ; 

d (C) = 2143 cm" 1 

; ui (C') = 2163 cm" 1 
o 


< CO > atl r 

PPM 

B-ll 

10 7 dv 

C 

B-26 

10 7 dv 

c 

B —45 “4 
7 

10 DV 

c 

B-23 

10 7 dv 

c 

B-24 

10 7 dv 

c 

B-45-3 

10 7 dv 

c 

V 

ol 

to 

.480 x 
DV(B-26) 
7 

X 10 

.705 x 
DV(B-45-4) 
7 

x 10 

0 

0 

0 

0 

0 

0 

0 

— 

— 

0 

0 

.025 

1.802 

3.60 

2.627 

3.773 

7.69 

5.38 

.4854 

.7013 

1.727 

1.852 

.050 

2.997 

6.05 

4.343 

6.124 

12.62 

8.71 

.4852 

.7031 

2.903 

3.062 

.075 

3.870 

7.88 

5.605 

7.765 

16.11 

11.03 

.4820 

.7040 

3.781 

3.952 

.100 

4.553 

9.33 

6.588 

9.010 

18.78 

12.78 

.4798 

.7050 

4.477 

4.645 

.125 

5.112 

10.53 

7.395 

10.010 

20.93 

14.17 

.4782 

.7064 

5.052 

5.213 

.150 

5.586 

11.56 

8.072 

10.843 

22.74 

15.36 

.4748 

.7059 

5.546 

5.691 

.200 

6.360 

13.25 

9.200 

12.181 

25.67 

17.28 

.4745 

.7050 

6.357 

6.486 

.300 

7.491 

15.79 

10.857 

14.082 

29.88 

20.00 

.4713 

.7041 

7.576 

7.654 

.400 

8.309 

17.66 

12.084 

15.405 

32.88 

21.97 

.4685 

.7012 

8.473 

8.519 



J L 

n 


J l 

n 





C(®),C>) 


180 cm * 

85 cm 


180 cm * 

85cm 



1.7% 

2.0 % 


T ° = 288. 2” K 

VJ 


T° = 313. 

2° K 




DV Erro 

r @ 0. 1 ppm 


T ~ 260° K 
atm 


T sc = 288. 

2° K 







Atmosphere 

Calibration 








5.3 IN-FLIGHT CALIBRATION 

In this section we show the results of the investigation about the influence of 

different calibration temperatures on AV. The ground temperature was kept at 

270° K, but in one case T = 270° K, in another case T = 300° K. The results for 

u C 

four different CO concentrations are as follows: 

Table 5-5. Signal Change of CO for the Calibration Temperatures 


T c 

300® K 

270°K 

Difference 

e GR 

1 

1 


AV @ 0 ppm 

1.363 X 10" 9 

0 

1.363 X 10" 9 

.05 

6.7216 X lo" 8 

-8 

6.5865 X 10 

1.351 X 10~ 9 

.1 

1.1515 X 10“ 7 

1.1381 X 10“ 7 

1.34 X lo" 9 

.2 

-7 

1.7772 X 10 

1.7639 X 10" 7 

1.33 X 10" 9 

.3 

2.2054 X lo" 7 

2.1922 X lo" 7 

1.32 X 10~ 9 

.4 

-7 

2.5390 X 10 

-7 

2.5259 X 10 

1. 31 X lo” 9 


As expected, in the case for which the ground and calibration temperatures 
are the same, AV = 0 for c qq = 0. But for a temperature difference of 30°K, 

AV° = 1. 363 X 10 W/cm sr. This offset is only about 1.4% of the signal for 

c co = 0-1 ppm. It is further found that the offset remains nearly constant for all 
CO concentrations, as shown in the last column. 
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The offset AV° as a function of ground temperature, emissivity and sun angles 
was further investigated. In Figure 5-8 we have plotted A V® for the three temperatures 
270, 300 and 330°K, whereby T Q = 300°K. The different lines are for the sun 
zenith angles 0, 30, 60 and 90°. One can see that for low angles, AV* increases due 
to the fact that the wavelength dependency of the reflected sun energy is that of a 
6000°K blackbody and thus quite different from the 300°K calibration blackbody. The 
offset AV° also increases with decreasing e , since then the reflectivity increases. 

In all cases, maximum AV° i s obtained for G = 0. The exception is the night time 
(9-90 ) and AV°- 0 for e = 0. Since for e = 0 only sun reflected energy is observed, 
the lines for a given 9 , but different T qr , must have the same value for e = 0. 

These points are convenient check values for the computer program as well as the 
hand plots. 

For the higher ground temperatures, the absolute magnitude of AV° appears 
to be large. However, when AV° is related to the signal obtained for 0.1 ppm CO, 
it is found that for the present case AV° is never more than 3% of AV . These 
results are shown in Figure 5-9, where =AV® x 100 /AV q l is plotted versus the 
ground temperature for e = 0. 8 and 0. 4 and 9 as a parameter. 

5.4 INVESTIGATION OF RATIO TECHNIQUES 
5 - 4,1 Patio of Signal Difference of One Cell Pair to Signal 
It was found that the ratio DV/V reduces atmospheric effects on the measure- 
ments of CO. As an example, this ratio has been plotted for different concentrations 
of CO and water in the atmosphere (see Figure 5-10). The four upper curves are the 
signal differences as a function of CO concentration for a dry and three different wet 
atmospheres. If the water concentration of the wet atmosphere is uncertain to 
within 1/2 and 2 times the normal Gutnick model, the determination of the CO 
concentration is uncertain to within 42%. The lower four curves are the values of 
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Figure 5-8. Sensitivity of the instrument to different ground temperatures 

and sun zenith angles as a function of ground emissivity at 4.6 M- 
for zero concentration of CO in the troposphere. 
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330° K 


Figure 5-9. Sensitivity of the instrument for the same conditions as in 

Fig. 5-8, but referenced to DV for 0.1 ppm CO concentration. 
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DV/V for the same atmospheres. As can be seen, the uncertainties in the 
determination of the CO concentrations are reduced to about 14%. The influence of 
varying ground temperature is also reduced when the ratio is used instead of AV 
alone. 

Figure 5-11 shows the ratio versus c for T from 315 to 260° . At c = l 

CO G CO 

the difference in the ratio for T Q = 300° and 285° is 5 x io out of 30. 3 x io , or 

16.7%, whereas the difference in AV for T = 300 and 285° is 3 x lo" 7 W/ cm 2 s ter 
-7 ^ 

out of 6. 13 x io , or ~ 50%, a considerably larger difference. The value of the ratio 
will decrease as T Q approaches the average atmospheric temperature. This is seen 
in Figure 5-12, a plot of AV/V versus T Q for different CO concentrations. The 
sensitivity of AV/V for changes in T and water vapor are indicated in Figure 5-ll 
by line lengths. 

From these calculations the changes in AV/V were converted to changes in CO 

concentration. Interpreting these as independent error sensitivities, and assuming 

the expected uncertainties in T^, T, and water vapor to be ± 2K, ± 2K, and 20% 

respectively, the expected errors in CO concentration at c = . l ppm are 

CO 

± 3.3% for A Tq = ± 2K 
± 4.4% for AT = 2K 
± 2.4% for A WV= ± 20% 

± 6.02% for all three error sources. 

The relationship between AV and V was further examined for different 
conditions of CO concentration, surface temperatures, atmospheric temperature 
and water vapor profile perturbations, and for different surface emissivities 
and sun zenith angles. Figure 5-13 shows the plot of Av and V for the temperature 
and water vapor. Several interesting observations can be made. For instance, the 


1?0 
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Figure 5-11. DV/V as a function of CO concentration in the troposphere 
for different ground temperatures and variations of 
atmospheric temperatures and water vapor concentration. 
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Figure 5-13. DV as a function of V for different conditions of CO concentration, surface 
temperature, atmospheric temperature and water profile perturbations, 
surface emissivities, and sun zenith angles. 



quantity Q = AV/ ( V-42), which is the approximate slope of 4V vs. V at different CO 
concentrations, appears to be independent of surface temperature. It is likely that 
the term with value 42 in the indicated ease depends in value on the average atmos- 
phere emission, and could be expected to differ for different atmospheres. Removal 
of the surface temperature dependence from the AV/v relation could result also in 
less sensitivity of the CO concentration to other uncertainties, beside the surface 
temperature dependence. It might even be possible to combine the gas cell radiance 
values from several instrument channels in order to better estimate the observed 
surface temperature. Figure 5-13 also shows the effect of different values of surface 
emissivities. The skpes of AV versus V as the sun zenith angle is changed appear 
to be the same as the: slope for the standard case implying independence. A trans- 
lation of the curves with surface emissivity, however, is in evidence. It may be 
possible that different term values in the denominator for Q could be found to reduce 
the effect of surface emissivity uncertainties on the Q function. 

Figure 5-14 shows superimposed on the standard cases the curves of AV versus 
V for CO of 0.4 and 0. 1 ppm, the perturbation on AV and V due to changes in the 
atmospheric temperature profile by ± 5K and in the water vapor profile by ± 50% in 
mixing ratio. The perturbations are shown for ground temperature of 300 and 285K 
The indicators at 300K and 0. 4 ppm CO equally hold for the other cases. 

The uncertainties of Q due to water vapor changes appear to be largely 

removed, because they line up with the AV vs. V slope. Changes in atmospheric 

temperature affect the Q values for different amounts of CO. The uncertainties in 

temperature are expected to be small, therefore Q will like, vise have small 
uncertainties. 

This approach has to be further explored, especially with respect to 

uncertainties in the atmospheric scattered cuivn^t • , . 

P vne scattered sunlight, important at the shorter wave- 


124 




Figure 5-14. Perturbation of DV and V due to changes in the atmospheric temperature 
profile by ± 5“K and in the water vapor profile by ± 50% in mixing ratio. 


lengths, and the effects due to the atmospheric emission. This may lead to a viable 
method of reducing the uncertainties due to some of the parameters affecting the 
instrument performance. 

5.4.2 Ratio of Signal Difference of Two-Cell Pairs 

Calculations were made for different cell conditions, in which the total pressure 
was changed from zero to 1 atm, keeping the CO concentration constant. The signal 
difference A V as a function of total cell pressure is plotted parametrically with 27 
different CO concentrations in Figures 5-4, 5-15 and 5-16. The distribution of CO 
in the 8 layers is given in Table 5-6. The numbers on the curves refer to the factor 
with which every layer is multiplied. 


Table 5-6. CO Concentration for 3 Distributions 


Layer 

B20 


B 37 

B36 

1 

!0- 7 


io" 7 

io' 5 

2 

io- 7 


io- 7 

io' 6 

3 

io -7 


io -7 

io -7 

4 

io- 7 


io“ 7 

io' 7 

5 

io -7 


io -7 

io -7 

6 

l(f 7 


io- 6 

io -7 

7 

Hf 7 


io- 6 

io' 7 

8 

io~ 7 


io~ 7 

io- 7 

Total Path Length 

56.6 X io” 3 


174.9 x io “ 3 

1065.9 X 10 


cm atm 


cm atm 


-3 

cm atm 
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Figure 5-15. Sensitivity of instrument as a function of total pressure in the instrument 
cell for different concentrations of CO in the troposphere, having a profile 
according to B37 (see Table 5-6). 



The results of DV versus the optical path (cm-atm) are plotted in Figure 5-17 
for the three distributions and two cell pressures. The shapes of the curves and 
ratios appear to be sufficiently different to give some information about the distri- 
bution. This is shown in Figure 5-18 , where the ratio of the DV for the channels at 
different pressures versus optical path is plotted. The two distributions B-37 and 
B-20 are very similar, butB-36 (build-up of CO near the surface) shows a different 
behavior. 
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Figure 5-17. Sensitivity of the instrument as a function of the total optical thickness of 

CO for the three profiles B20, B37 and B36, using cell pressures of 1 and 0.4 atm. 
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Figure 5-18. The ratio of DV as a function of total optical thickness for three CO con- 
centration profiles (see Fig. 5-17), using two cell pairs with pressures of 
1 and 0.4 atm. 



STUDY OF MISSION REQUIREMENTS 
(Task 5) 


In this section, the supporting data requirements are investigated. These are 
atmospheric temperature and pressure profile, the water vapor profile, and the 
extent of cloud cover and cloud top height. The effect of the uncertainties in these 
atmospheric parameters on the signal difference of the RGFC instrument are 
investigated for the measurement of CO. 

6.1 ATMOSPHERIC TEMPERATURE AND PRESSURE PROFILE 
6 * 1 * 1 Temporal and Spatial Variations of Temperature Profiles in the Real Atmosphere 
In order to gain some insight into the spatial and temporal variations of the vertical 
atmospheric temperature profiles, meteorological data for two different climatological 
regions were analyzed. Archival data from the National Climatic Center in Asheville, 

N. C. of the Northern Hemisphere and U. S. for July 1970 and February 1971 were 
utilized. 

(a) Temporal Changes 

The temperature and spatial temperature gradients at the 500 mb level were 
plotted from 12 hours data for a month’s interval. Figure 6-1 shows this data for 
San Diego, February 1971; Figure 6-2 is for Washington, D.C. , February 1971; and 
Figure 6-3 is for San Diego, July 1970. Temperature data was from station reports 
for Figure 6-1 and from the upper level charts for Figures 6-2 and 6-3. The 
temperature gradient was obtained from the upper level chart by estimating the 
distance between isotherms in the vicinity of the station and then converting to degrees 
centigrade per 100 nautical miles. It should be noted that the isotherms on the upper 
level charts are spatially smoothed data and therefore, do not represent the exact 
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Figure 6-2. Temperature and temperature gradient at 500 mb level for 
February 1971 in Washington, D.C. at 0000 and 1200 Z. 
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Figure 6-3. Temperature and temperature gradient at 500 mb level for 
July 1970 in San Diego, California at 0000 and 1200 Z. 
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values that may have existed at the time. The temporal variations of the spatial 
temperature gradient may give a clue to the error in the spatial temperature gradient 
estimate. Several observations were made from this limited data: 

1. From the data no apparent diurnal cycles were in evidence - at least not with 
significantly large values. Some diurnal variation is seen in the monthly 
averages. For 12 hours data there are variations from other causes masking 
the diurnal variation. The station data for San Diego show a 0. 3°C difference 
between 00 hour and 1200 hour monthly averages. 

2. The 12 hour sampling interval is inadequate to show variations with finer 
time resolution. This could lead to errors due to under-sampling. Something 
to watch out for. 

3. Changes in temperature and spatial temperature gradient appear to be 
associated with the passage of fronts and cyclones in a somewhat systematic 
manner, and with thermal activity and errors in a more random manner. For 
instance, February in San Diego is predominantly under the influence of the 
Pacific High with stable conditions, both in temperature and spatial gradient. 
Changes in temperature occur near times when lows are moving in with cold 
air. Spatial temperature gradients are higher also when lows or troughs are 
approaching. During July in San Diego the situation appears to be more stable with 
fewer large temperature spikes, but with increased fluctuations which may be 

due to the strong thermal effects in the Southwest, Data for Washington D. C. , in 
February show fluctuating temperature between warm and cold as lows are 
approaching. The spatial temperature gradients are fairly well behaved, 
except near the first of the month, when a very deep low existed. 

1. It appears from these data, assuming that the upper chart data represent 
a low error from the actual situation, and that no temporal under-sampling 
error exists, that the spatial temperature gradient is stable and low 
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enough to allow interpolation with little error over relatively short distances (100 
n.mi.) most of the time. The monthly average of spatial temperature gradient at 
Washington in February was 1. 67° C per 100 n.mi. , with maxima of 3-10° C 
per 100 n.mi occurring near lows. 

5. The spatial temperature gradients appear greatest when lows are approaching. 
These events are relatively rare events, considering all time and space, 
and occur near storms when GFC instrument data may also be less efficient 
for other environmental reasons such as clouds. 

(b) Spatial Changes 

The spatial distribution of temperature over the USA is shown in Figure 6-4 
for February 1, 1971. The (5 degree) isotherms are widely spaced (for low spatial 
temperature gradients) most everywhere. They are very closely spaced near the 
large low pressure area over eastern Canada implying large spatial temperature 
gradients. To illustrate the temperature gradients and changes in gradients, the 
upper level chart data was used to plot temperature against degrees latitude for 
every 5 degrees in longitude from 90° W to 120° W. These plots are shown in Figure 
6-5. 

1. From these data it appears quite feasible to interpolate over relatively long 
distances, even for large spatial temperature gradient areas using methods 
similar to those used operationally by NMC. Changes in temperature gradient 
are generally small with larger changes occurring more rarely at air mass 
boundaries. The larger values of finite differences lead to interpolation in 
errors, especially for larger spatial intervals. 

2. For February 1, 1971 the largest spatial temperature difference is about 
9. 82° C per 100 n.mi. with an average of about 2°C per 100 n.mi. The 
maximum value of the second temperature difference is estimated at about 
5° per (100 n.mi.) with an average value of zero. 
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Figure 6-5. Temperature spatial distribution for February 1, 1971. 




3. The error can be estimated for linearly interpolated data with finite values for 
second differences and zero higher order differences. This error is given by 
the values of the truncated Gregory-Newton interpolation formula. The error 
varies within the interval and, hence, the average error is obtained by inte- 
gration of the error over the interval. This leads to an interpolation error of 
A T/12, where A T is the second order finite difference. The maximum average • 
error for a maximum a 2 T of 5° per (100 n.mi.) 2 is then 5/l2°C over 100 n.mi. 
intervals. This is small compared to the expected value of 2°C measurement 
error. Smaller errors can be expected for the interpolation method used by 
NMC, because it is quadratic in two dimensions. 

6-1.2 The Accuracy of Temperature Data 

Existing verification procedures compare objective analysis of temperature to actual 
observations. Verification statistics grouped by geographical regions and by altitude 
were assembled for us by NMC. It appears that the existing meterological observation 
system produces objective analyses and forecasts which are sufficiently accurate for the 
measurements of trace contaminants. Current operational inclusion of data from satellite 
atmospheric soundings has improved the objective analysis and forecast in conventional 
data sparse areas (over oceans) as well. This has been verified experimentally with 
SIRS data. Temperature and water vapor profile soundings will be available from 
satellites, such as the NOAA series (ITOS-D 1972, ITOS-E 1973, ITOS-F 1974, ITOS- 
N 1975), the NIMBUS series (NIMBUS V 1972, NIMBUS VI 1974) and on geosynchronous 
satellites GOES (SMS-B 1973), and ATS-G(1974). The low altitude satellite sound- 
ings will take place at asynoptic times so that forecasting techniques must be used to 
incorporate these into the objective analysis and forecasts. The profiles from 

synchronous altitude satellites will be capable of giving soundings every three hours 
over the entire visible disc. 
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Objective analysis produces a best estimate of the state of the atmosphere at 
synoptic times (every 12 hours) using all available observational sources, including 
ship report, aircraft and SIRS soundings. These analyses are produced twice daily 
by the NMC analysis-forecast operation. The objective analysis is performed using 
corrections to a first guess field derived from actual data entered twice daily. The 
corrections are determined from a comparison of the data with the interpolated 
values of the given field at the observation points. For analysis of pressure heights 
the reported wind is taken into account in determining the lateral gradient of the 
correction. Several scans are made with successively smaller interpolation scales^. 
Forecasts are produced from the objective analyses. The methods used in objective 
analysis and forecasting are suitable for temporal and spatial interpolation to areas 
of air pollution observation from satellites. 

(a) Radiosonde Accuracy 

Analysis of paired radiosonde flights showed the RMS temperature error of 
the temperature sensor to be between 0. 3°C and 0.4°C. This error and an assumed 

RMS pressure error of 3 mb combine to give RMS height and density errors 
at specified pressure levels. At 100 mb pressure level, the RMS height error was 
reported at 105 ft (32 m)^ \ 

(b) Temporal and Spatial Temperature Changes 

( 2 ) 

Project Stormy Spring used a mesoscale rawinsonde network with ascents 
every 1.5 hours during the passage of five storms in April 1965. 
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Empirical study of the data confirmed the implication of the thermal wind relation 
that the mean temperature in a deep layer will not change rapidly in space or time 
and that windshears can be used as predictors of the changes. The thermal wind 
relation implies for each one degree change in mean temperature per 100 km a wind 
change of 39 m/sec (78kt) across the 7-16 km layer (similar for other layers). The 
observed extreme geostropic wind of less than 200m/sec from 400 to 100 mb gives 
a horizontal temperature gradient of 5. 1° /100 km. The Stormy Spring data showed 
that the temperature at a given level is considerably more variable than the average 
temperature between pressure levels. Changes in layer mean temperature were 
computed for intervals of 1.5 to 12 hours and space changes of 100 to 350 km 
(4. 1° C for 6 hours and 3. 5° C for 300 km). Large temperature changes over short 
time and space intervals were not found. This is an important characteristic of the 
atmosphere, and much smaller changes could be expected for clear air, non-storm 
situations which are of most interest to air pollution measurement. 

(c) Relative Accuracy of SIRS and Radiosonde Data 

The relative accuracy and utility of SIRS data relative to radiosonde data for 
objective analysis was determined using tropospheric height and thickness values 
derived from NIMBUS 3 SIRS soundings (thickness of pressure layers is related to 
layer temperature) (4) . Relative accuracy and utility depends strongly on cloud 
condition, pressure level and instrument status. The SIRS thickness in the upper 
troposphere compares well to radiosonde thickness in the absence of high clouds. 

Both data gave thickness rms errors of between 20 to 30 meters with the SIRS thick- 
ness error slightly greater. 

(d) Temperature Errors from SIRS 

The standard error of temperature was found to be between 1» C and 2°C 
for pressure levels from 850 to 30 mb from NIMBUS 3 SIRS data from May-June 
and October-November 1969^. 
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Correction for cloudiness and high terrain is necessary in the derived 
regression equations. A comparison with radiosondes with known surface tempera- 
ture shows good agreement at all levels up to the tropopause. In the Southern 
hemisphere profiles are degraded by lack of surface temperature* and 
inversions aloft cannot be retrieved. 

A comparison of geopotential heights from SIRS and objective analysis show 
reasonable agreements. 

A comparison of 500 mb analysis with and without SIRS soundings gave deeper 
lows over mid-Pacific with a maximum difference in height of 120 meters. 

Twenty-four, forty-eight hour and 3-day forecasts were made from the two analyses 
and compared to later analyses. The forecast made from the analyses which incor- 
porated the SIRS data showed marked improvement. 

(e) SIRS Sounding for Operational Use 

Temperature profiles can be retrieved from satellite radiance measurements 
with an algorithm to find the solution which differs least from the initial guess 
in order to satisfy the observed radiance observations^. The initial guess is from the 
dynamic forecast. SIRS-B data ha*e been used for NMC analysis-forecast operations. 
Verification statistics indicate SIRS derived temperatures are, in most cases, better 
than the 12-hour forecast. Partial cloudiness is a major problem source. 

(f) Sea Surface Temperature 

Data from HRIR on ITOS1 were used to determine global sea surface tempera- 
ture to an accuracy of between 2-3K (3. 5-4. 1M-) (standard deviation) as compared 
(7) 

to ship reports . 

6.1.3 Summary of Requirements for Atmospheric Temperature and Pressure 
Profile Data — — — 

The measurements in the thermal region depend upon the temperature and pressure 

distribution along the optical paths within the field of view existing at the time 
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data is taken. Subsequent processing of the RGFC data for data interpretation will 
require an estimate of the average temperature and pressure profile with height 
within the field of view. It appears from an examination of methods and techniques 
used that the existing meteorological observation system will be able to provide the 
required estimate of the atmospheric state. 

The National Meteorological Center (NMC) routinely produces objective analyses 
as a best estimate of the state of the atmosphere from all available observations. 
Forecasts are produced from the objective analyses. The methods employed in 
objective analysis and forecasting can be used for the estimation of the atmospheric 
temperature and pressure profiles for the locations and times applicable to the RGFC 
data. The NMC verifies objective analyses and forecasts by comparing these products 
with observations. The resulting verification statistics give the average error between 
the temperature observations at different pressure levels and the smoothed best 
estimate of the state of the atmosphere which was produced from the observations. 
Table 6-1 gives the mean absolute temperature error at different pressure levels 
for four months of 1970. The errors are less than IK with no essential difference 
between the statistics obtained from the 70 radiosondes in the northern hemisphere 
or the 19 radiosondes in the USA. 

The operational analysis and forecast operation includes all available 
operationally useful data from satellite temperature sounders. Reported experiments 
have shown that SIRS data has significantly improved the objective analysis in data 
sparse areas (judged from forecast performance with and without inclusion of SIRS 
data). However, it appears that a SIRS sounding does not give an independent 
observation of temperature as a radiosonde would, but relies on the operational 
objective analysis and short term forecasts (up to 12 hours) to provide an initializa- 
tion for SIRS analysis. The temperature profile from the objective analysis is 
changed if the SIRS radiances are inconsistent with it, until the computed and 
measured radiances agree within the expected errors of the instrument. This can 
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Table 6-1. Verification of NMC Operational Temperature Analysis 

MEAN AB SOLUTE TEMPERATURE ERROR (K) 


Pressure Level 
(mb) 

1970 

Jan 
a b 

April 
a b 

July 
a b 

Oct 
a b 

200 

.8 .9 

• 

• 

00 

.7 .7 

.7 .7 

250 

.7 .7 

.6 .7 

.6 .5 

.6 .7 

300 

.7 .6 

.7 .7 

.5 .5 

.6 .6 

500 

.6 .7 

.6 .7 

.5 .5 

.5 .6 

700 

.7 .7 

.6 .6 

.5 .6 

.6 .7 

850 

.8 1.0 

.7 .9 

.6 .7 

.7 .9 


a 70 radiosonde stations in Northern Hemisphere 
b 19 radiosonde stations in U.S.A. 

be interpreted to mean that SIRS type data can be expected to improve the analysis in 
data sparse regions, but that the objective analysis will remain the best estimate of 
the atmospheric state. 

Calculations have been made using the COLAYER program to obtain an under- 
standing of the instrument sensitivity to changes in the atmospheric temperature 
profile. The atmosphere was approximated by eight layers with average temperatures 
and pressures in each layer. The instrument sensitivity was determined by making 
different computer runs with the temperature in each layer changed by +5K, -5K 
from the standard average temperature. These runs were made for different CO 
concentrations, for an assumed surface temperature of 300K and for an assumed 
water vapor distribution according to Gutnick. The s aisitivity to temperature 
changes can be interpreted in terms of an error in CO concentration for ass um ed 
uncertainty in atmospheric temperature. In Section 5, the expected relative error in 
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CO centration (0.1 ppm) was reported as ± 4.4% for temperature uncertainties of 
± 2K. 

It seems reasonable from the foregoing to draw the preliminary conclusion 
that sufficiently accurate atmospheric temperature profiles can be obtained from 
operational objective analyses. An average error (from verification statistics) of 
± IK could lead to an average error in estimated CO concentration of ± 2 . 2% at 
0.1 ppm. The instrument sensitivities and expected errors would, of course, be 
different at other values of ground temperature, atmospheric relative humidity, 
other atmospheric pressure distributions, and for other species and concentrations. 
The effect of these have to be further examined. Additionally more detail examina- 
tions of the uncertainties of atmospheric temperature determinations are warranted. 

The following subsections present more detail discussions of atmospheric 
temperature uncertainties, and objective analysis and verification methods. 

(a) Objective Analysis 

The following was extracted in part from documents received from the NMC. 
Objective analysis of meteorological parameters are produced twice a day. "The 
height, temperature, and wind analyses are made for the 1977-point octagonal grid 
for the Northern hemispheres for each of nine mandatory pressure levels from 850 
to 100 mb. On the same grid, analyses are also made of the following parameters: 
height and temperature of the 70-, 50-, 30-, and 10 mb levels; sea level pressure; 
surface air temperature; temperature of four selected isentropic (potential tempera- 
ture) surfaces (for tropopause modeling); and the mean relative humidity of the surface 
to 500-mb layer. The sea level pressure and surface air temperature analyses are 
used to compute the 1000 mb height field. These analyses serve as input to the 
primitive equation (PE) forecast model." 

Tropical analyses of wind and temperature at different pressure levels are 
also produced on a grid extending between 48S to 48N. Current plans are to extend 
the analyses and subsequent forecasts on a global basis. Work is in advanced stages 
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of testing and verification. "The analysis techniques essentially apply corrections to 
a first-guess field. These corrections are determined from a comparison of the data 
with the interpolated value of the guess field at the observation point. For the analysis 
of the height of a pressure surface, the reported wind is taken into account in 
determining the lateral gradient of the correction to be applied. 

The analysis is performed in a series of cycles, usually 2 to 4 in number. 

The first-guess value at each grid point is adjusted by examining the reports in the 
vicinity. Once each point has been adjusted on a given cycle the newly assigned values 
are used in a succeeding cycle. For the first two cycles, the number of reports in 
the vicinity of a grid point is used to determine how far away an observation may be 
and still exert some influence on a grid point. For dense data regions, like the 
United States, this distance is less than for a sparse data region, like the Atlantic 
Ocean. " 

(b) Verification of Temperatures 

The first paragraph is from NMC descriptions of verification procedures. 
"Seventy radiosonde stations have been chosen as the locations at which the veri- 
fication results are computed. Analyses and forecasts based on 1200 Z data are 
verified. Values of wind and temperature at 850, 700, 500, 300, 250 and 200 mb are 
interpolated biquadratic ally to each of the verification stations from the grid point 
values of the operational analyses at those six levels. An identical interpolation 
procedure is applied to the grid point values of the operational forecasts at initial 
time and at 12, 24, 36 and 48 hours after initial time. Both the 0000Z and 1200Z 
radiosonde reports are examined and the wind and temperature values at the verifi- 
cation levels are saved for each of the seventy stations. " 

The verification statistics are computed monthly. Similar errors could be 
computed on a daily basis, but is not currently available. The average error from 
19 verification stations from the U. S. for a given objective analysis can be expected 
to be somewhat more variable than the monthly average errors because the error is 
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synoptic, i.e. , depends on the presence of storm systems. The average errors from 
70 verification stations, because they are distributed over the Northern hemisphere, 
can be expected to be closer to the monthly a\erage errors. 

(c) Uncertainty in Atmospheric Temperature Measurement 

The temperature measurements with a radiosonde appear to be quite accurate. 
This was determined by R. W. Lenhard at AFCRL from analysis of data obtained 
from paired simultaneous radiosonde flights released about 10 miles apart. These 
paired ascents show differences in heights of pressure surface which implies that 
the radiosonde system delivers temperature information accurate to 0.3K (RMS). 

The radiosonde measures temperature and pressure. The average temperature 
between two pressure levels is related to the thickness of the layer between the 
pressure levels by the hydrostatic equation: 


f Td( * nP) ~7 Ti "ir 

D 1 

where z are geopotential height of the pressure levels P t and P^, respectively, 

T is the average virtual temperature, R is the dry gas constant for air, and g is 
acceleration of gravity. This gives the height differential between pressure levels 
and hence the temperature and pressure profiles with height. 

An additional property of the atmosphere is that the spatial gradient of the layer 
thickness determines the wind velocity within the layer. A radiosonde, when tracked 
by radar, can yield the instantaneous wind vectors. These are smoothed to determine 
the wind velocity at pressure levels. The thermal wind equation is 
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where W is the geos trophic wind vector in the layer with thickness Az and average 
temperature T, V^T is the horizontal gradient of the mean temperature in the layer, 
f is the coriolis parameter and k is a unit vector in the local vertical direction. 

Spatial gradients of temperature tend to be high near the boundary of airmasses. 
At these boundaries large changes of temperature and pressure occur and are 
accompanied by high wind velocities. Geographic regions with spatial temperature 
gradients which are different from pressure level to pressure level then must have 
spatial gradients of layer thickness. 

Project Stormy Spring used a mesoscale rawinsonde network with ascents 
every 1.5 hours. Data was taken during the passage of five storms in April 1965 by 
C. Kreitzberg at AFCRL. Empirical study of the data confirmed the implication of 
the thermal wind relation that the mean temperature in a deep layer will not change 
rapidly in space or time and that windshears can be used as predictors of the changes. 
The thermal wind relation implies for each one degree change in mean temperature 
per 100 km a wind change of 39 m/sec (78kt) across the 7-16 km layer (similar for 
other layers). The observed extreme geostrophic wind of less than 200 m/sec from 
400 to 100 mb gives a horizontal temperature gradient of 5.1°/l00 km. The Stormy 
Spring data and other data show that the temperature at a given level is considerably 
more variable than the average temperature between pressure levels. Large 
temperature changes over short time and space intervals were not found. This is an 
important characteristic of the atmosplere, and much smaller changes could be 
expected for clear air, non-storm situations which are of most interest to air 
pollution measurement. 
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The variability of temperature at different pressure levels was examined 
from hourly radiosonde data collected by AFCRL and available in GRD Research 
Note No. 60. 

Figure 6-6 shows the variability of the 700 mb temperature and the thickness 
of the 700-500 mb layer from hourly radiosondes. The effect due to 12 hour sound- 
ings is illustrated in the figure by the straight lines joining 12 hourly soundings. The 
mean error and standard deviation was evaluated fi'om the soundings as the deviation 
from the straight line. This was done for the two 12-hour segments and for the total. 
The relative temperature error ct^/T is 1.34/(273 ± 8.5) or 0.5%, the relative 
thickness error cr /Az is 6.03/2580 or 0.233%. 

This indicates that indeed there is less vai’iability in the average temperatui'e 
within alayei than the temperature at a level. In the example a ratio of two was 
found. 


6. 2 WATER VAPOR PROFILE 

6.2.1 Temporal and Spatial Variations of Water Vapor Profiles in the 
Real Atmosphere 

Atmospheric relative humidity profiles are not normally included in objective 
analysis, instead an average relative humidity within the 1000 to 500 mb layer is 
included. No verification statistics for this are available. Water vapor soundings 
from satellites, however, have been found to give relative humidity to within 20% 
of radiosonde relative humidity. 

The station data obtained from the National Climatic Center was examined in 

order to gain a feel for the variability of water vapor profiles with time and in 

\ 

comparison with the standard atmosphere. In Figure 6-7 the constant pressure data 
for three groups of three days during February 1971 for San Diego are shown. The 
temperature and the mixing ratio (MR) are shown for the 850 and 500 mb pressure 
levels. The mixing ratio was derived from the relative humidity, which was given 
in the station data using the relationship 
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Figure 6-6. Variability of relative humidity, temperature at 700 mb 
and thickness of 700-500 mb layer for 48 hours. 
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Figure 6-7. Variability of temperature and water vapor mining 

ratio for selected days in February 1971 in San Diego, 
California. 








MR (gm/kg) - * e s ( mmH g) * RH 

The value of e , the saturation vapor pressure, is tabulated as a function of 
s 

temperature. The dashed lines indicate the respective values for the standard 
atmosphere and Gutnick's model. 

The actual data compared with the standard values of temperature and mixing 
ratio show the effect of shifting air masses. From February 1 to 3 a weak low 
appeared (visible on the upper air charts and not included in this report) with a drop 
of temperature at the 850 mb level and a rise of water content. The 500 mb level 
became drier with no significant change in temperature. A different air mass 
existed between February 12 to 14. It was recognized that the Pacific High was 
dominating the situation with warm temperatures and low values of spatial tempera- 
ture gradients. Now in Figure 6-7 this influence of the Pacific High is seen, 
because both levels show warmer temperatures than both the atmospheric model and 
the situation on February 1 to 3. The water vapor content is also considerably lower 
than the Gutnick's model (but not anywhere near zero). During February 17 to 19 
a low pressure system developed offshore, and rapidly passed eastward. The 850 mb 
temperature was much reduced with a high water vapor content at this level. However 
the 500 mb temperature showed a 2-day transient with a marked decrease of mixing 
ratio after the passing of the low. 

It appears reasonable to speculate that the temperature and water vapor 
content at any given time and location depend primarily on the characteristics of 
the prevailing airmass. Since airmasses tend to be conservative, the temperature 
and water vapor profiles tend to be conservative until the airmass changes or moves 
on, or except for localized heating effects . This implies considerable stability both 
temporal and spatial, which will help the temperature and water content interpolation 
problem. 
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6.2.2 Water Vapor from SIRS (NIMBUS IV) 

Comparisons of satellite derived and radiosonde observed water vapor profiles 
indicate that the error of SIRS-derived relative humidity in the mid-troposphere 
(400-600 mb) is less than 20%. Relative humidity error in the lower troposphere 
(1000 to 600 mb) is less than 30%. This RAOB-SIRS comparison was for coinci- 
dent data (within 3° latitude and longitude and 6 hours) for 21 November 197(/ 8 \ 

6.3 CLOUD COVER 

A partial cloud cover produces a temperature drop sensed from HRIR which depends 
on the cloud altitude and percentage cloud cover. A method is derived to compute a 
cloud temperature difference correction, from per cent cloud cover and cloud 
altitude guesses. Theory was compared to data processed from NIMBUS 3 and 

BOMEX showing a rather good correlation (within IK) in the respective temperature 
differences^. 

The COLAYER program was modified so that a thin cloud layer of varying 
amounts of cloud cover could be inserted between each atmospheric layer. The 
cloud was assumed to have a spatially uniform emissivity . Although 
different cloud emissivities could have been used, only clouds with unity emissivity 
were used. The program computes the radiance at the top of each atmospheric 
layer and the transmission through each layer in the same way as before. However, 
at the assumed cloud layer the radiance is the average radiance of the cloud and the 
radiance from below the cloud, according to the cloud amount. Runs were made with 
fractional cloud amounts of 0, . 3, . 6 and 1. 0 at 1 km cloud altitude and at 1, 2, 3, 4, 5 
and 6 km for full cloud cover. 

It was found that the instrument response, AV, is linear with cloud amount, 
Figure 6-8. The values of AV for solid cloud cover depend on the cloud surface 
temperature and hence give increasingly smaller values of AV for higher clouds as 
shown in Figure 6-9. The ratio AV/V is less sensitive than AV to the surface 
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Figure 6-8. Signal difference vs. fractional cloud cover for different 
amounts of CO concentration. Cloud altitude = 1 km. 
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Cloud Altitude (Km) 


Figure 6-9. Signal difference vs. cloud altitude for 100% cloud 
cover for different amounts of CO concentration. 
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temperature. In Figure 6-10 the ratio AV/V plotted against CO concentration is seen 
to closely compare to the shape of the curves for ground temperatures of 300, 285 
and 270K. 


6.4 THEORETICAL SENSITIVITY OF GFC INSTRUMENT TO UNCERTAINTIES 
IN ATMOSPHERIC PARAMETERS 

Calculations were made with the COLAYER program to determine the changes in 
instrument signal output with changes in atmospheric temperature and water vapor 
concentration profiles. The optimum gas cell conditions and an eight layer atmos- 
phere (to 10 km) were used. The standard atmosphere temperature, pressure and 
water vapor profiles which were used for the calculations during the last month 
are shown in Figure 6-11. The temperature and pressure profiles are for the U. S. 
Standard Atmosphere 1966. The water vapor profile is Gutnick's model. 

The temperature profile was changed from the standard profile by +5 and -5 
degrees at all layers simultaneously for an atmosphere with four different 
water vapor concentrations: zero, standard, 50% and 150% standard. Further- 
more, different amounts of ground temperature were assumed. The instrument 
signal AV and the radiometer value of the vacuum channel RR6 (RR6 is the symbol 
for J* E(tu) C ' (u))t duu ) are plotted as a function of CO concentration in Figures 

J\ 

6-12 and 6-13 for no water vapor. The ±5° temperature profile changes are also 

plotted. At c = .1 ppm and T = 300° K, an increase of 5° everywhere gives a 

. 53 ^ 

decrease in signal AV of 100% (or 1. 33% per degree) for a corresponding 

• • I 7*1 

uncertainty in CO concentration of 12% (or 2. 6% per degree). The instrument 
radiometer output (RR6) in Figure 6-13 is a weak function of atmospheric tempera- 
ture. It is seen however that both AV and RR6 are strongly affected by the surface 
temperature. The effect of different amounts of water vapor on the instrument 
signal is shown in Figure 6-14. A 50% reduction in water vapor gives a 7% increase 


in signal at c co = 


. 1 for a corresponding uncertainty in CO concentration of 14%. 
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Figure 6-10. Ratios of DV/V vs. CO concentration for A (ground temp. 

- 300K), B (0.6 cloud at 1 km, T * 291K), C (ground 
temp. = 285, cloud at 1 km, T = 284. 9K), D (cloud at 
2 km, T = 278. 4K), E (cloud at 3 km, T = 271. 9K), 

F (ground temp = 270K, cloud at 4 km, T ■ 265. 4K). 
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Figure 6-13. Radiance vs. CO concentration for same conditions as 
Figure 6-12. 
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The signal uncertainties due to independent layer temperature uncertainties 
were determined for each atmospheric layer. These signal changes per temperature 
change are not statistically independent. Therefore, it was not possible to simply 
combine them to obtain meaningful estimates of signal uncertainty due to random 
distributions of layer temperature errors. Figure 6-16 shows the rate of change 
of signal with temperature for the dry and the standard water vapor model. The 
temperature uncertainties at low atmospheric layers appear weakened by the presence 
of water vapor. The altitude of greatest signal sensitivity to temperature errors 
appears to be about 3 km. However the shape of the curve shown in Figure 6-15 
could change for other temperatures and water vapor distributions and temperature 
errors in other layers. 

It was found similarly unprofitable to perturb the water vapor concentration 
at each atmospheric layer in order to be able to statistically combine errors and to 
obtain an estimate of the expected signal uncertainty. This approach was abandoned 
in favor of changing the entire profile, whether temperature or water vapor, so that 
the instrument signal uncertainties can be maximally bounded. It appears reasonable 
that the uncertainties determined in this manner are maximum with respect to 
random errors in temperature and water vapor profile and one might reasonably 
expect the actual expected signal uncertainties to be less for truly random errors. 
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Figure 6-15. Signal difference change per degree change per layer vs. 

altitude for no water and with water (Gutnick) for CO 
concentration of 0. 2 ppm and ground temp. = 300K. 
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